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Abstract 

We review the underpinnings of the standard Newton-Einstein theory of gravity, and identify 
where it could possibly go wrong. In particular, we discuss the logical independence from each 
other of the general covariance principle, the equivalence principle and the Einstein equations, and 
discuss how to constrain the matter energy-momentum tensor which serves as the source of gravity. 
We identify the a priori assumption of the validity of standard gravity on all distance scales as the 
root cause of the dark matter and dark energy problems, and discuss how the freedom currently 
present in gravitational theory can enable us to construct candidate alternatives to the standard 
theory in which the dark matter and dark energy problems could then be resolved. We identify 
three generic aspects of these alternate approaches: that it is a universal acceleration scale which 
determines when a luminous Newtonian expectation is to fail to fit data, that there is a global 
cosmological effect on local galactic motions which can replace galactic dark matter, and that to 
solve the cosmological constant problem it is not necessary to quench the cosmological constant 
itself, but only the amount by which it gravitates. 



■ 1 Introduction 

c$ | Following many years of research in cosmology and astrophysics, a picture of the universe has emerged 
[H El EI III E] which is as troubling as it is impressive. Specifically, a wide variety of data currently 
support the view that the matter content of the universe consists of two primary components, viz. dark 
matter and dark energy, with ordinary luminous matter being relegated to a decidedly minor role. The 
nature and composition of these dark matter and dark energy components is not at all well-understood, 
and while both present severe challenges to the standard theory, each presents a different kind of 
challenge to it. As regards dark matter, there is nothing in principle wrong with the existence of non- 
luminous material per se (indeed objects such as dead stars, brown dwarfs and massive neutrinos are 
well-established in nature). Rather, what is disturbing is the ad hoc, after the fact, way in which dark 
matter is actually introduced, with its presence only being inferred after known luminous astrophysical 
sources are found to fail to account for any given astrophysical observation. Dark matter thus seems to 
know where, and in what amount, it is to be needed, and to know when it is not in fact needed (dark 
matter has to avoid being abundant in the solar system in order to not impair the success of standard 
gravity in accounting for solar system observations using visible sources alone); and moreover, in the 
cases where it is needed, what it is actually made of (astrophysical sources (Machos) or new elementary 
particles (Wimps)) is as yet totally unknown and elusive. 
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Disturbing as the dark matter problem is, the dark energy problem is even more severe, and not 
simply because its composition and nature is as mysterious as that of dark matter. Rather, for dark 
energy there actually is a very good, quite clear-cut candidate, viz. a cosmological constant, and the 
problem here is that the value for the cosmological constant as anticipated from fundamental theory 
is orders of magnitude larger than the data can possibly permit. With dark matter then, we see 
that luminous sources alone underaccount for the data, while for dark energy, a cosmological constant 
overaccounts for the data. Thus, within the standard picture, arbitrary as their introduction might be, 
there is nonetheless room for dark matter candidates should they ultimately be found, but for dark 
energy there is a need not to find something which might only momentarily be missing, but rather to 
get rid of something which is definitely there. And indeed, if it does not prove possible to quench the 
cosmological constant by the requisite orders and orders of magnitude, one would have to conclude that 
the prevailing cosmological theory simply does not work. 

In arriving at the predicament that contemporary astrophysical and cosmological theory thus finds 
itself in, it is important to recognize that the entire case for the existence of dark matter and dark 
energy is based on just one thing alone, viz. on the validity on all distance scales of the standard 
Newton-Einstein gravitational theory as expressed through the Einstein equations of motion 



for the gravitational field g^. Specifically, the standard approach to cosmology and astrophysics is to 
take the Einstein equations of motion as given, and whenever the theory is found to encounter observa- 
tional difficulties on any particular distance scale, modifications are to then be made to T^ v through the 
introduction of new, essentially ad hoc, gravitational sources so that agreement with observation can 
then be restored. While a better understanding of dark energy and explicit observational detection of 
dark matter sources might eventually be achieved in the future, at the present time the only apparent 
way to avoid the dark matter and dark energy problems is to modify or generalize not the right-hand 
side of Eq. (0) but rather its left. In order to see how one might actually do this, it is thus necessary to 
carefully go over the entire package represented by Eq. (|T|). to determine whether some of its ingredients 
might not be as secure as others, and investigate whether the weaker ones could possibly be replaced. 
To do this we will thus need to make a critical appraisal of both of the two sides of Eq. ([!}. 

2 The underpinnings of the standard gravitational picture 

2.1 Massive test particle motion 

Following a year of remarkable achievement, achievement whose centennial is currently being celebrated, 
at the end of 1905 Einstein found himself in a somewhat paradoxical situation. While he had surmounted 
an enormous hurdle in developing special relativity, its very establishment created an even bigger hurdle 
for him. Specifically, with the development of special relativity Einstein resolved the conflict between 
the Lorentz invariance of the Maxwell equations and the Galilean invariance of Newtonian mechanics 
by realizing that it was Lorentz invariance which was the more basic of the two principles, and that it 
was Newtonian mechanics which therefore had to be modified. While special relativity thus ascribed 
primacy to Lorentz invariance so that observers moving with large or small uniform velocity could then 
all agree on the same physics, such observers still occupied a highly privileged position since uniform 
velocity observers form only a very small subset of all possible allowable observers, observers who could 
move with arbitrarily non-uniform velocity. Additionally, if it was special relativity that was to be the 
all-embracing principle, all interactions would then need to obey it, and yet what was at the time the 
accepted theory of gravity, viz. Newtonian gravity, in fact did not. It was the simultaneous resolution 
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of these two issues (viz. accelerating observers and the compatibility of gravity with relativity) at one 
and the same time through the spectacular development of general relativity which not only established 
the Einstein theory of gravity, but which left the impression that there was only one possible resolution 
to the two issues, viz. that based on Eq. To pinpoint what it is in the standard theory which leads 
us to the dark matter and dark energy problems, we thus need to unravel the standard gravitational 
package into what are in fact logically independent components, an exercise which is actually of value 
in and of itself regardless in fact of the dark matter and dark energy problems. 

In order to make such a dissection of the standard picture, we begin with a discussion of a standard, 
free, spinless, relativistic, classical-mechanical Newtonian particle of non-zero kinematic mass m moving 
in flat spacetime according to the special relativistic generalization of Newton's second law of motion 

m^ = , iW = , (2) 

where dr = (— yapd^d^) 1 / 2 is the proper time and r) af3 is the flat spacetime metric, and where we 
have indicated explicitly that the Riemann tensor is (for the moment) zero. As such, Eq. (J2J) will be 
left invariant under linear transformations of the coordinates but on making an arbitrary non-linear 
transformation to coordinates x^ and using the definitions 

A _dx x a 2 e d^oe_ 

we directly find (see e.g. jB], a reference whose notation we use throughout) that the invariant proper 
time is brought to the form dr = (—g^dx^dx") 1 / 2 , and that the equation of motion of Eq. (|2j) is 
rewritten as 

'D 2 x x \ (d 2 x x . dx^dx yS 

\ — m \ 1_ F A 

Dr 2 / \ dr 2 dr dr 



M^r\=™[ — +rlu-r--r- =° ' *W = > ( 4 ) 



with Eq. (j3J) serving to define D 2 x x / Dt 2 . As derived, Eq. (£Q) so far only holds in a strictly flat spacetime 
with zero Riemann curvature tensor, and indeed Eq. (J3J) is only a covariant rewriting of the special 
relativistic Newtonian second law of motion, i.e. it covariantly describes what an observer with a non- 
uniform velocity in flat spacetime sees, with the T x u term emerging as an inertial, coordinate-dependent 
force. The emergence of such a T x u term originates in the fact that even while the four- velocity dx x /dr 
transforms as a general contravariant vector, its ordinary derivative d 2 x x jdr 2 (which samples adjacent 
points and not merely the point where the four- velocity itself is calculated) does not, and it is only the 
four-acceleration D 2 x x /Dt 2 which transforms as a general contravariant four- vector, and it is thus only 
this particular four- vector on whose meaning all (accelerating and non-accelerating) observers can agree. 
The quantity is not itself a general coordinate tensor, and in flat spacetime one can eliminate it 
everywhere by working in Cartesian coordinates. Despite this privileged status for Cartesian coordinate 
systems, in general it is Eq. (0J) rather than Eq. (J2J) which should be used (even in flat spacetime) since 
this is the form of Newton's second law of motion which an accelerating flat spacetime observer sees. 

Now while all of the above remarks where developed purely for flat spacetime, Eq. (0J) has an 
immediate generalization to curved spacetime where it then takes the form 

(D 2 x x \ (d 2 x x ,dx»dx»\ n 

-(7^j--(^ + rJ— -j=o , v^o • (5) 

In curved spacetime, it is again only the quantity D 2 x x / Dt 2 which transforms as a general contravariant 
four-vector, and the Christoffel symbol T x u is again not a general coordinate tensor. Consequently at 
any given point P it can be made to vanish, 1 though no coordinate transformation can bring it to zero 

Under the transformation v {T x w )pi the primed coordinate Christoffel symbols (T'^)p will vanish at 

the point P, regardless in fact of how large the Riemann tensor in the neighborhood of the point P might actually be. 
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at every point in a spacetime whose Riemann tensor is non- vanishing. As an equation, Eq. (J5J) can also 
be obtained from an action principle, since it is the stationarity condition 5It/Sx\ = associated with 
functional variation of the test particle action 

It = —mc J dr (6) 

with respect to the coordinate x x . To appreciate the ubiquity of the appearance of the covariant 
acceleration D 2 x x /Dr 2 , we consider as action the curved space electromagnetic coupling 



r(2) 



mc J dr + e J dr— — A\ , (7) 



to find that its variation with respect to x leads to the curved space Lorentz force law 



t d x dx^ dx \ ,_, \ dx r .. . 



Similarly, the coupling of the test particle to the Ricci scalar via 

4?) = -mc fdr-K [ drR a a , (9) 



leads to 

(TO+ ^.,(£ +I *££) = _^(^ + ££) , (10) 

(Rfj, uaT necessarily non- zero here), while the coupling of the test particle to a scalar field S(x) via 

4 4 ) =-kf drS{x) , (11) 

leads to the curved space 

^ n fd 2 x x \ dx^ dx v \ _ ( XR dx x dx p \ n , n , . 

an expression incidentally which reduces to Eq. (jSJ) when S(x) is a spacetime constant (with the mass 
parameter then being given by mc = kS). In all of the above cases then it is the quantity D 2 x x /Dr 2 
which must appear, since in each such case the action which is varied is a general coordinate scalar. 



2.2 The equivalence principle 

Now as such, the analysis given above is a purely kinematic one which discusses only the propagation of 
test particles in curved backgrounds. This analysis makes no reference to gravity per se, and in particular 
makes no reference to Eq. ([TJ at all, though it does imply that in any curved spacetime in which the met- 
ric g^y is taken to be the gravitational field, covariant equations of motion involving four-accelerations 
would have to based strictly on D 2 x x /Dt 2 , with neither d 2 x x /dr 2 nor T x u (dx tl / dr){dx v / dr) having any 
coordinate independent significance or meaning. As such, we thus recognize the equivalence principle 
as the statement that the d 2 x x /dr 2 and T^i^dx^ / dr){dx u / dr) terms must appear in all of the above 
propagation equations in precisely the combination indicated, 2 and that via a coordinate transformation 
it is possible to remove the Christoffel symbols at any chosen point P, with it thereby being possible to 

2 This thereby secures the equality of the inertial and passive gravitational masses of material particles. 
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simulate the Christoffel symbol contribution to the gravitational field at such a point P by an acceler- 
ating coordinate system in flat spacetime. We introduce this particular formulation of the equivalence 
principle quite guardedly, since in an equation such as the fully covariant Eq. (jl(Jj) . one cannot remove 
the dependence on the Ricci scalar by any coordinate transformation whatsoever, and we thus define the 
equivalence principle not as the statement that all gravitational effects at any point P can be simulated 
by an accelerating coordinate system in flat spacetime (viz. that test particles unambiguously move on 
geodesies and obey Eq. (JHJ) and none other), but rather that no matter what propagation equation is 
to be used for test particles, the appropriate acceleration for them is D 2 x x /Dr 2 . 3 With Eq. (JH) thus 
showing the role of coordinate invariance in flat spacetime, with Eqs. (jSJ), ©, and (fTT)|) exhibiting the 
equivalence principle in curved spacetime, and with all of these equations being independent of the Ein- 
stein equations of Eq. (0), the logical independence of the general covariance principle, the equivalence 
principle and the Einstein equations is thus established. Hence, any metric theory of gravity in which 
the action is a general coordinate scalar and the metric is the gravitational field will thus automatically 
obey both the relativity principle and the equivalence principle, no matter whether or not the Einstein 
equations are to be imposed as well. 4 

2.3 Massless field motion 

Absent from the above discussion is the question of whether real, as opposed to test, particles actually 
obey curved space propagation equations such as Eq. (jSJ) at all, and whether such a discussion should 
apply to massless particles as well since for them both m and the proper time dr vanish identically, 
with equations such as Eq. (jSJ) becoming meaningless. Now it turns out that these two issues are 
not actually independent, as both reduce to the question of how waves rather than particles couple to 
gravity, since light is described by a wave equation, and elementary particles are actually taken to be 
the quanta associated with quantized fields which are also described by wave equations. And since the 
discussion will be of relevance for the exploration of the energy-momentum tensor to be given below, 
we present it in some detail now. For simplicity we look at the standard minimally coupled curved 
spacetime massless Klein-Gordon scalar field with wave equation 



where S' ,fl denotes the contravariant derivative dS/dx^. If for the scalar field we introduce a eikonal 
phase T{x) via S(x) = exp(iT(x)), the scalar T(x) is then found to obey the equation 



3 This particular formulation of the equivalence principle does no violence to observation, since Eotvos experiment type 
testing of the equivalence principle is made in Ricci-flat Schwarzschild geometries where all Ricci tensor or Ricci scalar 
dependent terms (such as for instance those exhibited in Eq. 1)10(1 ) are simply absent, with such tests (and in fact any 
which involve the Schwarzschild geometry) not being able to distinguish between Eqs. I@J and 1(10(1 . 

4 To sharpen this point, we note that it could have been the case that the resolution of the conflict between gravity 
and special relativity could have been through the introduction of the gravitational force not as a geometric entity at 
all, but rather as an analog of the way the Lorentz force is introduced in Eq. JHJ. In such a case, in an accelerating 
coordinate system one would still need to use the acceleration D 2 x x /Dt 2 and not the ordinary d 2 x x /dr 2 . However, if 
the gravitational field were to be treated the same way as the electromagnetic field, left open would then be the issue of 
whether physics is to be conducted in flat space or curved space, i.e. left open would be the question of what does then 
fix the Riemann tensor. There would then have to be some additional equation which would fix it, and curvature would 
still have to be recognized as having true, non-coordinate artifact effects on particles if the Riemann tensor were then 
found to be non-zero. Taking such curvature to be associated with the gravitational field (rather than with some further 
field) is of course the most economical, though doing so would not oblige gravitational effects to only be felt through 
D 2 x x I Dt 2 , and would not preclude some gravitational Lorentz force type term (such as the one exhibited in Eq. 1(10(0 
from appearing as well. 



(13) 



T ;/i T ;M - iT'^ = 



(14) 
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an equation which reduces to T ;A *T ;M = in the short wavelength limit. From the relation T ;M T ;At;I/ = 
which then ensues, it follows that in the short wavelength limit the phase T(x) obeys 

T^, u = T^[T.^ + d v T.^ - d„T. u ] = T'^T,^ + d v d„T - d^T] = T^T.^ = . (15) 

Since normals to wavefronts obey the eikonal relation 

T ;m = ^! = p ( 16 ) 
dq 

where q is a convenient scalar affine parameter which measures distance along the normals and k^ is 
the wave vector of the wave, on noting that (dx^ 1 / dq) {d / dx^ 1 ) = d/dq we thus obtain 

_ , , _ \ Lb Jb „ \ Lb Jb LbJb „ , > 

« A , = ^ + ^-- = 0, (17) 

a condition which we recognize as being the massless particle geodesic equation, with rays then precisely 
being found to be geodesic in the eikonal limit. Since the discussion given earlier of the coordinate 
dependence of the Christoffel symbols was purely geometric, we thus see that once rays such as light 
rays are geodesic, they immediately obey the the equivalence principle, 5 with phenomena such as the 
gravitational bending of light then immediately following. 

Now while we do obtain strict geodesic motion for rays when we eikonalize the minimally cou- 
pled Klein-Gordon equation, the situation becomes somewhat different if we consider a non-minimally 
coupled Klein-Gordon equation instead. Thus, on replacing Eq. (fT3j) by 

S'^ + ^SR a a = , R^ T ^0 , (18) 
the curved space eikonal equation then takes the form 

r ,^_i E a a = Q ; (19) 

with Eq. (fTTjl being replaced by the Ricci scalar dependent 

d x „ \ dx^ dx £ , % . \ / \ 

in the short wavelength limit. 

A similar dependence on the Ricci scalar or tensor is also obtained for massless spin one-half and 
spin one fields. Specifically, even though the massless Dirac equation in curved space, viz. ^^(x)^ + 
F^ijj^x) = (where T^(x) is the fermion spin connection) contains no explicit direct dependence 
on the Ricci tensor, nonetheless the second order differential equation which the fermion field then 
also obeys is found to take the form [<9 M + T M ][<9 M + T^]ip(x) + (l/4)R a ol i(:(x) = 0. Likewise, even 
though the curved space Maxwell equations, viz. F^ u . u = 0, F^ U] \ + F\^ u + F u \- ifl = also possess no 
direct coupling to the Ricci tensor, manipulation of the Maxwell equations leads to the second order 
g ap F fiV . a . p +Fi MX R v a -F' Ja R^ a = 0, and thus to the second order equation g a f 3 A l ,. a .p-A a . a . fi +A a R fMX = 
for the vector potential A 11 . Characteristic of all these massless field equations then is the emergence 
of an explicit dependence on the Ricci tensor, and thus of some non-geodesic motion analogous to 
that exhibited in Eq. (J2U|) when we eikonalize. 6 Now equations such as Eq. ()20|) will still obey the 

5 According to Eq. I|17() . an observer in Einstein's elevator would not be able to tell if a light ray is falling downwards 
under gravity or whether the elevator is accelerating upwards. 

6 We note that none of these particular curved space field equations involve the Riemann tensor, but only the Ricci 
tensor. This is fortunate since Schwarzschild geometry tests of gravity would be sensitive to the Riemann tensor, a tensor 
which in contrast to the Ricci tensor, does not vanish in a Schwarzschild geometry. 
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equivalence principle (as we have defined it above), since the four-acceleration which appears is still one 
which contains the non-tensor Christoffel symbols. Moreover, despite the fact that all of these particular 
curved space equations are non-geodesic, nonetheless, in the flat space limit all of them degenerate 
into the flat space geodesies, with the eikonal rays travelling on straight lines. Now in general, what is 
understood by covariantizing is to replace flat spacetime expressions by their curved space counterparts, 
with Eq. (j3J) for instance being replaced by Eq. (J3J). However, this is a very restrictive procedure, since 
Eqs. (jSJ) and (fTOj) both have the same flat space limit. The standard covariantizing prescription will 
thus fail to generate any terms which explicitly depend on curvature (terms which in some cases we see 
must be there), and thus constructing a curved space energy- momentum tensor purely by covariantizing 
a flat space one is not at all a general prescription, a issue we shall return to below when we discuss 
the curved space energy-momentum tensor in detail. 



2.4 Massive field motion 

A situation analogous to the above also obtains for massive fields. Specifically for the quantum- 
mechanical minimally coupled curved space massive Klein-Gordon equation, viz. 

S^___ 5 = , R pi/aT ^0 , (21) 

the substitution S(x) = exp(iP(x)/h) yields 

P ;M P ;M + mV = ihP'% . (22) 

In the eikonal (or the small K) approximation the ihP 11 .^ term can be dropped, so that the phase P(x) 
is then seen to obey the purely classical condition 

g tlu P' fX P' v + m 2 c 2 = . (23) 

We immediately recognize Eq. (J23j) as the covariant Hamilton- Jacobi equation of classical mechanics, 
an equation whose solution is none other than the stationary classical action J p^dx^ as evaluated 
between relevant end points. In the eikonal approximation then we can thus identify the wave phase 
P(x) as / p^dx^, with the phase derivative P ;M = d^P then being given as the particle momentum 
= medx^/dr, a four- vector momentum which accordingly has to obey 

V pV + mV = , (24) 

the familiar fully covariant particle energy-momentum relation. With covariant differentiation of Eq. 
immediately leading to the classical massive particle geodesic equation 



u \ d x \ dx^ 1 dx „ , 

^ = ^ + rJ— - = (25) 



(as obtained here with the proper time dr appropriate to massive particles and not the affine parameter 
q), we thus recover the well known result that the center of a quantum-mechanical wave packet follows 
the stationary classical trajectory. Further, since we may also reexpress the stationary / p^dx^ as 
—mc J dr, we see that we can also identify the quantum-mechanical eikonal phase as P(x) = —mc J dr, 
to thus enable us to make contact with the It action given in Eq. ©. 7 Though we have thus made 
contact with It, it is important to realize that we were only able to arrive at Eq. f!23|) after having 
started with the equation of motion of Eq. (|2*T|) . an equation whose own validity requires that stationary 

7 We make contact with the It of Eq. 10 since we start with the minimally coupled Eq. (|2TT) . 
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variation of the Klein-Gordon action from which it is derived had already been made, with only the 
stationary classical action actually being a solution to the Hamilton- Jacobi equation. The action It = 
—mc J dr as evaluated along the stationary classical path is thus a part of the solution to the wave 
equation, i.e. the output, rather than a part of the input. 8 Thus in the quantum-mechanical case we 
never need to assume as input the existence of any point particle action such as It at all. Rather, 
we need only assume the existence of equations such as the standard Klein-Gordon equation, with 
eikonalization then precisely putting particles onto classical geodesies just as desired. To conclude then, 
we see that not only does the equivalence principle hold for light (even though it has no inertial mass or 
gravitational mass at all) and hold for quantum-mechanical particles, we also see that the equivalence 
principle need not be intimately tied to the classical test particle action It at all. Given this analysis, 
we turn now to consideration of the curved space energy-momentum tensor. 



3 Structure of the energy-momentum tensor 
3.1 Kinematic perfect fluids 

At the time of the development of general relativity, the prevailing view of gravitational sources was 
that they were to be treated like billiard balls, i.e. as purely mechanical kinematic particles which 
carry energy and momentum; with the advent of general relativity then requiring that such energy and 
momentum be treated covariantly, so that the way the energy-momentum tensor was to be introduced 
in gravitational theory was to simply covariantize the appropriate flat spacetime expressions. Despite 
the subsequent realization that particles are very far from being such kinematic objects (elementary 
particles are now thought to get their masses dynamically via spontaneous symmetry breaking), and 
despite the fact that the now standard SU (3) x SU (2) x U (1) theory of strong, electromagnetic and weak 
interactions ascribes primacy to fields over particles, the kinematically prescribed energy-momentum 
tensor (rather than an SU(3) x SU(2) xU(l) based one) is nonetheless still used in treatments of gravity 
today. Apart from this already disturbing shortcoming (one we shall remedy below by constructing the 
energy-momentum tensor starting from fields rather than particles), an additional deficiency of a purely 
kinematic prescription is that since flat spacetime energy-momentum tensors had no need to know where 
the zero of energy was (in flat spacetime only energy differences are observable), their covariantizing left 
unidentified where the zero of energy might actually be; and since gravity couples to energy itself rather 
than to energy differences, this kinematic prescription is thus powerless to address the cosmological 
constant problem, an issue to which we shall return below. 

Historically there was of course good reason to treat particles kinematically, since such a treatment 
did lead to geodesic motion for the particles. Thus for the test particle action It of Eq. ©, its functional 
variation with respect to the metric allows one to define an energy-momentum tensor according to 9 

~ T =7 vF/9 / ^ \ x ~ y( T ))~] — J- > (26) 



(-9) 1/2 8g, u {-g) l ' 2 J " yy " dr dr 

with its covariant conservation {T^ u . u = 0) leading right back to the geodesic equation of Eq. (JSJ)- 10 
Moreover, an analogous situation is met when the energy-momentum tensor is taken to be a perfect 

8 Thus we cannot appeal to It to put particles on geodesies, since we already had to put them on the geodesies which 
follow from Eq. ill'MI) in order to get to It in the first place. While the use of an action such as It will suffice to obtain 
geodesic motion, as we thus see, its use is not at all necessary, with it being eikonalization of the quantum-mechanical 
wave equation which actually puts massive particles on geodesies. 

9 With our definition here and throughout of as T^ v — 2(SIt / Sg^) / (—g) 1 ^ 2 , it is cTqq which then has the dimension 
of an energy density rather than T 00 itself. 

10 While our ability to impose a conservation condition on T MI/ would follow from Eq. since the Einstein tensor 
Qfiiy _ jyiv _ (\l2)g^ v R a a obeys the Bianchi identity, the use of a conservation condition in no way requires the validity 
of Eq. QJ. Specifically, in any covariant theory of gravity in which the pure gravitational piece of the action, I gray, 
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fluid of the form 

= -[(p + pWU y + pg^} (27) 

with energy density p, pressure p and a fluid four- vector which is normalized to U^U^ = — 1. Specifically, 
for such a fluid, covariant conservation leads to 

[(p + p)EW + = [(P + P)U%U» + (p + p)U%U» + p»* = , (28) 

and thus to 

-[(p + p)U% + U^ = -(p + p)U%U^ = , (29) 
with the last equality following since U^U^. V = 0. With the insertion of Eq. (j2Sj) into (J25|) then yielding 

(p + p)^^ + p-Ag^ + u»u v ] = o , (30) 

viz. 

_ = _^ + ^_^_, (31) 

we see that the geodesic equation then emerges whenever the right-hand side of Eq. (|31j) is negligible, 
a situation which is for instance met when the fluid is composed of non-interacting pressureless dust. 

Since Eqs. (l2l)j) and (}3*T]) do lead to geodesic motion, it is generally thought that gravitational 
sources should thus be described in this way. However, even within an a priori kinematic perfect fluid 
framework, Eq. (|31j) would in no way need to be altered if to the perfect fluid of Eq. (J27|) we were 
to add an additional Tj^ XTRA which was itself independently covariantly conserved. Since a geometric 
tensor such as the Einstein tensor G^ u is both covariantly conserved and non-existent in flat spacetime, 
a curved space generalization of a flat space energy-momentum tensor which would include a term of 
the form T^ XTRA = G^ v would not affect Eq. (|31|) at all. Moreover, since the metric tensor g^ v is also 
covariantly conserved, the inclusion in T^ XTRA of a Kg^ v term with A constant would also leave Eq. 
(13 1|) untouched, and while such a term would not vanish in the flat spacetime limit, its presence in flat 
spacetime would only lead to a non-observable overall shift in the zero of energy. Restricting to the 
perfect fluid of Eq. (}2*Tj) is thus only sufficient to recover Eq. (}3~T]) . and not at all necessary. 

Within the kinematic perfect fluid framework, the use of such fluids as gravitational sources is greatly 
facilitated if some equation of state of the form p = wp can be prescribed where w would be a constant. 
To see what choices are suggested for w from flat space physics, we consider a relativistic flat space 
ideal iV particle classical gas of particles of mass m in a volume V at a temperature T. For this system 
the Helmholtz free energy A(V, T) is given as 

e -A(V,T)/NkT = yj d 3 pe -^+m^/kT ? (32) 

so that the pressure takes the simple form 

'~(£) T -T. 

while the internal energy U = pV evaluates in terms of Bessel functions as 

U = A-t(^) =3 NkT + N,n I ±MEl . (34) 
\ST) V K 2 (m/kT) 1 ' 



is a general coordinate scalar function of the metric, the quantity A^ v = (2/ '(—g) ) ' )(SIgrav I '^9^v) will, because of 
covariance, automatically obey A^ v '. v = 0, and through the gravitational equation of motion — T^ v then lead to 
the covariant conservation of T^ v . The use of the Einstein equations is only sufficient to yield T M!y .„ = 0, but not at all 
necessary, with the conservation ensuing in any general covariant pure metric theory of gravity. 
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In the high and low temperature limits (the radiation and matter eras) we then find that the expression 
for U simplifies to 



U 
V 

U Nm 3NkT Nm 



V V 2V V 



3NkT 


= 3P . 


m 






V 


uf 


3P 


Nm 


m 


V 





oo . (35) 



Consequently, while p and p are nicely proportional to each other in the high temperature radiation 
and the low temperature matter eras (where w(T — > oo) = 1/3 and w(T — > 0) = 0), we also see that 
in transition region between the two eras their relationship is altogether more complicated. Use of a 
p = wp equation of state would at best only be valid at temperatures which are very different from 
those of order m/K, though for massless particles it would be of course be valid to use p = p/3 at all 
temperatures, a point to which we shall return below. 



3.2 Perfect Robertson- Walker fluids 

In trying to develop equations of state in curved space, one should replace the partition function in 
Eq. (jH2J) by its curved space generalization (i.e. one should covariantize it, just as is proposed for 
T^ u itself), 11 and then follow the steps above to see what generalization of Eq. ()35|) might then ensue. 
However for curved backgrounds of high symmetry, the use of the isometry structure of the background 
can greatly simplify the discussion. Thus, for instance, for the Robertson- Walker geometry of relevance 
to cosmology, viz. 

f dr 2 \ 
ds 2 = -c 2 dt 2 + R 2 {t) L 1 _ kj . 2 ^ + r 2 dn 2 j , (36) 

the maximal 3-symmetry of the background entails that any rank two tensor such as the energy- 
momentum tensor itself must have the generic form 

T» v = [C(t) + D(t)]U^U u + D{t)gT , (37) 

to thus automatically be of a perfect fluid form with comoving fluid 4- vector = (1, 0, 0, 0) and general 
functions C and D which can only depend on the comoving time t. With the energy-momentum tensor 
being covariantly conserved, C and D have to obey 

| (fl . (0 + (38, 

with D and C thus being related according to 

D = -±{*C) . (39, 

While D and C must be also directly proportional to each other in the Robertson- Walker case since 
each only depends on the single parameter t, nonetheless, even though we therefore can set D/C = w 
in such cases, in general the quantity w could still be a function of t and need not necessarily be a 
constant. 12 In the above we have purposefully not identified C with a fluid p and D with a fluid p, since 
should the energy-momentum tensor actually consist of two types of fluid (this being the conventional 
dark matter plus dark energy picture), even if both of them have their own independent w according 



11 Typically one replaces (p 2 + m 2 ) 1 ' 2 / kT by (dx* 1 / dr)U u g ^ / kT \U V is a four-vector] and replaces J d 3 p by a sum over 
a complete set of basis modes associated with the propagation of a spinless massive particle in the chosen background. 
12 When w is a constant, we can then set C = 1/R 3 ^ W+1 \ 
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to pi = u>ipi, P2 = W2P2, the sum of their pressures would then obey pi + P2 = W\Pi + W2P2 and would 
not in general be proportional to pi + p 2 , with the total pi + p 2 and pi + p 2 obeying 

Pi+P2 = -^(R 3 (Pi+P2)) , (40) 

and with neither pi + p 2 or p± +p 2 then scaling as a power of i?. 13 However, if the full energy- momentum 
tensor is to describe radiation fluids alone, then no matter how many of them there might actually be 
in total, in such a situation the full must additionally obey the tracelessness condition = 0, to 
then unambiguously fix w = D/C to the unique value w = 1/3. 



3.3 Schwarzschild fluid sources 

While the isometry of the Robertson- Walker geometry does automatically lead us to the perfect fluid 
form given in Eq. (p?Tj) (though not necessarily to any particular equation of state unless = 0), the 
situation for lower symmetry backgrounds is not as straightforward. Specifically, for a standard static, 
spherically symmetric geometry of the form 

ds 2 = -B[r)c 2 dt 2 + A{r)dr 2 + r 2 dQ 2 , (41) 

the three Killing vector symmetry of the background entails that the most general energy-momentum 
tensor must have the generic diagonal form 

To = p{r)B{r) , T rr = p{r)A{r) , T ee = q{r)r 2 , T H = q{r)rhm 2 6 (42) 

with conservation condition 

dp (p + p)dB 2 

d-r + -2B-^ + r {p - q)= ° ' (43) 
with the function q(r) not at all being required to be equal to p(r). 14 Thus while a flat space static, 
spherically symmetric perfect fluid would have p = q and equation of state p — wp, it does not follow 
that a curved space one would as well, with this being a dynamical and not a kinematic issue whose 
resolution would require an evaluation of the covariant partition function in the background of Eq. 

AH}. 15 

Despite these considerations, it turns out that a quite a bit of the standard phenomenology associated 
with static, spherically symmetric sources still holds even if q and p are quite different from each other. 
Specifically, for the geometry of Eq. (}4*Tj) the components of the Ricci tensor obey 



Rqo _|_ Rrr Ree 1 



- (- 1 I 



2B 2A r 2 r 2 
while the components of the energy-momentum tensor obey 



r (I dB 1 dA\ 1 
R ™ = - 1 + 2a[b^-A^ + A ' (44) 



Too-^) = §( P + P + 2q) 



13 As well as needing to require that w\ and W2 both be constant, to secure the conventional p\ = 1 / R 3 ( Wl+1 \ P2 — 
1 ^ R 3( W2 +i) in the twQ fluid 

case additionally requires the separate covariant conservation of the energy-momentum tensor 
of each fluid, so that the two fluids do not then exchange energy and momentum with each other. A dynamics which is to 
secure this would have to be identified in any cosmological model (such as a dark matter plus quintessence fluid model) 
which uses two such fluids unless one of the two fluids just happens to have w = — 1, since this would correspond to a 
cosmological constant whose energy momentum tensor T^ v — —Kg 11 " is in fact independently conserved. 

14 Unlike the Robertson- Walker case, a static, spherically symmetric geometry is only spherically symmetric about a 
single point and not about all points in the spacetimc. 

15 Sincc q would be equal to p in the flat space limit, one might anticipate that for weak gravity the difference between 
q and p would still be small, though q could differ radically from p in the strong gravity black hole limit. 
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T rr - 9 -^p) = ^(p + p-2q) 
g ee T a a \ r 2 



277 { ' 



^0 n) + ^ (T rr - ^QrrT^ + i - ^T Q Q ) = p , (45) 



with the last expression conveniently being independent of q. If we now do impose the Einstein equations 
of Eq. (Q) , then in terms of the quantity 



M(r) = 47r / drr 2 p(r) (46) 



we immediately obtain 

to then find that B(r) obeys 



1C 

A-\r) = 1 - —Mir) , (47) 
r 



1 dB 2G 



B dr r (r - 2GM) 



(M + 4vrr 3 p) . (4f 



We recognize Eqs. (|47|) and (|48|) as being of precisely the same form as the standard expressions which 
are obtained (see e.g. 0) for A and (1 / B)dB / dr when q(r) is equal to p(r), though the substitution of 
these expressions into Eq. (}43|) leads to 

t + / P+ , P ^n (* + 4 " r3 ") +-&>-») = ». < 49 > 

dr r (r - 2GM) v ' r 

an equation which differs from the standard expression by the presence of the 2(p — q)/r term. If the 
matter density terminates at some finite r = R, then outside of the fluid the geometry is a standard 
exterior Schwarzschild geometry with metric 

IMC 

B(r> R) = A~Hr> R) = 1 -. (50) 

r 

Matching this exterior solution to the interior solution at r — R then yields for M the standard 

r R 

M = Arc drr 2 p(r) , (51) 
J o 

with the integration constant required for Eq. (}4~Sj) then being fixed to yield for B(r) the standard 
expression 

B(r)=e W (-2G f dr^±^] . (52) 
1 1 r(r~2GM)J 



As we thus see, the functional forms of Eqs. (|47|). (|51|) and (|52j) are completely unaffected by whether 
or not q is equal to p, with any difference between q and p only showing up in Eq. ([49)1 . As far as 
the geometry outside of the fluid is concerned, the structure of the exterior Schwarzschild metric is the 
standard one with the standard form for the total mass M as given in Eq. (|5ip. However, inside the 
fluid the dynamics could be different from the conventional treatment because of modifications to the 
equation of state. However, for weak gravity where p is already of order G, should the term of order 
G in dp/dr be given exactly by GpM/r 2 , the quantity p — q would then only begin in order G 2 and 
the standard lowest order in G hydrostatic treatment of sources such as stars would not be affected. 
Nonetheless, even if that is to be the case (something which is not immediately clear), for strong gravity 
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inside sources (where there are currently no data), any difference between p and q could have substantial 
consequences, and so we should not in general expect a static, spherically symmetric source to possess 
an energy-momentum tensor of the form T^ LU = [p(r) + p(r)]U tJ- U u + p{r)g^ v just because it does so in 
flat space. 16 



3.4 Scalar field fluid sources 

As we have thus seen, wisdom gained from experience with kinematic particle sources in flat space 
serves as a quite limited guide to the structure of gravitational sources in curved spacetime. However, 
that is not their only shortcoming, with their connection to the structure of the energy-momentum 
tensor which is suggested by fundamental theory being quite remote. Thus we need to discuss what 
is to be expected of the energy-momentum tensor in a theory in which the action is built out of fields 
rather than particles, and in which the fields develop masses by spontaneous symmetry breakdown. 
However, before going to the issue of dynamical masses, we first need to see how we can connect our 
above analysis of geodesic motion of eikonalized fields to the structure of the energy-momentum tensor. 

To illustrate what is involved, it is convenient to consider a massive complex flat spacetime scalar 
field, and to get its energy-momentum tensor (and to subsequently enforce a tracelessness condition 
for it when we restrict to massless fields), we take as action the non-minimally coupled curved space 
action 17 



-S'^S* + -m 2 SS* - ^SS*R> 1 U 
2 , M 2 12 



(53) 



Its variation with respect to the scalar field yields the equation of motion 



S^ + ^SR^-m 2 S = , (54) 
while its variation with respect the metric yields the energy-momentum tensor 

rp | 1 £ | ( Q Q* i ri Q* \ (3 2£) a £ / 

J-fxv — I 2 g I ^ >M ,f ' \ u ,fJb) g D ;a g ° ;n;v ° 

P / \ 1 P / 1 



+pr (S*S> a . a + SS*^. a ) - -g^m'SS* - ^SS* (i^ - -g^R a a ) . (55) 
With the trace of this energy-momentum tensor evaluating to 

T% = (£ - 1) (s^S^ + \S*S^ + ISS*'^ - m 2 SS* (56) 

in field configurations which obey Eq. ([54)1 . we see that the choice £ = 1, m = will enforce the 
tracelessness of T^ v . Bearing this in mind we thus set £ = l, 18 so that the flat space limit of the £ = 1 
theory then takes the form 

= ~ (dpSdrS* + d^Sd^S*) - ^ u d a Sd a S* - ~ (Sd^S* + Fd^S) 



16 In passing we additionally note that once q is not equal p, for a traceless fluid with T^„ = — p + p + 2q = 0, Eq. (|43(l 
reduces to dp/dr + (p + p) /2B) (dB / dr) + (3p — p)/r = 0, and dependent on how p depends on p, it may be possible for 
radiation to support a static, stable source. (When p = 2>p the only solution is p ~ 1/B 2 which would require B(r) to be 
singular at the point r = R where p vanishes.) 

17 To construct the correct energy-momentum tensor in flat space, it is necessary to first vary the curved space action 
with respect to the metric and then take the flat limit. 

18 When £ = 1, the coupling of the scalar field to the geometry is conformal, with the massless action Im = 
- J d^xl-g) 1 / 2 [S^S? /2 — 55*i? M M /12] being invariant under the local conformal transformation S(x) — > e' a ^S(x), 
g„ v (x) - e 2a ^ 9fiU (x). 
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+\g»» (S*d a d a S + Sd a d a S*) - l -g^m 2 SS* 



(d^Sd u S* + d v SdpS* 



\g^d a Sd a S* - \ (Sd,d v S* + S*d^d v S) - \g^m 2 SS* . (57) 





In a plane wave solution to the d^d^S = m 2 S wave equation of the form S(x) = e lk ' x jV 1 l 2 (E k ) x l 2 where 



k^k^ = —m 



E k = (k 2 + m 2 ) 1 / 2 and V is the 3- volume, T^ v then readily evaluates to 



T, 



fJIV 



k^k v 



(5f 



and even though the wavefront to the massive plane wave is geodesic, this particular energy-momentum 
tensor does not look anything like a perfect fluid form, since for a single massive plane wave k^ = 
(E k , 0, 0, k) propagating geodetically in the z-direction T^ v evaluates to 



T = — 



1 

V 



-k \ 









\-k k 2 /E k J 



(59) 



To get a perfect fluid we need to incoherently add an entire family of these plane waves. 19 Thus if we 
take a set of six plane waves moving in the ±x, ±y, ±z directions, all with the same k = \k\ and E k , 
viz. k» = (E k , k, 0, 0), k» = (E k , -k, 0, 0), — (E k , 0, k, 0), k** — (E k , 0, —k, 0), k» = (E k , 0, 0, k), and 
k^ = (Ek, 0, 0, —k), and then incoherently add up their individual contributions to T MZ ,, we obtain 



■oo 



V 



- yu 



2k 2 
EkV 



rpfj, 



Qm 
EkV 



(60) 



We recognize Eq. (JoUj) as being of precisely the perfect fluid form with p = QE k /V, p = 2k 2 /E k V. As 
such, we not only get a perfect fluid form, we even see that in the event that the mass is introduced 
kinematically as in the action of Eq. (|53|L then for small k m we obtain p p, with the effective 
p/p = w being zero, while for very large k (or equivalently for zero m) we obtain p = p/3. As 
such the above procedure shows how to obtain a perfect fluid form starting from field theory, with 
its generalization to curved space requiring an equivalent incoherent averaging over the basis modes 
associated with the curved space wave equations appropriate to the curved space backgrounds of interest, 
a procedure which we indicated earlier is not guaranteed to automatically yield the straightforward 
T xx = T yy = T zz condition found in flat space. 



3.5 Dynamical fluid sources 

While the above analysis shows how the perfect fluid form for massive particles can be obtained from 
field theory, and shows why the test particle action and energy-momentum tensor of Eqs. (JHJ) and 
fl20(l are not at all germane to the issue, this analysis still fails to take into consideration the fact 
that elementary particle masses are not kinematic, but rather that they are acquired dynamically by 
spontaneous breakdown. To investigate what is to happen in the dynamical mass case, it is convenient 
to consider a spin one-half matter field fermion ip(x) which is to get its mass through a real spin- zero 
Higgs scalar boson field S(x). In order to illustrate the difference between dynamic and kinematic 
masses in the sharpest way possible, we take the action to possess no intrinsic mass scales and the 
energy-momentum tensor to be traceless. We thus give neither the fermion nor the scalar field any 

19 This is equivalent to using (the zero temperature limit of) the partition function discussed earlier. 
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kinematic mass at all, and in order to secure the tracelessness of the energy-momentum tensor couple 
the scalar field conformally to gravity, to thus yield as curved space matter action 20 



- / d A x{-g) 



1/2 



" 12" " M 



(61) 



where h and A are dimensionless coupling constants. 21 Variation of this action with respect to ip(x) and 
S(x) yields the equations of motion 



tY(x)[d fl + T ll (x)}ilj-hS^ = 



(62) 



and 



1 



S^„ + ~SR» U - AXS 3 + h^ = 



(63) 



while variation with respect to the metric yields (without use of any equation of motion) an energy- 
momentum tensor of the form 



T, 



2 111 

i%l)<~1^x)[d u + T v {x)]i) + -S-f.S-y - -g^S' a S- a - -SS lK „ + -g^SS ? 



> 2 1 '•' 



— n R a 
/if ' ^''/"' 1 a 



\S A + iipj a (x)[d a + T a (x)}ip - hS-M 



(64) 



with use of the matter field equations of motion then permitting us to rewrite T^ v as 

2 11 

T^ v = ii)j^(x)[d v + T v (x)]ip + -S- ifl S. u - -g^ u S> a S. a - -SS ;(t .„ 

+-^9^SS' a . a - ^S 2 (r^ - -g^R a ^j - -g^hS^ijj . 



(65) 



Additional use of the matter field equations of motion then confirms that this energy-momentum tensor 
is indeed traceless. 

In the presence of a spontaneously broken non-zero constant expectation value Sq for the scalar 
field, the energy-momentum tensor is then found to simplify to 



T^v = iij)lvi(x)[d v + T u (x)]iJ) - jg^hSoipip - (r^ - ^g fMU R a a 



with flat space limit 



T, 



fJ,U 



^7/AV> - -7]^hS ipip 



(66) 



(67) 



20 In Eq. I jfiip the general- relativistic Dirac matrices 7 M (x) and the fermion spin connection r M (x) are defined as 7 M (x) = 
V < j 1 (x)j a and = [y l '(x),d fl ^(x)]/8—[ r Y' / (x),^ (T (x)]rZ iV /8wheveV^'(x) is a vierbein and 7 a is a special-relativistic Dirac 
gamma matrix which, with the ds 2 



{dx ) 2 + {dx 1 Y + {dx 2 Y + {dx i Y = r] ab dx a dx b metric, obeys 7 Q 7 6 + 7 6 7 a = -2 v ab 



while ip is given by ip = ip'D where I? is a Hermitian flat spacetime Dirac matrix which effects D^ a D 1 = 7 Q t. In our 
notation, here and throughout Hermiticity in the fermion kinetic energy sector is implicit, with i-0 7 ' i (x)[9 M + T^x)]^ 

denoting (i/2)ip^(x)[d IJ , + T^(x)]ip - (i/2)Jj[d IJ , +T ^(x)]^ {x)xp = (i/2,)^j t "(x)[d fl + F^x)]^ + Hermitian conjugate, 
and i-i/5 7M (x)[^ + Y v {x)]ij) denoting (i/4)-0 7M (a:)[9y + T u {x)]^ + (i/4)-0 7l/ (a;)[c} At + T^(x)]ijj + Hermitian conjugate in 
the associated energy-momentum tensor constructed below as the symmetric T^ v — 2g lia g v p{8lM / 8g a p) / {—g) 1 ^ 2 = 
-{l/2)[{V m {x)/\V{x)\){8I M /5V^{x)) + (V ua (x)/\V(x)\)(5I M /5V£(x))} where \V(x)\ = {-g) 1 ' 2 . In this construction 
we note that Sg^ = V^ a SV v a + SV^Vf, S\V(x)\ = -\V(x)\V£(x)8V£(x), and 5V ac = -V^VfSVf, to thereby 
yield (y va {x)/\V{x)\){5/5V^{x))[j d^x{-gY l2 ^{x)[d a +Y a {xM ^ ^(x)[d^ +T^(xM + (1/8)^ (x)^,^}).^ - 
9v^l a (x)[d a + r a (x)]ip, from which Eq. (|64() then follows. (The author wishes to thank Dr. R. K. Nesbet for helpful 
comments on the material presented here.) 

21 As such, the action of Eq. (|61[) is the most general curved space matter action for the ip{x) and S(x) fields which 
is invariant under the local conformal transformation S(x) — > e~"( x ) S(x), ^p(x) — > e~ 3a ^ x >'^ip{x), ^{x) — > e~ 3Q ^)/ 2 -0(x), 
g^{x) - e 2a ^g^ v (x). 
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With the fermion now obeying 



i^d^ - hS^ = 



(6£ 



the tracelessness of the flat space energy-momentum tensor of Eq. (J67|) is manifest. The quantization 
of the flat space limit theory is straightforward since Eq. (|68|) just describes a free fermion with mass 
m = hSo, and yields one particle plane wave eigenstates \k) of four momentum = (Ek,k) where 
Ejt = (k 2 + m 2 ) 1//2 . For a positive energy, spin-up, Dirac spinor propagating in the z-direction, in analog 
to Eq. (|59|) we obtain the matrix elements 



T = — 



with trace 



(E k -k \ 





V— A; k 2 /E k J 



— 



1 

V 



( 



-m 2 /AE k 




k' 2 





m 2 /AE k 








m 2 /AE k 










\ 



(69) 



m 2 /AE k j 



m 



-E k + — + — 







(70) 



E k Ek 

In Eq. (I69|) we recognize a two-component structure to the fermion energy-momentum tensor, a standard 
kinematic piece in which (k\T 0Q \k) = E k /V, (k\T 33 \k) = k 2 /EkV, and a dynamic part coming from the 
symmetry breaking in which (k\Too\k) = —m 2 /4E k , (k\Tn\k) = (k\T22\k) = (A: | X33 1 A;) = m 2 /AE k . On 
incoherently averaging over the directions of k, the energy momentum tensor is then found to take the 
form 



Tp> = (p + p)U^U v + pr/^ + Ar] 



where 



P 



2k 2 



p 



A 



3m 



V ' r VE k 2VE k 
and where the tracelessness of T^ u is enforced by the relation 

3p - p + 4A = . 



(71) 



(72) 



(73) 



As an energy-momentum tensor, Eq. ()71|) provides an explicit example of a phenomenon we referred 
to earlier, namely that it is possible to add on to a kinematic energy-momentum tensor an additional 
tensor which is itself covariantly conserved without affecting the covariant conservation of the kinematic 
energy-momentum tensor itself. As well as showing how a cosmological constant type term A can 
naturally arise in dynamical mass theories, the great virtue of using Eq. ()71|) is that the tracelessness 
condition of Eq. (|T3*j) constrains the value of A to be neither smaller nor larger than p — 3p, an issue 
which we shall revisit below when we discuss the dark energy problem, with the tracelessness condition 
thus being seen to give us control of the cosmological constant. 

It is instructive to compare this dynamical theory of fermion masses with a strictly kinematic fermion 
mass theory in which the action is given by 



I M = - d\{-g) 1 ' 2 ^(x)^ + r M (x)]V - mW> 



with the flat space fermion wave equation being given by 

i^d^ip — mip = 
the flat space energy-momentum tensor being given by 



(74) 



(75) 



(76) 
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and the trace being given by the non-zero 



(77) 



This time, an incoherent averaging over plane wave states yields 



6E k 



2k 2 



(78) 



T, 



p = 



V 



VE k 



with trace = 3p — p ^ 0. With both Eq. (fTTj) and Eq. ((75)) leading to the covariant conservation 
of the kinematic (p + p)U fl U u + pr/^, we see that both energy- momentum tensors lead to eikonalized 
geodesic motion, with the validity of the equivalence principle in no way requiring the use of Eq. (J78|) . 22 
The key distinction between Eqs. (fTTj) and (|75j) is that Eq. (fTTj) contains not just the energy on the 
fermion field, but also that of the Higgs field that gave it its mass, an energy which couples to gravity. 
In dynamical theories of mass generation, gravity is thus sensitive to the Higgs field associated with the 
fermion, while in kinematic case it of course is not. 

3.6 Implications of elementary particle physics 

While it is conventional to use Eq. ()78)1 and its analogs in current standard model cosmological studies, 
it is actually Eq. (fTTj) and its analogs which is suggested by elementary particle physics. Specifically, 
in the standard SU(3) x SU(2) x U(l) theory of strong, electromagnetic and weak interactions, the 
only mass scale present in the fundamental Lagrangian is the wrong sign, tachyonic mass term in 
the Higgs potential V(S) = XS 4 — p, 2 S 2 /2 with all fermions and gauge bosons getting their masses 
dynamically through a non-vanishing Higgs field expectation value, and all couplings between fields 
being dimensionless. As regards the Higgs field, it could be a fundamental field with a bona fide 
fundamental tachyonic mass, or it could be a dynamical manifestation of an underlying symmetry 
breaking through bilinear fermion condensates, with S(x) then only being a Ginzburg-Landau type 
long range order parameter (and V(S) its effective Ginzburg-Landau action) which arises when the 
fermion bilinear takes a non-zero expectation value in a spontaneously broken vacuum. In such a case 
the underlying theory of fermions and gauge bosons would possess no intrinsic mass scales at all and 
would then have a traceless energy-momentum tensor. Beyond these two possibilities for the Higgs field, 
inspection of the action of Eq. (foT|) reveals yet another origin for the wrong sign mass term, namely that 
it could even arise from curvature in a theory with no fundamental mass scale at all, with W 1 serving 
as p? in a non-flat background. Thus either the energy-momentum tensor of SU(3) x SU(2) x £7(1) is 
traceless (dynamical Higgs) or its trace is given (fundamental Higgs) by the non-zero = [i 2 Sl (c.f. 
Eq. (jHHI)) when 5* is constant. At best then the trace can only be related to the Higgs field itself, and 
would thus not be given by the T M = 3p — p expected of the kinematic fluid of Eq. (J7HJ). The use of 
the perfect fluids commonly employed in standard gravity and cosmology would thus appear to be at 
variance with the standard SU(3) x SU(2) x U(l) model of particle physics. 23 Hence, even before we 

22 Analogously, if we take the scalar S(x) to be a constant in Eq. Eq. H12|) then reduces to Eq. (JSJ, with test 

particle motion being geodesic whether mass is kinematic or dynamic. 

23 In the standard particle physics model radiative corrections will lead to the generation of a trace anomaly in theories 
whose starting Lagrangian has no fundamental mass scales at all. However, such anomalies while then making the trace 
of the renormalized energy-momentum tensor non-zero, do not themselves introduce any new mass scale (the dimension 
four trace is proportional to terms which are quadratic in the curvature tensor) and would not in any way make the trace 
equal to the kinematic 3p — p required of Eq. I|78|l . Moreover, it is possible to actually cancel the trace anomaly altogether, 
by a judicious choice of fields, a judicious choice of geometric background, or by a renormalization group fixed point at 
which the coefficient of the trace anomaly is then zero. That the trace anomaly is to be cancelled would appear to be 
part of the standard cosmological wisdom anyway, since no such anomaly term contribution is ever included alongside 
the kinematic ordinary matter, kinematic dark matter and dark energy components which make up the energy content 
of the standard cosmological model, and so we shall not consider the trace anomaly any further here. 
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enter into the issue of the fact that the use of standard gravity in astrophysics and cosmology leads to 
the dark matter and dark energy problems, we see that the very use of such dark matter and dark energy 
sources at all in the kinematic way in which they are commonly used is already not in accord with the 
standard model of particle physics. With this word of caution in mind we shall now discuss the status 
of standard gravity with its conventional kinematic perfect fluid sources, and in looking for departures 
from standard gravity then allow for theories in which the exact energy-momentum tensor is traceless. 
Since the use of a traceless energy-momentum tensor in the Einstein equations of Eq. {1} would lead to 
a Ricci scalar which would have to vanish in every conceivable situation (to rapidly then bring standard 
gravity into conflict with data), when we do consider the source of gravity on the right-hand side of the 
gravitational equations of motion to be traceless, we will have to modify the left-hand side as well in 
a way which would make it automatically traceless as well. As we shall see, such a departure from the 
standard theory will readily allow us to resolve both the dark matter and the dark energy problems. 
Having now explored constraints on the right-hand side of Eq. (JTJ in some detail, we turn next to an 
exploration of its left-hand side, and begin first with the Newtonian gravity which was its antecedent. 



4 Newtonian gravity 
4.1 The Newtonian potential 

The Newtonian prescription for determining the non-relativistic potential <f>(r) at any point r due to a 
set of static, mass sources, rrij, at points fl is to sum over them according to 

N m-G 

i r ' t\ 

where G is Newton's constant, with the motions of material test particles then being determined via 

cP'r -> 

3P — V* . (80) 

As such, Eq. (|79|) contains the full content of Newton's law of gravity, and for any given set of gravi- 
tational sources, any candidate theory of gravity needs to recover the associated <p(r) in any kinematic 
region in which Eq. (J79|) has been confirmed, to the level of precision required by available data. To 
determine what are the appropriate kinematic regions requires first testing Newtonian predictions in a 
candidate region using sources which are already known and prescribed in advance of data, and after 
success has been achieved, and only after, is one then free to use the law again in that same region to 
infer the existence of other, previously unknown, sources. At the present time one can say that in this 
way Newton's law of gravity has been well-established on distance scales from the order of millimeters 
(the smallest distance scale on which there has been testing) out to distances of the order of 10 15 cm 
or so (viz. solar system distance scales). At much larger distances all tests which use only the sources 
which were known in advance of observation have been found to fail, with additional (dark matter) 
sources always having to be invoked after the fact, with the validity of Newton's law yet to be con- 
firmed on those distances. 24 The establishing of a given law in one kinematic region cannot be regarded 
as evidence for its validity in others. 

As regards solar system tests of Newton's law of gravity, these tests are spectacular, with the 
planets being found to move according to Eq. (fKUj) to very high accuracy when the sun and the planets 

24 Testing the validity of Newton's law on a distance scale such as galactic would require knowing the dark matter 
distribution in a galaxy in advance of measurement of the orbital velocities of the material in it. 
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themselves are taken to be the sources. To a very good approximation the planets move around the 
sun with a centripetal acceleration of the form 



V 2 _ GM r . 
r r 





2 



and have Keplerian fall-off of velocity with distance of the form 



G 1/2 M^! 2 

,1/2 



v = - ;; ° (82) 



where M & is the mass of the sun. In fact so reliable was this law found to be in the solar system when 
only known visible sources were used, that when the orbit of Uranus was found to depart slightly from 
the Newtonian expectation, a perturbation due to a then undetected nearby planet was proposed, with 
the subsequent detection of Neptune giving dramatic confirmation of Newton's law of gravity at that 
distance scale. 

While Newton's law had always been thought to hold on all solar and sub-solar distance scales, 
recently it has come in to question on very small laboratory scales of order millimeters or so. While the 
meter or so distance region had actually been quite throughly searched (though to no avail) when it was 
thought that there might be a so-called fifth force operative on those distance scales (fifth since it would 
be in addition to the strong, electromagnetic, weak and gravitational forces), recent studies of large 
extra dimension physics had led to a reopening of the issue with a focus on the even smaller millimeter 
region distance scale. While there has been a longstanding theoretical interest in the possible existence 
of additional spacetime dimensions beyond the four established ones, it had generally been presupposed 
that such dimensions would be truly microscopic, possibly being as small as the 10 -33 cm Planck length. 
In an attempt to resolve the longstanding hierarchy problem of understanding why there was such a 
huge disparity between the Mew = 10 3 GeV electroweak and the Mpi = 10 19 GeV gravitational mass 
scales, Arkani-Hamed, Dimopoulos and Dvali |Zj found a candidate extra dimension based solution in 
which rather than be microscopic, the additional dimensions beyond four would instead need to be of the 
order of millimeters. In such a situation gravitational flux lines could then spread out to such distances 
in the extra dimensions, and therefore lead to modifications to standard gravity at the millimeter level. 
Following on the work of Arkani-Hamed, Dimopoulos and Dvali, Randall and Sundrum [HJ E] found an 
alternate way to address the hierarchy problem which required the geometry in the extra dimensions 
not to be the Minkowski one that had always been assumed in higher dimensional theories, but rather 
to be anti-de Sitter. In such a situation the extra dimensions could then not only be large but even 
infinite in size, with the curvature of the anti-de Sitter space acting as a sort of refractive medium 
which would sharply inhibit the penetration of gravitational flux lines into the higher dimensions. In 
such theories there could again be departures from Newton's law at the millimeter level, which thus 
prompted a renewed search for possible millimeter region departures. Since for the moment, none has 
yet actually been detected [TU], the issue should be regarded as open, and so we shall assume here 
the validity of Eq. (J75J) on all solar and sub-solar distance scales including millimeter ones and below, 
with the deriving of the Newtonian phenomenology on such scales thus being set as a requisite for any 
candidate gravitational theory. 



4.2 The second order and fourth order Poisson equations 

While Eq. (J75j) contains the full content of Newton's law of gravity, in order to actually perform the 
needed sum over sources, it very convenient to recast Newton's law in the form of a second order Poisson 
equation. Thus if we set 

V 2 0(r") = g(r) , (83) 
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we can write the potential as 



^ = -H d3p iSi • (84) 



with the potential exterior and interior to a spherically symmetric static source of radius R then being 
given by 

M r > R) = -- [ R dr'r' 2 g(r') , Mr < R) = -- f dr'r' 2 g(r') - [* dr'r'gir') . (85) 
r Jo r Jo J r 

However, Eq. ()8Hj) is not the only Poisson type equation which will lead to a 1/r potential. Thus if we 
for instance consider the fourth order 

VV(f) = h{r) , (86) 

we can write the potential as 

0(f) = / d 3 r'h(r')\f-r'\ , (87) 
on J 

with the potential exterior and interior to a spherically symmetric static source of radius R then being 
given by 

Mr > R) = -- [ R dr'r' 2 h(r') - — [* ' dr'r'^hir') , 
2 jo 6r Jo 



r 



(r<R) = -- f dr'r' 2 h{r') - — f dr'r' A h{r') - - [* dr'r' 3 h(r') - — [* dr'r'h{r') 
2 Jo 6r jo 2 Jr 6 Jr 



Moreover, recovering a 1/r potential is not even restricted to the above choices of Poisson equation, 
since the sixth order V 6 0(f) = k(r) would yield an exterior potential of the generic form <p ~ 1/r + r+r 3 , 
with the pattern repeating for all higher even number of derivative Poisson equations. Characteristic of 
all of these possible Poisson equations is the fact that we recover a 1/r term for each and every one of 
them, and in all of them the terms which depart from a pure 1/r form are all only important at larger 
distances (viz. just the distances on which the dark matter problem is encountered), with there being 
no modifications to the 1/r law at very small distances. Providing that the r, r 3 and so on type terms 
are all negligible on solar system distance scales, all of these Poisson equations will thus reproduce the 
standard Newtonian solar system phenomenology. 

Now at first glance the expressions for the coefficients of the 1/r terms in Eqs. ()85|) and (jHHj) appear 
to differ since they are given as different moments of the source. However, this is not of concern since 
measurements of <f>(r) in the r > R region where there are data can never uncover any information 
regarding the behavior of the integrand in Eq. (J85)) in the r < R region. What is required of the 1/r 
term in the r > R region is only that it conform with the large r behavior required by the original 
Newton law of Eq. (|79J) (and thus the all r > R behavior since 1/r is the unique behavior of the Newton 
term in the entire exterior region), viz. that for a set of sources all with the same mass m the potential 
behave as 

NmG , , 

Mr > R) = , (89) 

r 

i.e. that the coefficient of the 1/r term be linear in the number of sources N. To achieve such a linearity 
in the second order Poisson equation case one can choose to set 

N S(r - r ) 

g(r < R)=mGY y „ l > (90) 
£i r 

in the r < R region, though as we had just noted, this is not mandated by Eq. (|79j). (While use of this 
particular g(r < R) is sufficient to give Eq. (jSHJ, it is not necessary.) What is necessary is only that 
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there be iV discrete sources, with Eq. (179)1 only counting the number of them in the r > R region. To 
achieve exactly the same result through the use of the fourth order Eq. ()88|) . it is convenient, though 
not mandatory, to consider as source 25 



h(r < R) 



with its insertion in Eq. 



i=i ' 
yielding 



3(3c 



<f)(r > R) 



2 N 

E 



Nf3c 2 N<yc 2 r 



(91) 



(92) 



As we see, the coefficient of the 1/r term is again linear in the number of sources, with Newton's law of 
gravity only counting the number of discrete sources that are present in the r < R region. Comparing 
Eqs. (jSUjl and (j9"2~|) . we see that we can once and for all define (3 = mG/c 2 at the individual discrete 
source level (a level which could be microscopic), without ever needing to decompose (3c 2 into separate 
m and G pieces, with only the Schwarzschild radius of the source (viz. IraGjc 2 in the second order 
case and 2(3 in the fourth order) ever being measurable gravitationally. 26 

It is important to stress that it is the presence of two independent singularities in the source of Eq. 
(l9*Tj) which leads to the logical independence of the (3 and 7 coefficients, since if we could approximate 
the source h{r) by a constant source h{r) = h, we would instead have to conclude that the coefficients 
of the two potential terms would be related according to 



(r > R) 



NhR 5 NhRh 



30r 



6 



(93) 



and thus have a ratio of order the radius squared of the source. This is to be contrasted with the choice 
of a constant source g(r) = g in the second order case where the potential is then given by 



(r > R) 



NgR 3 
3r 



(94) 



Since a comparison of Eqs. ()89|) and (pUJ) would entail that one can make the identification mG = gR 3 /3, 
a view of gravitational sources has developed that macroscopic gravitational sources can be treated as 
being continuous rather as being a collection of independent discrete microscopic sources. However, 
apart from not at all being mandated by the successful use of Eq. (|5U|l . the issue is in fact not even 
addressable if one only makes measurements in the exterior r > R region. However, in the fourth order 
Poisson case the issue does become relevant, because here one does measure more than just one moment 
of the source, with the higher derivative theory thus being able to probe deeper into the source. Thus 
within the fourth order theory, an experimental determination of the ratio of the coefficients of the 1/r 
and r terms which would indicate that the ratio of the two coefficients is very far from being of order the 
square of the radius of a given macroscopic source, would indicate that the macroscopic source would 



25 We choose the particular [V 2 — (r 2 /12)V 4 ](<5(r — ri)/r 2 ) source rather than the more straightforward V 2 (<5(r — ri)/r 2 ) 
one, since unlike the latter, the former one happens to be positive definite. An additional virtue of this particular choice of 
source is that it only couples to the fourth moment integral and not to the second moment one; and with the S(r — rj)/r 
source only coupling to the second moment integral and not to the fourth moment one, the logical independence of the 
two moments is thus established. 

26 Even if the microscopic f3c 2 were not to be equal to the product of the conventionally defined m and G (not that we 
know what the gravitational coupling of a microscopic source is), the use of the coefficient N(3c 2 would then lead to a 
estimate for the number of atoms in a macroscopic source such as the sun which would be different from the conventional 
one, a point which is not of concern since without appealing to the mG coefficient, we do not know how many atoms there 
anyway are in the sun, with the standard estimate of the number being based on the a priori use of the mG coefficient. 
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have to be composed of microscopic components which are discrete. Quite remarkably then, higher 
order derivative theories open up the possibility of establishing discreteness at the microscopic level 
from macroscopic measurements alone. 27 Moreover, below we shall actually apply the potential of Eq. 
(|92|) to galactic rotation curves, to find that the f3* and 7* coefficients of the fourth order theory stellar 
potential 

0*(r>£O = -— + ^ (95) 
r 2 

actually are very far from possessing a ratio (3* /j* which is anything like the square of a typical stellar 
radius (numerically we find (P*/^*) 1 ^ 2 ~ 10 23 cm, with the j*c 2 r/2 term first becoming competitive 
with the /3*c 2 /r term on galactic distance scales). Within the framework of a gravity which is based on 
fourth order equations then we can conclude that matter must be discrete at the microscopic level. 

While the second order Poisson equation of Eq. (J83)l does yield the Newtonian potential of Eq. (|55j) 
as solution, there is actually a critical difference between the use of the second order Newton potential 
and Poisson equation on the one hand and the use of their fourth order counterparts given as the 
potential of Eq. ()92)) and the Poisson equation of Eq. (jHfij) on the other. Specifically, while the fourth 
order potential of Eq. ((QSj) does nicely reduce to the second order potential of Eq. (fKH|) in the small 
r limit, the fourth order Poisson equation itself never reduces to the second order one in any limit at 
all. It is thus possible for solutions to two different equations to approximate each other very closely 
in some kinematic region even while the equations to which they are solutions never approximate each 
other at all. The equations which describe any candidate alternate theory of gravity thus do not at 
all need to reduce to those which describe the standard theory. Rather, it is only the solutions to the 
alternate theory which have to recover the solutions to the standard theory in any of the kinematic 
regions in which the standard theory solutions have themselves been tested. 28 Outside of such regions 
the alternate theory can then provide for departures from the standard theory expectations, something 
we will capitalize on below to explain galactic rotation curve systematics without the need to invoke 
dark matter. The essential difference between the second order Poisson equation and the fourth order 
one, is that once the second order Poisson equation is specified as being the fundamental equation of 
motion, there is no possibility to ever get any departure from Newton's law on any distance scale at 
all, something which is allowed for in a higher order theory. Since the Einstein equations of Eq. (|T|) 
were explicitly constructed so that their non-relativistic weak gravity limit would be the second order 
Poisson equation, the Einstein equations thus lock non-relativistic weak gravity into being Newtonian 
on all distance scales, and to thus be so even on distances well beyond the solar system ones where 
Newton's law of gravity was originally established. In this sense then the dark matter problem could 
originate in the lack of reliability of the extrapolation of standard solar system wisdom to altogether 
larger distance scales. 



4.3 Universal acceleration scale and the MOND theory 

While the above discussion has focused on the standard Newtonian picture and those allowed departures 
from it which are associated with a critical distance scale, an alternate form of allowed departure has 
been suggested by Milgrom ^U^lEI], one which has proven to be remarkably instructive and fruitful. 
Specifically, on noting that the centripetal accelerations of particles in orbits around galaxies were 

27 As a historical aside, we note that since the use of higher order gravity can thus provide for a macroscopic manifestation 
of the existence of atoms, it thus meets an objective which was first sought by Boltzmann, one which was part of his 
motivation for developing kinetic theory (and which was subsequently achieved by a determination of a finite value for 
Avogadro's number by macroscopic means). 

28 If the alternate equations of motion did reduce to the standard ones in some limit, then of course so would their 
solutions. However, the converse is not true, since the alternate theory solutions are still able to reduce to the standard 
theory solutions even if there is no such parallel for the equations of motion themselves. 
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typically much smaller than those associated with the orbits of planets around the sun or satellites 
around the earth, Milgrom suggested that the determining factor was not in fact a distance scale at 
all but rather an acceleration one, with departures from the standard wisdom occurring whenever 
accelerations of particles dropped below some critical universal acceleration scale a®. Quite remarkably, 
Milgrom found that it is actually possible to choose such an ao so that its effects would be of significance 
for galaxies and yet leave solar system phenomenology untouched. 29 In contrast to proposals which seek 
to modify the right-hand side (viz. the gravitational side) of Eq. 1)80)1. Milgrom instead proposed to 
modify the left-hand side (viz. the inertial side) with the proposal thus being known as MOND (modified 
Newtonian dynamics). Specifically Milgrom replaced Eq. (|KT)|) by 

n(-)s=f (96) 
\a J 

where a = d 2 f/dt 2 is the ordinary acceleration, / = — V0 is the gravitational force and n(a/ao) is the 
modification. In order for this modification to reduce to the standard Newton law when a 3> a , the 
function fj,(x) has to behave as fi(x ^> 1) — > 1. Guided by the fact that many prominent spiral galaxies 
possess rotational velocity curves which are flat (i.e. the rotational velocity v is independent of the 
distance R from the center of the galaxy), and that the velocities of these particular galaxies obey the 
empirical Tully-Fisher law which relates v to the total luminosity L of the galaxy as v 4 ~ L, Milgrom 
further suggested that at very small x the function /x(x) should behave as fi(x <C 1) — > x, since the 
centripetal accelerations a = v 2 /R <C ao of orbits in a galaxy of mass M would then obey 



a 2 v 4 MG 



ao a R 2 R? 



(97) 



and thus yield the R independent v 4 = a MG. Then, if the mass involved was just the luminous mass 
of the galaxy, it would be related to the luminosity of the galaxy by M = (M/L)L where M/L is the 
mass to light ratio of the galaxy, with both the flatness of galactic rotation curves and the Tully-Fisher 
relation thereby being obtained from luminous matter alone. 

As a proposal, the relation of Eq. (jHEJ) is perhaps a little bit too strong a departure from the standard 
picture, since as such it would modify the inertial side of Newton's second law of motion not just for 
gravitational processes but for non-gravitational ones as well, processes for which there are no known 
problems. Additionally, since the proposal is a modification of the inertial side of Newton's second law 
of motion, in an accelerating coordinate system Eq. 1)96)1 would not readily recover Eqs. (jlj and (0). To 
get a sense of what might instead happen, one possible way to effect a covariant generalization of Eq. 
dHIIJ) is to replace the acceleration a by the contravariant D 2 x x /Dr 2 of Eqs. (jlj) and (0), to introduce a 
four-vector V x which reduces to V x = (1, 0, 0, 0) in a non-accelerating coordinate system, and to take 
n{x) to be a general coordinate scalar function of 

VxD2%X 

x = — 7T 98 

a Dr 2 v ; 

where the parameter ao is then a general coordinate scalar with the dimensions of acceleration. While 
the presence of such a V\ would lead to frame- dependent effects, such departures from standard wisdom 
need not necessarily be in conflict with observation if fi(x) is sensibly close to one in the kinematic 
regime where tests of the equivalence principle and preferred frame effects have so far been made. 

Nonetheless, since it would be more economical not to have to make so radical a conceptual departure 
from standard wisdom, an alternative procedure for introducing a universal acceleration ao would be 



29 Phcnomcnologically ao was found to be of order 10 8 cm/sec 2 , to be compared with an acceleration at the outer 
limits of the solar system of order 10~ 3 cm/sec 2 due to the pull of the sun. 
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to have it instead appear on the gravitational side of Eq. ([80)1 . with Eq. (JHUj) then being replaced by 




(99) 

instead. To get a sense of what the function v(y) might look like, we consider a particularly simple 
choice for fi(x) which meets its required large and small x limits, viz. 

M^)= ( 1+ ^2)l/2 ' ( 10Q ) 

to then find that v(y) and a are given by 

U + 2/ ° } j / ■ (102) 

With the two relevant limits of the function u(y) being u(y ^> 1) — > 1, z^(y < 1) -» l/y 1 ^ 2 , we thus 
recover Eq. (jSUj) when / ^> a , and come right back to Eq. (J9~T|) when / < a , with Eqs. (J5fij) and 
f!99|) leading to the same gravitational phenomenology. Having now presented an analysis of Newtonian 
gravity and some possible allowed departures from it, we turn next to a discussion of Einstein gravity. 





5 Einstein gravity 

5.1 Einstein gravity for static sources and the Newtonian limit 

In his construction of a covariant theory of gravity, even though Einstein used a fundamental principle, 
viz. the equivalence principle, to identify the metric as the gravitational field, his choice for the 
dynamical equation of motion which was to then fix the metric, viz. Eq. (|T|). appears to have been based 
on a phenomenological consideration, namely that the theory recover the standard second order Poisson 
equation in the non-relativistic, weak gravity limit. As such, the second order Einstein equations can be 
viewed as being a covariant generalization of the second order Poisson equation of Eq. (|83p. replacing 
it by ten equations (reducible to six by gauge invariance) of which only the (00) component equation 
would be of order one in the standard non-relativistic, weak gravity limit in which the components of 
Tju, other than T 00 are taken to be negligible compared to T 00 itself. The great virtue of the Einstein 
equations is that in prescribing a covariant generalization of the second order Poisson equation, it not 
only incorporated Newtonian gravity, but also provided it with general relativistic corrections which 
could then be tested, leading to the three classic tests of general relativity, the gravitational red shift, 
gravitational bending of light, and the precession of planetary orbits. With the Einstein equations not 
only satisfying these three tests, but doing so in truly spectacular fashion, by and large the consensus 
in the community is that the issue is therefore considered settled, with the Einstein equations being the 
true equations of gravity which no longer need to be open to question. 

However, it is important to note that while Einstein's original construction was motivated by the 
second order Poisson equation, the application of Eq. (JTJ) to the three classic solar system tests is not 
actually sensitive to this particular aspect of the Einstein equations at all. Specifically, what is needed 
for the three tests is only knowledge of the metric exterior to a static, spherically symmetric source, viz. 
the Schwarzschild metric as given earlier as Eqs. (j4*Tj) and (|5U|). a metric which obeys R^ v = in the 
source-free region. As such, the exterior Schwarzschild solution is parameterized by a parameter MG 
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associated with a given source, and while this parameter can be determined from a matching of the 
exterior solution of Eq. (JoTty to the interior solution of Eq. (|4"7j) to yield Eq. (JoTj) . no use of Eq. (JoTj) is 
actually made in the analysis of the classic tests. While the = condition does indeed follow from 
Eq. (JTJ) in the source-free T^ v = region, Eq. (JTJ) does not follow from this fact, i.e. knowing that R^ 
vanishes in the source-free region is not sufficient to determine what is to be equal to in regions 
where it does not vanish. To emphasize the point, we note that while functional variation of the second 
order Einstein-Hilbert action 



c 3 



16%G 

with respect to the metric yields the Einstein tensor according to 



J d A x(-g) 1/2 R a a (103) 



16ttG 5I EH G » v = W v_\ rR * a j (104) 



variation of the fourth order general coordinate scalar actions 

i Wl = J dH-g) 1/2 [R a a } 2 (105) 



and 



respectively yield 



! w 2 



d 4 x(-g) 1/2 R a(3 R a P (106) 



(■^1/2^ = W W = 2 9" u (R a a f-,, - 2(iTJ^ - 2R\R^ + \g^{R\) 2 , (107) 

and 

(-^1/4^! = W & = \ g "^ Ra ^^ + ***** - ***** - - 2R " Prv ? + l^Ra P R aP ■ (108) 

Then, since the Schwarzschild solution is one in which the Ricci tensor vanishes, in the Schwarzschild 
solution it follows that covariant derivatives of the Ricci tensor vanish as well, with both Wff{ and 
W(2) thus vanishing when R^ v = 0. Moreover, this same argument immediately generalizes to general 
coordinate scalar actions based on even higher powers of the Ricci scalar or Ricci tensor. The exterior 
Schwarzschild solution is thus an exterior solution to any pure metric theory of gravity which uses any 
Ricci tensor based higher order action in place of the Einstein-Hilbert action, 30 with all of them thus 
satisfying the three classic tests. 31 Thus, in complete parallel to our earlier discussion of the difference 
between the second and higher order Poisson equations, we see that even though there is no limit in 
which any linear combination of and WftX reduces to the Einstein tensor of Eq. (J 104)) . nonetheless 
the use of any such combination as a gravitational tensor which could replace G^ v in Eq. ((TJ would 
still lead to the Schwarzschild solution. The Schwarzschild solution thus bears the same relation to the 
Einstein-Hilbert action and its higher order generalizations as the Newton potential does to the Poisson 
equation and its higher order generalizations, with use of the Einstein equations only being sufficient to 
secure the three classic tests but not at all necessary. Hence, as had first been noted by Eddington as 
long ago as only shortly after he himself led an expedition in 1919 to confirm that light did indeed bend 



30 This would not necessarily be true for any higher order Riemann tensor based action since the Schwarzschild metric 
is not Riemann flat, only Ricci flat. 

31 While the vanishing of the Ricci tensor implies the vanishing of its derivatives as well, the vanishing of the derivatives 
can be achieved without the Ricci tensor itself needing to vanish, with the vanishing of tensors such as W^X or W[tX thus 
being achievable by additional, non-Schwarzschild type solutions as well, a point we shall return to below. 
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as it passed by the sun, the success of the three classic tests does not secure the validity of the Einstein 
equations, with other options being possible. Indeed, at that time Eddington issued a challenge to the 
community to tell him which gravitational action then was the correct one to use, a challenge that the 
community has never answered even though it goes directly to the heart of gravitational theory. 

Beyond the three classic tests of general relativity, the only other test on solar system sized distance 
scales is the decay of the orbit of a binary pulsar, an effect which again only involves an exterior metric, 
this one due to the presence of two stellar sources rather than just the one involved in the solar system 
tests. The decay of the orbit is due to gravitational radiation reaction, with each star in the binary 
responding to retarded gravitational signals emitted by the other and not to instantaneous ones. Such 
retardation effects will exist in any covariant metric theory of gravity since in all of them gravitational 
information cannot be communicated faster than the speed of light, and in all of them there will thus 
be some decay of the binary pulsar orbit. While there thus will be orbit decay not only in Einstein- 
Hilbert based gravity but also in its higher order alternatives, nonetheless all such alternatives face 
the (calculationally daunting) challenge of actually obtaining the specific amount of decay which has 
actually been observed, with calculations having so far only been carried through (and with stunning 
success) is the second order Einstein theory itself. 



5.2 The Einstein static universe and the cosmological constant 

The one situation in which the full content of the Einstein equations is felt is in cosmology, since in that 
case both sides of Eq. (JTJ are non-zero, with cosmological observations being made inside the universe 
and thus in the region where is explicitly non-zero. In order to prepare for the discussion of the dark 
energy problem and the accelerating universe to be given below, it is useful to recall Einstein's effort to 
produce a static universe solution to the Einstein equations through the introduction of a cosmological 
constant. In order to see the structure of the Einstein model it is instructive to start not with a static 
model but rather with the general time dependent Robertson- Walker geometry of Eq. (J36j) and the 
generic source of Eq. (|37|) with its two separate time-dependent functions C(t) and D(t). For such a 
set-up the Einstein equations of Eq. (JTJ) yield 

3cR= -4irG(C + 3D)R (109) 

and 

+ kc * = ^ C R 2 . (110) 
3c 

A solution in which C, D and R are all constant is then obtained when 

On taking C to be positive, we see that k would need to be positive too (viz. topologically closed 
3-space), while D would have to be negative, with the Einstein static universe thus being supportable 
by a source with D/C = —1/3. Fluids in which D/C is negative are known as quintessence fluids 
and have recently come into prominence jT^j through the discovery of the accelerating universe, with 
a quintessence fluid with C + 3D < then leading to a net cosmic acceleration (R > 0) in Eq. ()109|) . 
though its role in the Einstein model is to put the acceleration precisely at the R = borderline where it 
it is neither positive nor negative. In his actual construction, Einstein himself did not use the language 
of quintessence fluids. Rather, Einstein used an energy- momentum tensor containing just an ordinary 
non-relativist ic matter fluid with a density p m and zero pressure p m , and actually modified Eq. ([TJ by 
adding a cosmological constant term to its left-hand side according to 

( ( R? V - -g^R a \ + KsT = . (112) 
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With the cosmological constant term itself having the form of a perfect Robert son- Walker geometry 
fluid with energy density p^ = cA > and pressure p^ = — cA < (so that the cosmological constant 
term itself can be thought of as being a quintessence fluid with p\ = —p\ and p\ + 3p^ < 0), we can 
relate p m and A to C and D according to C = p m /c + A, D = —A. The requirement that C and D obey 
C + 3D = then requires that p m and A be related according to the very specific fine-tuned relation 

2cA = p m . (113) 

Since on its own the ordinary p m fluid would lead to R < in Eq. ()109j) while on its own a positive 
cosmological constant would lead to R > 0, as we thus see, to get the Einstein model to be static we 
need to fine-tune the attraction generated by the energy density p m and the repulsion generated by the 
cosmological constant A so that they precisely cancel each other. Interestingly, we shall encounter a 
reflection of such a fine-tuning below in our discussion of the dark energy problem. 

With Slipher and Hubble's discovery of the expansion of the universe Einstein realized that there 
was no further need for a static cosmology, 32 and so he abandoned the cosmological constant. However, 
with the recent discovery of the acceleration of the universe (i.e. not just an expanding R > 0, but 
also an accelerating R > 0), the cosmological constant has come back into prominence, not just as a 
fundamental term of the form inserted on the left-hand side of Eq. ()112|) but also as a dynamically 
generated spontaneous symmetry breakdown contribution to the energy-momentum tensor on the right- 
hand side of Eq. (|112|) . with current cosmology requiring these two types of cosmological constant to 
cancel each other to unbelievable precision. 

From a theoretical perspective, the very fact that one is able to introduce the Ag^ u term on the 
left-hand side of Eq. (|112|) at all is a reflection of the fact that Einstein's very choice of the original 
Eq. (JIJ did not derive from a fundamental principle but only from the phenomenological need to 
recover Newton's law. 33 Indeed, it was phenomenology which lead him first to introduce and then 
subsequently to remove the cosmological constant term. It is this lack of a principle which renders the 
choice of gravitational action non-unique, with there then being no principle with which to control the 
contribution of the cosmological constant. In addition, the use of Eq. (Q) dictates that the low energy 
limit of the theory is the second order Poisson equation not just on solar system distances but on all 
distances. And so we turn now to a study of the motions of particles on larger distance scales to see 
just how well the second order Poisson equation and its Newtonian potential solution then fare. 



6 The dark matter problem 

6.1 The galactic rotation curve problem 

It is really quite remarkable how early in the study of galaxies and then clusters of galaxies that mass 
discrepancies started to appear. Specifically, not all that long after Shapley's work on the shape and 
size of the Milky Way galaxy and Hubble's demonstration that the spiral nebulae were distant galaxies 
in their own right, Oort's measurements ^H] of the velocities of stars in our local solar neighborhood 
and Zwicky's and then Smith's JH] measurements of the velocity dispersions of galaxies in clusters 
of galaxies all pointed to discrepancies. Since the discrepancy that Zwicky detected was based on the 
use of the Newtonian gravity virial theorem, Zwicky was able to conclude that unless the virial theorem 
itself did not hold for clusters (either because Newton's law did not hold on those distances or because 
the cluster had not yet virialized), there then had to be more mass in the cluster than he could detect. 
As such, it was thought at the time that any such missing mass would be in the form of astrophysical 

32 Even in his model, any departure from the 2cA = p m fine-tuning condition would lead to a non-static cosmology. 
33 Phenomcnological that is unless there is some as yet unknown principle which would actually require the validity of 
Newton's law on all distance scales, since for the moment Newton's law itself is only phenomenological. 
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objects which were too faint to be detectable, with the notion that the missing matter might be of an 
entirely different form that could never be detected optically at all coming only altogether later when 
it was realized that modern elementary particle physics could actually provide such so called Wimp 
(weakly interacting massive particle) objects. 

While there continued to be further indications of mass discrepancies over time (for a history see 
e.g. JH]); the situation started to become critical following the HI 21 cm radio studies of the spiral 
galaxies NGC 300 and M33 by Freeman [20] and of M31 by Roberts and Whitehurst [2H which found 
no sign of the expected Keplerian fall-off of rotational velocity v with distance R from the center of the 
galaxy. With hydrogen gas being distributed in spiral galaxies out to much further distances than the 
stars themselves, the HI studies were thus the ideal way to probe the outer reaches of spiral galaxies, 
and thus the ideal probe to test for the Keplerian fall-off expected of the rotation velocities of particles 
which were located far beyond the visible stars (viz. the visible mass) of the galaxies. Moreover, the 
measurement of such rotation velocities would be unhindered by projection effects if studies were made 
of spirals that were close to being edge on along our line of sight. 

The continuing persistence of the cluster problem and concern in general about mass discrepancies 
prompted detailed exploration of outer region galactic rotation curves, with early systematic surveys 
being made by Bosma [22] and Begeman j2H| as part of the University of Groningen Westerbork Synthesis 
Radio Telescope survey. What was found in these and in many related studies was that the outer 
region rotation velocities were systematically larger than the Keplerian expectation, with none giving 
any indication of a fall-off at all. In Fig. 1 we display the rotational velocity curves of a sample of eleven 
spiral galaxies which have been identified by Begeman, Broeils and Sanders [21] as being particularly 
reliable and characteristic of the general pattern of behavior of galactic rotation curves which has been 
found, with their paper giving complete data references. With the stellar surface matter density of each 
spiral galaxy in the sample behaving as 

E(R) = E e- R/Ro , 

N* = 2vr / dRRZ(R) = 2nZ R 2 , (114) 
Jo 

where Rq is the scale length and N* is the total number of stars, 34 in Fig. 1 we have plotted the 
rotation velocities (in units of km/sec) in each given galaxy as a function of R as measured in units 
of that particular galaxy's own scale length, with the rotation curves being displayed in the figure in 
order of increasing luminosity from the lowest (DDO 154) to the highest (NGC 2841) in the sample. In 
Table 1 we provide some characteristics of the galaxies, and note that the luminosities of the galaxies 
in the selected sample range over a factor of order a thousand. 



To get a quick sense of what the luminous Newtonian expectation for galaxies might look like, we 
note that for an infinitesimally thin disk of stars with distribution = £ e ~ fi//it ° and potential 

V*{r) = —j3*c 2 /r = —M Q G/r, the total centripetal acceleration is then given by [2U| 

(In Appendix A we provide a straightforward derivation of this relation.) From the behavior of the 
modified Bessel functions, the v(R) velocity profile associated with Eq. ()115|) is found to consist of an 
initial rise from R = up to a peak at R — 2.2R which is then followed by a steady fall-off. With T,(R) 

34 We assume that light traces mass, so that the surface matter density is proportional to the measured luminous surface 
brightness with a mass to light ratio which is independent of R. 
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Figure 1: Some typical measured galactic rotation velocities (in km/sec) with their quoted errors as 
plotted as a function of R/Ro where i?o is each galaxy's own optical disk scale length. The dashed 
(falling) curve shows the velocities that would be produced by the luminous Newtonian contribution 
alone. The dash-dotted (rising) curve shows the velocities that would be produced by the linear poten- 
tials of the alternate conformal gravity theory which is discussed below, and the full curve shows the 
overall predictions of that theory, with no dark matter being assumed. 
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Table 1: Characteristics of the eleven galaxy sample 



Galaxy 


Distance 


Luminosity 




(Mpc) 


(1O 9 £b ) 


DDO 154 


3.80 


0.05 


DDO 170 


12.01 


0.16 


NGC 1560 


3.00 


0.35 


NGC 3109 


1.70 


0.81 


UGC 2259 


9.80 


1.02 


NGC 6503 


5.94 


4.80 


NGC 2403 


3.25 


7.90 


NGC 3198 


9.36 


9.00 


NGC 2903 


6.40 


15.30 


NGC 7331 


14.90 


54.00 


NGC 2841 


9.50 


20.50 



R (v 2 /c 2 R) last (M/L) 
(kpc) (lO^W 1 ) {M & L- B l J 



0.48 1.51 0.71 

1.28 1.63 5.36 

1.30 2.70 2.01 

1.55 1.98 0.01 

1.33 3.85 3.62 

1.73 2.14 3.00 

2.05 3.31 1.76 

2.72 2.67 4.78 

2.02 4.86 3.15 

4.48 5.51 3.03 

2.39 7.25 8.26 



typically reaching to about four scale lengths before becoming negligible, the visible stellar material 
in the galaxy is thus contained almost entirely within this region (the optical disk), with Eq. (|115j) 
asymptoting to the Keplerian 



R W (U6) 

far outside the optical disk. To treat the luminous Newtonian contribution more accurately we need 
to allow for the fact that the optical disk actually has a thickness (with its own scale height z ), that 
the HI gas comes with its own surface brightness distribution and can have a total mass as high as 
15 per cent or so of the luminous mass, and that the two largest galaxies in the sample also have a 
central bulge. Using the calculational tools presented in Appendices A, B and C the contributions of 
all of these various components can readily be determined for each galaxy in the sample, with the exact 
luminous Newtonian expectation for v (R) reported in the typical calculation of [23] being exhibited 
as the dashed curves in Fig. 1. In generating this total luminous Newtonian contribution, the total 
number of stars N* in each galaxy was adjusted so that the luminous Newtonian contribution would 
account for the entire rotational velocity at the R = 2.2Rq peak (the maximum disk prescription), with 
the obtained values for N* being presented in Table 1 as galactic mass to light ratios M/L which thus 
express iV* in terms of the measured total luminosity of each galaxy. As we see, normalizing to this 
maximum disk prescription then completely accounts for the entire initial rise in the rotation curve, 
to therefore indicate that the inner region is well described by the luminous Newtonian contribution 
alone. However, as we also see from Fig. 1, in the outer regions of the galaxies the luminous Newtonian 
contributions totally underaccount for the measured velocities. 35 This then is the mass discrepancy 
problem in galaxies, and as such it is actually far more severe than the cluster problem, since it is based 
not on any statistical averaging over the velocities of many galaxies in a cluster or on any assumption 
regarding their virialization, but only on looking at individual orbits in galaxies one by one just as was 
done for planetary orbits in the solar system. Not only do we find such discrepancies in galaxies, we 
wryly note if the galactic data were the only data that were available to us and we had no solar system 



35 The maximal disk prescription takes the luminous Newtonian contribution to be as large as it possibly can be without 
overshooting the inner region data. Since this simultaneously fixes the normalization of the Newtonian contribution to 
the outer region, it thus represents the best that the luminous Newtonian contribution can possibly do in the outer region. 
And even if we were to change the overall normalization of N*, this would anyway not change the shape of the luminous 
Newtonian contribution, with its behavior at large R still falling off even then. 
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data at all, then from the measured galactic rotational velocities alone it would not be possible to infer 
an inverse square gravitational law at all. 36 



6.2 The dark matter solution for galaxies 

Given the prior experience with first the prediction and then the discovery of the planet Neptune, it was 
most natural to try to explain the galactic mass discrepancies by positing the presence of non-luminous 
dark matter in galaxies. In fact, there was already some theoretical basis for doing so since Ostriker and 
Peebles |26 j had suggested that there should be non-luminous spherical halos (non-luminous since there 
was no sign of any luminous such halos of the size they envisaged) surrounding unbarred spiral galaxies, 
due to the fact that in Newtonian gravity such disk shaped structures would not be gravitationally 
stable, though spherical ones would. However, any such stabilizing halo would only need to envelop the 
optical disk in order to stabilize it and would thus have no immediate need to extend beyond four disk 
scale lengths or so. Moreover, since both elliptical galaxies and clusters of galaxies are also found to have 
mass discrepancies, it would appear that whatever is causing such spherically shaped systems to possess 
these discrepancies, it could not be stability since such spherical distributions of luminous Newtonian 
matter are perfectly capable of being stable on their own. Nonetheless, the difficulty with not taking 
galactic dark matter to be spherically distributed is that if it were also to lie in a disk then it would 
not be stable either. With the measured rotational velocities far exceeding the luminous Newtonian 
contributions in the outer galactic regions, within Newtonian gravity any proposed dark matter which 
would be capable of affecting the shape of the rotation curves in those outer regions would explicitly 
have to be located in those selfsame regions itself. Moreover, with the measured velocities in the outer 
galactic regions being at least twice as big as the luminous Newtonian contributions, the dark matter 
contribution to the v 2 /R centripetal accelerations would thus have to be three to four times as big as 
the luminous one. Thus not only would dark matter have to be located in precisely the region where 
there are very few stars, the total amount of dark matter in galaxies would have to greatly exceed the 
total amount of visible matter, matter which itself is predominantly located in the inner region optical 
disk. In this sense then the galactic dark matter solution is not entirely equivalent to the Neptune 
solution to the solar system Uranus problem, since the Neptune modification to solar system orbits was 
a very minor modification to the Newtonian gravity produced by the sun and not an effect which was of 
the same order as it. From the point of view of gravity theory, a modification to the galactic luminous 
Newtonian expectation which is comparable with it in magnitude is a quite drastic modification, while 
from the point of view of observational astronomy, the notion that such a huge amount of matter is 
not only not optically observable but is not even predominantly located in the region occupied by the 
visible stars is equally drastic, to thus underscore how very serious the dark matter problem is. 

Beyond the actual assumption of the existence of galactic dark matter at all, its explicit use in 
galaxies still leaves much to be desired. Specifically, in order to actually fit the galactic rotation curve 
data explicitly using dark matter, it is conventional to take the spherical dark matter density a(r) to 
be the distribution associated with an isothermal Newtonian sphere in hydrostatic equilibrium, viz. 



as modified by the introduction of a core radius r to prevent the distribution from diverging at r = 0. 
Given such a distribution, the associated centripetal accelerations in the galactic plane are given by 



36 It is amusing to ponder how gravity theory might have evolved if the earth had been the only object in orbit around 
the sun. 




(117) 



(see e.g. Eq. 
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While Eq. ()118|) leads to rotation curves which are asymptotically flat, an actual fitting to an observed 
flat rotation curve (such as the high luminosity ones displayed in Fig. 1) is not automatically achievable 
simply by adding g^ ark to gl um of Eq. (J115)) . Rather, the two free parameters in g^ ark have to be adjusted 
galaxy by galaxy. As we see from a typical galaxy such as NGC 3198, the observed rotation curve is 
pretty flat from the inner region peak near R = 2.2Rq where (c.f. Eq. fll 1 5J) ) 



2 0AN*p*c 2 4 
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0.327rE /rV N* , (119) 



all the way out to the last detected point at R ~ 10i?o- While the asymptotic limit of Eq. ()118|) leads 
to 

v 2 daik - 47T/r C V , (120) 

there is no apparent reason why the dark matter <Jq should be related to the luminous parameters iV* 
and Rq via 

N* 

°° = was • (121) 

and yet without such a fine-tuning between the dark halo and optical disk parameters, the asymptotic 
halo contribution would not match the inner region peak velocity at all, being either larger or smaller 
than it. Moreover, even if this particular fine-tuning is invoked, at best that would only match the 
velocity at R — 2.2Rq to the velocity at r = 10i?o and not require the rotation curve to be flat at 
all points in between. To achieve such intermediate region flatness, in addition it is necessary to also 
fine-tune the halo core radius parameter ro as well (so as to match at R = 6Rq or so), and only after 
this is also done is the flatness of the rotation curve over the entire 2.2Rq < R < 10i?o region then 
secured. Unsatisfactory as this use of two fine-tuned parameters is, the situation is even worse than 
that since such a fine-tuning has to be made individually for each and every galaxy; and while a two 
parameter per galaxy prescription is basically found to work well phenomenologically for the rotation 
curves that have so far been observed, at two halo parameters per galaxy this would already lead to 22 
free halo parameters for just the eleven galaxy sample of Fig. 1 alone. 

Beyond the isothermal sphere shape for o"o, numerical iV-body cold dark matter (CDM) simulations 
of hierarchical clustering have also yielded halos, of which the typical one found by Navarro, Frenk and 
White (23 |2H] is quite close in form to the profile 

*lr)= a ° . (122) 

r[r + ro) 



Unlike the halos of Eq. (|117|) . with the halos of Eq. (|122|) being cuspy ones which diverge at r = 0, 
they make a quite substantial (though still finite) contribution to rotational velocities in the inner 
region, and have been criticized on this score |29 j since they are found to have some difficulty in fitting 
galaxies such as dwarfs and low surface brightness galaxies 37 where the inner region luminous Newtonian 
contribution is small, a region where the true nature of the halo should then be manifest. Moreover, 
even though these cuspy dark matter halos do provide good fitting to high surface brightness galaxies, 
since these particular halos would contribute in the inner region, getting good fitting in the high surface 
brightness case entails giving up the nice accounting of the inner region which the luminous Newtonian 
contribution provides, and instead requiring an interplay between dark and luminous matter even in 
that region. Very recently attempts have been made to remedy the cuspiness problem |3(J| I3T] leading 
to halos close in form to 



a[r) = aoexp 



1 



(123) 



37 



Viz. galaxies with a small N* or a small central surface brightness Eq. 
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and the issue should be regarded as being open at the present time. However, apart from issues such 
as this, there is actually a far more serious concern for essentially all dark halo models, namely even 
if one can generate dark matter halos dynamically, it is not clear what fixes which halo is to go with 
which given luminous matter distribution, i.e. if the only information available is the measured surface 
brightness of a given galaxy, from that how do we then determine which values of <To and r$ to use with 
each such given galaxy so as to be able to predict the structure of the galactic rotation curve in advance 
of its measurement, and thereby render dark matter theory falsifiable. Forms such as those of Eqs. 
(I122j) and (J123)) only give the generic shapes of halos, and do not provide any a priori determination of 
values for the do and r$ parameters or their relation to the luminous E(i?) = Soe -7 ^ profile, leaving 
do and r to be free parameters to be fit to galactic data galaxy by galaxy, to then quickly generate 
large numbers of free parameters as more and more galactic rotation curves are considered. 

6.3 Clues from the data 

In fact perhaps the most serious challenge to the whole dark matter idea is something which is intrinsic 
to the very motivation which gave rise to it in the first place. Specifically, since galactic rotation 
curves could not be fitted by luminous Newtonian matter alone, any additional matter which was to 
be invoked would then have to be non-luminous, and thus be essentially decoupled from the luminous 
matter. However, as we see from relations such as the v 4 ~ L Tully-Fisher relation which holds for the 
prominent high luminosity spiral galaxies, the velocity which is fixed by the total gravitational potential 
is in fact correlated with the luminosity, and thus whatever is producing this total gravitational potential 
has to know about the luminous matter even as it is to totally dominate over it in the potential. 
Moreover, the Tully-Fisher relation is not the only regularity in the data, and in order to get some 
guidance from the data as to what is needed of dark matter theory (or of any alternate theory for that 
matter), we now look at some particularly instructive regularities in the galactic data themselves. 

As regards the surface brightness, as first noted by Freeman, as well as obey the Tully-Fisher relation, 
high surface brightness spiral galaxies are found to all have the same central surface brightness E^, with 
all other types of spiral galaxy typically being found to have central surface brightnesses which are lower 
than this Eq 7 , with the Freeman value for E F serving as an upper limit on the general So- Moreover, 
from Eq. (jll9|) we see that Freeman limit galaxies have an inner region peak velocity which obeys 

v? um = 0.327rE F /rViV* , (124) 

with all galaxies with this common E^ then immediately obeying the Tully-Fisher relation, and doing 
so entirely because of the behavior of the luminous matter distribution alone, to thus not immediately 
point in the direction of dark matter. 

While most attention in galactic rotation curves has focussed on the flat ones since they are so very 
striking, nonetheless, as we see from Fig. 1 not all rotation curves are in fact flat. Rather, as first 
recognized by Casertano and van Gorkom [32], the data fall into three broad classes. Specifically, the 
rotation curves of the low luminosity dwarf galaxies are typically found to be rising at the last detected 
data points, the rotation curves of the intermediate to high luminosity galaxies are found to be flat, 
while the rotation curves of the very highest luminosity galaxies are found to be mildly falling. While 
there are many flat rotation curves, there are also many which are seen to be rising, and if these curves 
are to eventually flatten off, this has yet to be observed. Of these various classes, it is only the high 
luminosity, flat rotation curve galaxies which are Freeman limit, Tully-Fisher galaxies, with the lower 
luminosity, non-flat galaxies not enjoying this particular form of universality at all. 

To explore the rotation curves of these various classes of galaxies in more detail, and to try find a 
universality which both Freeman limit and sub-Freeman limit galaxies can enjoy, it is very instructive 
to focus not on the total rotational velocities, but rather on the difference between the total velocity 
and the luminous Newtonian expectation. With the luminous Newtonian contributions all falling far 
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outside the optical disk, for galaxies for which the total velocities are flat, the difference between the 
total and the luminous Newtonian expectation can not itself flat, but must instead actually be rising in 
the detected region. Inspection of the flat rotation curves in the sample of Fig. 1 shows no galaxy for 
which the difference between the total and the luminous Newtonian expectation is itself flat, so if the 
effect of galactic halos is to produce contributions to rotational velocities which are to be asymptotically 
flat (as per Eq. (|12U|) ). no sign of any such asymptotic behavior is yet manifest in the data. A similar 
shortfall pattern is also found to obtain for the dwarf galaxies, galaxies whose rotation curves are not flat, 
with the shortfall between the measured rotational velocities and the luminous Newtonian expectation 
again being seen to be rising. In fact the situation for dwarfs is actually more severe than for bright 
spirals, since while the ratio of the total velocity to the luminous Newtonian expectation at the last 
detected data point is a factor of two for the bright spirals, for dwarfs the ratio is a factor of three (and 
thus a factor of nine in the centripetal acceleration). Mass discrepancies in dwarfs are thus even more 
pronounced than in the bright spirals. 

As we see, the characteristic feature of all the eleven galaxies in the sample is that the shortfall 
between the measured rotational velocities and the luminous Newtonian expectation is a shortfall which 
explicitly rises with distance from the center of each galaxy. It is thus of interest to see whether this rise 
might have any universal structure which is common to all the galaxies in the sample. To this end it is 
very instructive to evaluate the centripetal acceleration at the last detected data point in each galaxy, 
with the obtained values being exhibited in Table 1. As we see from the table, the quantity (w 2 /c 2 i?)i ast 
is found to basically increase with luminosity. With the asymptotic luminous Newtonian contribution 
to the centripetal acceleration being given by Eq. (jll6|) where N* has been fixed once and for all by 
the inner region maximum disk procedure, its subtraction from (f 2 /c 2 -R)i ast is found to yield a function 
which satisfies a two parameter formula containing one term which is linear in iV* and a second term 
which, quite extraordinarily, is independent of it altogether, with the full (v 2 /c 2 -R)i ast and tf ast being 
given by 

(v*_\ _ (3*N* | 7* N* | to 
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vL t = + + (125) 

where 7* and 70 are two new universal parameters with values 

7* = 5.42 x KT^cirT 1 , 7o = 3.06 x lCT^cm^ 1 . (126) 

To get a sense of the relative importance of these various terms, we note that with the solar (3* being 
given by (3* = 1.48 x 10 5 cm, the quantity ((3* /7*) 1 / 2 evaluates to ((3* /7*) 1 / 2 = .52 x 10 23 cm, with 
the f3*N*/R 2 and / y*N*/2 terms thus being competitive with each other on none other than galactic 
distance scales. With the large galaxies typically having of order 10 11 stars, for them the , y*N*/2 and 
7o/2 are of the same order of magnitude, with all three contributions in Eq. (|125|) being of importance, 
and with the falling luminous Newtonian contribution having its most significant effect in the most 
luminous (largest iV*) of the galaxies. However, for small N* the 7oC 2 i?/2 term completely dominates 
in vf^ with the last data points then rising with R, just as exhibited in Fig. 1 for the dwarf galaxies, 
and with Eq. (jl25|) precisely capturing the essence of the pattern identified in [3*2"] . 

In order to assess the significance of Eq. (j!25J) . we note that there is nothing at all special about the 
last data point in each galaxy. Such data points are fixed purely by the observational limit of the radio 
detectors used to survey the galaxies, a limit which is not fixed by the galaxies themselves but by their 
distances from us. And yet we are obtaining a universal formula. It appears to us that the only way 
to achieve such an outcome is to have shortfalls which are rising not just at the last currently detected 
R but at all possible R beyond the optical disk, so that no matter where an individual detector just 
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happens to cut off, it wiii be doing so on a curve that is universally rising, to thus always yield the same 
Eq. (|125[) . Since the eleven galaxy sample is large enough to be characteristic of galactic rotation curve 
systematics and since it was by no means chosen by the authors of [21] in order to obtain a formula 
such as that of Eq. (|125j) . we thus believe we can identify Eq. (|125|) as a universal characteristic of 
galactic rotation curve systematics, one which is of interest in its own right and which dark matter 
theory is therefore required to reproduce. With Eq. (j!25J) reducing to the iV*-independent relation 
v \ast = 7oc 2 -R/2 in the small N* limit, the small N* limit exposes the departure from the luminous 
Newtonian expectation in the starkest terms. And it could be that ACDM (CDM with a cosmological 
constant A) generated dark matter halo models are having difficulty fitting dwarfs because relations 
such as v 2 ast = 7oc 2 .R/2 might not be natural to them. 

In regard to Eq. (jl25J) . we additionally note that the quantity which measures the departure of the 
centripetal acceleration (f 2 /-R)i as t from the luminous Newtonian expectation is a universal acceleration 
7 c 2 /2 (7o + 7*A r * is never smaller than 7 ) whose magnitude is given by 7oC 2 /2 = 1.4 x 1CT 9 cm/sec/sec, 
a value which is close to that found for the universal acceleration parameter ao of the MOND theory 
discussed above. Further support for a role for a universal acceleration was noted by McGaugh [3*3*] 
who studied the behavior of the quantity Md yn (-R) / M\ um (R) as a function of the measured centripetal 
acceleration v 2 (R)/R at points R within galaxies. Here M dyn (R) = Rv 2 (R)/G is the total amount of 
matter interior to R as would be required by Newtonian gravity to produce the measured v 2 (R)/R, 
while Mi um (i?) is the total amount of luminous matter which is detected within the same region. In his 
study of a quite extensive sample of galaxies McGaugh found that mass discrepancies [viz. M lum (i?) 
less than the needed M dyn (R)} systematically occurred in galaxies whenever the measured v 2 (R)/R 
fell below a universal value of order 10~ 8 cm/sec/sec or so, a value which is in accord with that of 
both the acceleration ao associated with MOND and the 7oC 2 /2 acceleration associated with Eq. (|125[) . 
Thus from both McGaugh's phenomenological study and from the phenomenological Eq. (j!25j) . we see 
that it is indeed an acceleration scale (or equivalently an inverse distance scale times the square of the 
speed of light) which determines exactly when luminous Newtonian matter fails to account for observed 
data. As such, the acceleration criterion which Milgrom originally identified has to be regarded as 
being a bona fide property of nature. As formulated, the MOND theory actually possesses not one 
but two logically independent components. It comes not just with this acceleration criterion which 
indicates when the luminous Newtonian prediction is to fail, but also it provides a model (Eq. ([<""*"])) for 
how Newtonian gravity is to then be modified in the region (the MOND region) where the luminous 
Newtonian expectation does in fact fail. However, regardless of how the theory is to behave in the 
MOND region, i.e. regardless of what particular form the function u(f /a ) might take, 38 and regardless 
of whether one likes any particular form (such as the one of Eq. ([1(J1|) with its square root structure) 
that might be used to fit data, it remains true that Milgrom's acceleration criterion is the criterion which 
determines when departures from the luminous Newtonian expectation are to first occur, a criterion 
which fundamental theory is therefore obliged to explain. As we thus see, from Milgrom's work, from 
McGaugh's work and from Eq. (|125[) . there is a lot of universality in the rotation curve data, universality 
which does not immediately appear to be natural to dark matter and which therefore challenges it, a 
universality which we will use below to guide us when we consider alternatives to dark matter. 



38 If in Eq. I|125|) we ignore the j*N* /2 term and set ao = 7oc 2 /2, we can then write Eq. I|125|) in the generic form of 
Eq. (|99|l with the function v(f/ao) being given by u(f /qq) = 1 + ao//, an issue we will elaborate on further below. 
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7 The dark energy problem 



7.1 The cosmological dark matter and flatness problems 

As we go to scales even larger than galaxies, mass discrepancies are not only found to persist, two new 
types of mass discrepancy begin to appear, one for clusters of galaxies and the other for the entire 
universe itself. While non-relativistic luminous Newtonian expectation shortfalls are found not just for 
galaxies but for clusters of galaxies as well, for such clusters general-relativistic shortfalls are also found, 
with the use of the Schwarzschild metric for a cluster leading to far less gravitational lensing by the 
cluster than is observed if only the luminous inventory of the cluster is used as the Schwarzschild source. 
Thus for clusters one can say that either the dark matter idea is internally consistent since both non- 
relativistic and general-relativistic shortfalls are found in one and the same system, or one can say that 
absent dark matter, with clusters it is not just Newtonian gravity which is failing but Einstein gravity 
too, with the cluster geometry probed by light from a distant quasar not being the Schwarzschild metric 
associated with the visible sources. In making a statement such is this it is important to note that the 
detection of gravitational lensing is generally regarded as being a great triumph for Einstein gravity. 
However the fact of lensing as opposed to the amount of lensing only requires that gravity be a metric 
theory to which light couples. It is the amount of lensing which is sensitive to the specific make-up of 
the lens and the specific geometry which it produces, and unless there is dark matter in clusters, then 
the standard theory would simply not be predicting enough. As with our discussion of ACDM generated 
galactic halos, what is needed to make cluster dark matter into a falsifiable theory is a prediction of the 
amount of lensing given only a knowledge of the luminous content of the cluster. Currently the amount 
of dark matter in clusters is inferred only after the lensing measurements have been made. While 
making lensing predictions based only on the known luminous content of a cluster thus serves as a goal 
for dark matter, at the present time it also remains an objective for the alternate gravitational theories 
we shall discuss below, including those for which galactic rotation curve predictions can be made from 
luminous information alone. Progress on cluster gravitational lensing should thus be definitive for both 
dark matter and its potential alternatives. 

The other large distance scale on which mass discrepancies are found is the largest distance scale 
there is, viz. that of the universe itself, with the standard theory not only again encountering the 
issue of dark matter, this time the problem is compounded by a need for dark energy as well. To 
heuristically see why such problems might be anticipated, it is instructive to consider the implications 
for a simple Newtonian cosmology of the existence of a universal acceleration. Specifically, we recall 
that in a critical density Newtonian cosmology of size the Hubble radius L = c/H (viz. a cosmology 
where pc = 3Mc 2 /4:7iL 3 = 3H 2 c 2 /8nG), the Hubble radius happens to be equal to the Schwarzschild 
radius of the matter within it (viz. 2MG/c 2 = L). For such a cosmology the centripetal acceleration of 
a test particle located at the Hubble radius is given by v 2 /L = MG/L 2 = cH/2. Using a value of 72 ±8 
km/sec/Mpc for H [33], we find that cH/2 evaluates to cH/2 = 3.5 x 10~ 8 cm/sec/sec, and with this 
value being very close to the values found for both the MOND ao and the acceleration 7oc 2 /2 associated 
with Eq. (|125|) . we see that the centripetal accelerations of Newtonian cosmology occur precisely in the 
acceleration regime where luminous Newtonian matter falls short for galaxies. We this anticipate, and 
shall in fact shortly see, that a pure luminous matter based standard Einstein gravity cosmology will 
not adequately fit available cosmic data either. Moreover, in addition this time we will find that the 
situation can not even be rectified by invoking dark matter, with it being this very failure which has 
led to the introduction of yet another ad hoc dark object, viz. dark energy. 

In standard Einstein gravity a full general-relativistic treatment of cosmology is based on the Fried- 
mann equations of Eqs. (ITU3D . ffTTTT^ and flU 



R 2 + kc 2 



8ttG 



CR 2 , 3R 




(127) 



3c 



36 



as written for the moment in terms of the generic source T^ u = [C(t) + D(t)]U^U u + D(t)g^ u given 
in Eq. (J37j) . Given the Robert son- Walker metric of Eq. (J36|) . the objective of cosmology is to then 
determine the parameters k and R(t) which appear in the metric and determine the matter content 
of the universe as expressed through C(t) and D(t). As a result of many years of observational effort, 
values for the Hubble parameter R/R = H have now converged on an agreed value for H as given 
above. With H and G being known, it is convenient to introduce the dimensionless quantities 

= 8*GC(t) = _kJ_ 

cv ) 3cH 2 , fcV ) R2 ^ , V ) 

and the deceleration parameter 

. „ RR 

= — ^ > ( 129 ) 
since in terms of them the Friedmann equations of Eq. (|127j) can be written very compactly as 

r 2 + kc 2 = R 2 n c (t) , n c (t) + fifc(t) = i , q(t) = ^(i + ^p\nc(t) . (130) 

Historically, the cosmological source was taken to be ordinary matter with C = p m /c, D = p m /c, 
with equation of state p m = p m /3 for radiation and high energy (kT ^> mc 2 ) matter and equation of 
state p m = for low energy (kT <C mc 2 ) matter, with an inventory of the visible universe using galaxy 
counts leading to a current era value for 

n-W - (13D 

of order 0.01 at the current time t . 39 Such a value for fl M (t ) is quite close to the critical value M = 1, 
the value at which the Friedmann equation would require the parameter k to be zero (flat universe). 40 
As first noted by Dicke, this closeness of Qm^o) to one creates a severe problem for the Friedmann 
equation, the so-called flatness problem. Specifically, regardless of the value of k, the use of the above 
normal matter equations of state in the Friedmann equation leads to the occurrence of a big bang 
singularity in the early universe, with R(t = 0) diverging at t = 0. Near to such a singularity the same 
Friedmann equation then requires that M (i = 0) be incredibly close to one, again regardless of the 
choice of k. Then, with Qm (t = 0) being close to one in the early universe, and with the universe then 
expanding for an entire Hubble time (1/H(t ) ~ 10 10 years), it is very hard to understand why the 
current era (to) would be anywhere near unity today instead of having redshifted to close to zero by 
now. In fact for the Friedmann equation to lead to a cosmic evolution in which the current era value 
of Q M (t ) would be close to one today requires fine-tuning the initial value of fi M (t = 0) at a typical 
early universe time to be one to within one part in 10 60 or so. As an evolution equation, the Friedmann 
equation is highly unstable, with it not being able to easily reconcile its initial singularity with the 
current era value of f^M^o)- More naturally, it would have ^m(^) fall to zero in the early universe itself. 



7.2 Fine-tuning and the inflationary universe 

A beautiful solution to this conflict was proposed by Guth (S3], the so-called inflationary universe, in 
which a de Sitter phase was to precede the current Robertson- Walker phase, with the rapid exponential 
expansion of the de Sitter phase precisely providing a dynamics which would lead to the needed initial 
condition for the subsequent Robertson- Walker phase, no matter what the value of k. An additional 

39 Since such inventorying does not undercount, 0.01 thus serves as a lower bound on Qjvf (to)- 

40 It is called the critical value since the sign of the 3-curvature k of the universe has the same sign as fijif (io) — 1- 
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virtue of having such an early universe de Sitter phase was that it also provided a solution to the horizon 
problem associated with the observed uniformity of the cosmic microwave background across the sky, 
with the de Sitter phase phase being able to causally connect regions which a luminous (f2 A /(t ) = 0.01) 
or even a luminous plus dark (Qm(^o) > 0.01) matter Friedmann universe could not. 

Geometrically it is possible to describe a de Sitter phase using the comoving Robertson- Walker 
coordinate system, and for it we can take the Friedmann equation to be of the form associated with a 
cosmological constant source, viz. 

R 2 {t) + kc 2 = ac 2 R 2 {t) , (132) 

where a = 8nGA/3c 3 . For the allowable choices of sign for a and k solutions to Eq. ()132|) can readily 
be found, and they are given as 



R(t, a < 0, k < 0) = I - I sin[(-a) 1/2 ct] 



a 



R(t,a = 0, k < 0) = (-k) 1/2 ct , 

( k\ 1 ' 2 

R{t, a > 0, k < 0) = sinh(a 1/2 ct) 



R(t, a > 0, k = 0) = R(t = 0)exp(a 1/2 ct) 



ot . 

1/2, 



On defining 



R(t, a > 0, k > 0) = \^-J cosh(a i/2 ct) . (133) 
„ . . ac 2 8ttGA 

we find that in the de Sitter phase the quantities q(t) and Q A (t) are related by q{t) = —Q A (t), with q{t) 
and Q A (t) then being found to be given by 

tt A (t,a < 0, k < 0)) = -q(t, a < 0, k < 0) = - tan 2 [(-ct) 1/2 ct] , 

Q A (t,a = 0, k < 0)) = -q(t, a — 0,k<0) = , 

tt A (t,a > 0, k < 0)) = -q(t, a > 0, k < 0) = tanh 2 (a 1/2 ct) , 

Q A (t, a > 0, k = 0)) = -g(t,a > 0, fc = 0) = 1 , 

fi A (t, a > 0, fc > 0)) = -g(t, a > 0, k > 0) = coth 2 (a 1/2 ct) (135) 

in the various cases. As we see from the solutions, when a is taken to be positive, then no matter 
what the value of k, at late times the expansion radius will grow as exp(a 1 / 2 ct) with ^a(^) asymptoting 
to one and ^(t) = 1 — ^a(^) to zero. What a period of de Sitter inflation thus achieves is that it 
quenches the contribution of curvature to cosmic expansion, with the universe inflating so much that 
curvature is no longer able to contribute to the expansion. If at the end of the inflationary era all of 
the energy density contained in Q A (t) can be converted into a matter QM{t), Q A {t) will then play no 
further role and flM{t) will then enter the Robertson- Walker era with a value that will be extremely 
close to one in a universe in which curvature is no longer capable of influencing cosmic expansion. With 
Qk{t) thus remaining negligibly small throughout the Robertson- Walker era, and with ^a(^) have been 
disposed of, ^Ai(t) will then continue to remain close to one, and will thus still be extremely close to 
one today. The inflationary solution to the flatness problem is thus to make an ^Af(^) which is of order 
one actually be equal to one (in each and every Robert son- Walker epoch), with the shortfall between 
a luminous contribution of 0.01 and one itself then requiring a 0.99 contribution through matter that 



L That a de Sitter phase would solve the horizon problem had also been noted in [201 and [77). 
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will contribute to fU/(£) with exactly the same equation of state as ordinary matter but whose visual 
absence would require it to be non-luminous. Inflation thus leads to the need for cosmological dark 
matter in an amount which would make the current era ^m(^o) be equal to one, viz. to precisely the 
critical density value. With such a value for ^m(^o), and with D/C being zero for the non-relativistic 
fluids (dark or luminous) of relevance at the three degree current era temperature, inflation would thus 
entail that the current era q(to) would be equal to 1/2 (c.f. Eq. ([130)) ). so that even while inflation itself 
leads to a negative q(t) (c.f. Eq. ([135)1 ) in the inflationary era, following the transfer of energy density 
from Qh(t) to Oji^(t) at its end, we transit into a Robertson- Walker phase in which q(t) is positive. 
While inflation thus provides a nice solution to both the flatness and the horizon problems, the very 
fact that it produces an era of exponential expansion no matter what the value of k means that inflation 
inflates away any dependence on k, with inflation thus not being sensitive to the sign or the magnitude 
of k. The fact that O^(t ) is to be negligible today does not mean that k itself is zero. Rather inflation 
enables the current era fifc(£o) to be negligible even if k is non-zero, though this very insensitivity to k 
means that the physics of what is to actually fix k is thus not addressed by inflation. Since the universe 
does have a global topology (flat, open or closed), something does have to fix it, with physics either 
beyond inflation (or perhaps instead of it) still being required, a point we return to below. 

7.3 The accelerating universe and the cosmological dark energy solution 

While the development of the inflationary universe model immediately prompted a vigorous search 
for dark matter candidates, a search which so far has only yielded massive neutrinos (though with 
masses which are too small to be of significance for the needed CIm^o)), curiously, at no time did any 
astrophysical observations themselves ever indicate that f^M^o) was i n fact equal to one or that q(t ) 
was equal to 1/2. In fact the best estimates for f^A/(£o) as obtained from using dark matter to fit 
the luminous shortfalls found for clusters of galaxies, lead to an Qm^o) of no more than 0.3 As 
regards determining q(to), we note that since the deceleration parameter involves a second derivative 
with respect to time, its measurement requires a measurement of the change in the first derivative 
(the Hubble parameter) over cosmologically separated distances, and thus requires an extension of the 
Hubble plot out to very high redshift. A sufficient such extension out to z — 1 has recently been made 
using type la supernovae as standard candles, and has led (TJ [2f to an epochal discovery, namely that 
far from being equal to 1/2, the current era deceleration parameter is not even positive at all, with the 
current universe actually being in an accelerating rather than a decelerating phase. However, before 
addressing this remarkable finding of the supernovae data, we note first that a fit to the 54 most reliable 
of the original 60 type la supernovae sample of [2] using the luminosity function 



associated with a pressureless f2 M (t ) = 1, &>k(to) — universe is found to yield a \ 2 value of 92. 9, 43 
a value which thereby actually renders an f2jv/(t ) = 1, ^fc(^o) = universe unviable. Consequently, 
recent supernovae data rule out the ^m(^o) = 1 inflationary universe solution to the flatness problem. 

To get a sense of what is required to fit the supernovae data, we note that as the matter content 
is reduced in this same pressureless fW(^o)> fifc(to) universe, a reduction which requires a concomitant 
increase in Q k (t ), the associated \ 2 is found to decrease, reaching a value of x 2 — 61.5 when Q M (t ) 
reaches its minimum value of zero, a situation in which Qk{to) has then risen to one (viz. k now 
negative), the deceleration parameter has dropped to zero, and the luminosity function can be written 

42 That a positive cosmological constant a leads to a negative q(t) is something we shall return to below. 
43 The various x 2 values reported here and below are taken from [351159) . 
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in the closed form 
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Since this particular model is very close to a universe with just regular luminous matter (f2 M (t ) = 0.01) 
and no other components, we see that while this model still gives a x 2 which is somewhat larger than the 
number of independent degrees of freedom, the model does not fare all that badly, suggesting that it may 
have captured some of the essence of the Hubble plot data. In fact paradoxically, a pure luminous model 
with Q M (t ) = 0.01, flk{t ) = 0.99 actually does better with the supernovae data than a k < luminous 
plus dark matter model with f^M^o) = 0.3, flkito) = 0.7. 44 Within a pressureless Qm^o), ^k{to) model, 
the only way to reduce the \ 2 even lower than its Qm^o) = 0, flkito) = 1 value is to make either G or p m 
negative - not that one of course can do so in the model, though it does show the needed trend, with the 
data favoring an even more reduced deceleration parameter, viz. one which would actually be negative, 
a trend which having large amounts of dark matter actually works against. As it stands then, we see 
that the supernovae data exclude the only value of fi^^o) f°r which the Friedmann evolution equation 
has so far been shown to not have a fine-tuning problem (viz. f2jvf(£ ) = 1), and leaves us with a universe 
which has even more cosmic repulsion than that achievable by the smallest permissible f^M^o)- We 
stress this point now to show that even before invoking the cosmological constant (to get some additional 
cosmic repulsion akin to that which it provided in the original Einstein static universe discussed earlier), 
once f^M^o) is not equal to one we are right back in the flatness problem, with the Friedmann equation 
again requiring fine-tuning. Thus even while the introduction of the cosmological constant will also 
entail fine-tuning, a fine-tuning of its value down from that anticipated from fundamental physics, 
this particular fine-tuning is quite distinct from that associated with the Friedmann cosmic evolution 
equation, with the standard cosmological evolution equation already having fine-tuning problems even 
without a cosmological constant, fine-tuning problems which do not disappear following its introduction. 

As regards the possible addition of a cosmological constant to the model, an addition in which Eq. 
(I130j) is then replaced by 



we now see that a positive ^(i) will nicely serve to reduce q{t) just as we would want. To get a sense 
of how big we would need the VL A {t) contribution to be, we note that for a pure Vt A {t) = 1, Vt M {t) = 0, 
&k{t) = universe the luminosity function can be written in the closed form 



with it leading to x 2 = 75.8 for the same 54 supernovae data points. While this value is still too large 
to be acceptable, it is closer to the number of degrees of freedom than the \ 2 associated with a pure 
(flM{t), ^a(^); Qk{t)) — (1)0,0) universe. Hence within the family of universes with Qk{to) = 0, we 
can anticipate a best fit which is weighted more to fl\(t) than to Q M (t), with the values Om^o) — 0.3, 
^a(^o) = 0.7 being found to yield \ 2 — 57.7 and the explicit data fit displayed as the lower curve in Fig. 
2. With these values of Qm(^o) an d ^a(^o) the deceleration parameter given by Eq. (|138|) evaluates to 
<?(£()) = —0.55, and is thus seen to be quite negative. 

44 Of course both of these models will still have to be excluded since neither of them can account for the uniformity of 
the cosmic microwave background, and both of them suffer from the flatness problem. 
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7.4 The freedom inherent in fitting the current supernovae data 



To get a sense of the fitting freedom inherent in the supernovae data, it is instructive to look not just at 
the curvature-free Ofc(t) = family of solutions, but also at a particular class of matter-free Om(^) = 
solutions, viz. the a > 0, k < ones as given in Eqs. (j!33|) and (jl35j) with R(t, a > 0, k < 0) = 
(-/c/a) 1/2 sinh(a 1/2 c)f:) and Q\(t, a > 0, k < 0) = -q(t,a > 0, k < 0) = tanh 2 (a 1/2 ct). This particular 
family of solutions is one which is accelerating in all epochs [q(t, a > 0, k < 0) is negative at all t for 
any positive a], and on parameterizing a in terms of the current era deceleration parameter q(to) = qo, 
its luminosity function is given [3SJ EH] as the one parameter 
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Using this luminosity function a best fit of \ 2 = 58.62 is obtained with q = —0.37, a decelera- 
tion parameter which again is negative, and is exhibited as the upper curve in Fig. 2. 45 With the 
(n M (t ),^A(to),^k(t )) = (0.3,0.7,0.0) and the (0 M (^o), A (£o), ^fe(^o)) = (0.0,0.37,0.63) fits being of 
equal and indistinguishable quality, the best fits in the (Qm(^o), ^a(^o)) plane will to good approxima- 
tion thus lie along the line 

n A {t ) = l.m M {t ) + 0.37 , fi fc (to) = 0.63-2.m A/ (t ) (141) 

which runs through these particular two fits. Along this line we obtain 

q(t ) = -0m M (t )- 0.37 , (142) 

with the current era deceleration parameter thus being negative everywhere along the best fit line. 
With ^a(^o) being zero on the best fit line when Qm^o) is given by Qm^o) = —0.34, we recover our 
earlier observation that in pure matter universes with no cosmological constant, the best fits are actually 
obtained with fl M (t ) negative. 



7.5 The cosmic coincidence problem 

From the point of view of the standard model, there will be a problem no matter where on the best fit 
line it finds itself, since along the best fit line Om^o) and ^a(^o) are of the same order of magnitude. 
In the presence of a cosmological constant, the occurrence of an early universe initial singularity in the 
Friedmann equation of Eq. (J 138)) this time translates into the requirement that at the initial time the 
quantities Qm(^o) and ^a(^o) must add up to one, viz. 

fijif (t = 0) + tt A (t = 0) = 1 . (143) 

However, since p m redshifts while A does not, in order for these two quantities to be of the same order 
today, in the initial universe Om(£ = 0) must have been incredibly close to one (typically to within one 
part in 10 60 or so) and fl\{t = 0) must have been of order 10~ 60 . There is thus again a fine-tuning 
problem, with the cosmic evolution Friedmann equation again not being able to naturally reconcile its 
initial singularity with the current era values of ^m(^o) and Q\(to). With the most natural expectation 
for a Friedmann universe with both matter and a cosmological constant being that entirely in the early 
universe itself f2^(i) should go from one to zero and fl\{t) should go from zero to one, the closeness of 
their values some ten billion years later is often referred to as the cosmic coincidence problem. 

45 In Fig. 2 we have labelled this fit as the conformal gravity fit, as it will be associated with the conformal gravity 
alternative to standard gravity which will be discussed below. However, for the moment it need only be regarded as just 
one of the possible phenomenological fits obtainable by treating f^M^o) and f^A^o) as free parameters in Eq. (|138|l . 
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Figure 2: The go — —0.37 conformal gravity fit (upper curve) and the ^jw(^o) — 0.3, I^a^o) — 0.7 
standard model fit (lower curve) to the z < 1 supernovae Hubble plot data. 
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While inflation does not solve this cosmic coincidence problem, an early universe inflationary era 
can nonetheless still occur (and indeed its occurrence would still solve the horizon problem), with 
then still being negligible throughout the Robertson- Walker era. Within the class of best fit solutions 
to the supernovae data this would then single out the fUf(to) — 0.3, fi\(£o) — 0.7 one, and thus bring 
inflation into agreement with the ^m(^o) — 0.3 dark matter cluster estimate. Further support for 
this fUf^o) — 0.3, fU(£o) = 0.7 realization of inflation has recently been obtained from study of the 
temperature anisotropy of the cosmic microwave background, since such anisotropics can actually be 
generated in a period of rapid early universe inflation in which initially microscopic fluctuations could 
be amplified to macroscopic size. The size to which such inflationary fluctuations will grow at the time 
of recombination will thus act as a theoretical standard yardstick which will be imprinted on the cosmic 
microwave background, with the apparent size with which it will then appear to us then reflecting the 
geometry through which photons have had to travel to us since recombination. Through use of this 
technique the cosmic microwave background anisotropy measurements are then found to have a best fit 
in the (Om(^o)i ^a(^o)) plane which scatters around the line [HE] 



with the compatibility of Eq. (jl41|) with Eq. (|144|) remarkably bringing us right back to ^m(^o) = 0.3, 
^a(^o) — 0.7. Consequently, these values are now regarded as being the ones both required of and 
established by the standard cosmological theory, with the standard theory thus having achieved its 
primary objective of determining the matter content of the universe. 

However, even while the standard theory has now ostensibly done its job, and while one should 
not understate how remarkable it is that there actually is a specific single set of Q M (t ), Q\{t ) values 
which fits currently available cosmological data at all, 46 the universe which emerges is still a highly 
perplexing one which not only possesses huge amounts of still poorly understood dark matter, it needs 
to augment this dark matter with an even more poorly understood f2\.(£o)- However, before looking at 
the implications of such a value for f^^o), it is important to note that there could be an alternative 
to it within standard gravity, namely the so-called quintessence model of ^H]- The quintessence model 
was actually developed prior to the discovery of the accelerating universe, with its objective at the time 
being to try to reconcile an fW(to) which would be less than one with inflation in a way that would 
not require any early universe fine-tuning of the Friedmann equation. To this end a quintessence source 
was added to the dark matter source, one which would then have an Qq (to) which would be equal to 
1 — Q M (t ) today, and thus be of the same order as Q M (t ) today, to initially appear to replace the 
^m(^o); ^a(^o) cosmic coincidence by an Q M (t ), £lcj(t ) cosmic coincidence instead. However, it was 
found possible to construct a tracking potential for the quintessence fluid which would have the effect 
of driving £)q(£) and ^m(^) to comparable current era values today regardless of what early universe 
initial conditions were presupposed, to thus render the cosmic evolution equation free of early universe 
fine-tuning problems. In order for the quintessence fluid to be able to achieve this objective, it would 
have to have an equation of state in which p/p = w is negative, and in order for it to have escaped 
visual detection it would need to be non-luminous. Non-luminous negative pressure fluids such as these 
are known as dark energy, 47 with a cosmological constant being a special case of such fluids, one with 
w = — 1. With quintessence fluids having negative w, they can thus lead to cosmic repulsion, and if w is 
negative enough can provide for fitting to the supernovae data of the same quality as the ^m(^o) = 0.3, 
^a(^o) = 0.7 cosmology fit shown in Fig. 1. While quintessence is thus offered as a candidate solution 
to the cosmic coincidence problem, it leaves open the question of what is to happen to the cosmological 

46 While it might be thought that the lines of Eqs. I|141|) an d H 144(1 have to cross somewhere, this is not quite the 
case since the lines are not actually of infinite extent, with both of them only describing their respective data in limited 
regions, regions which remarkably do in fact overlap. 

47 Since it is their negative pressure which is abnormal and not their energy density (which remains positive), a possibly 
better designation for quintessence fluids would be dark pressure. 
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constant, with quintessence only being able to succeed if some mechanism can be found which would 
naturally set the cosmological constant to zero. For quintessence models to be valid then still requires 
a solution to the cosmological constant problem. 

7.6 The cosmological constant problem 

Indeed the real problem with the cosmological constant is that rather than provide an ^a(^o) as small 
as 0.7 or even one equal to zero, the a priori expectation for fl\(t ) is that it will be huge, being 
as large as 10 120 if fixed by the quantum gravity Planck scale, and being as large as 10 60 or so if 
fixed by the electroweak phase transition scale, with use of either of these values leading to absolute 
phenomenological disaster for the standard theory. Moreover, even if Q\(t) is quenched to zero at the 
end of the inflationary era, as the universe then cools down it will go through phase transitions such as 
the electroweak one at around T v = 10 15 degrees, at which point a brand new cosmological constant 
will be generated by the change in free energy due to the phase transition, with the difference between 
the free energy just above the phase transition and its value at current temperatures way below it 
being of order aT v where a is Stefan's constant. Not only is such a aT v huge, since the free energy is 
lowered in a phase transition, one would even expect it to produce a negative contribution to A, with 
A then differing both in sign and magnitude from the required scale of the order of few degrees which 
is associated with the positive £l\(to) = 0.7. Apart from the electroweak scale, there might even be 
a supersymmetry breaking scale as well, and while supersymmetry could provide some possible dark 
matter wimp candidates, the very same mechanism which would put their masses into at least the 
TeV region (the current lower mass bound associated with their lack of detection to date) would at 
the same time generate an even bigger contribution to the cosmological constant than that associated 
with electroweak symmetry breaking itself. Thus mechanisms which will generate the cold dark matter 
particles needed for Qm^o) = 0.3 will also generate huge contributions to Q\(t ) in a dynamics which 
would leave us nowhere near the cosmic coincidence values at all. This then is the cosmological constant 
problem. 

At the present time within the standard model no convincing resolution of the cosmological constant 
problem has been identified, and while there are some suggestions, none of them are anywhere near being 
at the stage where they could be falsified. Enormous efforts have been made to find a dynamics which 
would naturally quench Q\(to) down from either 10 60 or 10 120 to a current value of 0.7, so far to no avail. 
Additionally, one can consider the possibility that in addition to the overall dynamical cosmological 
constant which is generated on the right-hand side of the Einstein equations as the matter density 
goes through its various phase transitions, there could also exist an a priori fundamental cosmological 
constant on the left-hand side, one which will then almost, but not quite cancel the dynamical one so 
that the current era observer then sees a net Qx(t ) = 0.7. Such a proposal requires that a fundamental 
cosmological constant exist as part and parcel of the fundamental gravitational equations, that it be 
included in gravitational theory in all epochs, and that it be preset with just the value needed to lead 
to a net fl\{to) = 0.7 today. It is not clear how one might test such an idea, with such a cancellation 
having to be regarded as being a cosmic coincidence itself. 

Another proposed mechanism is the so-called anthropic principle which takes the view that there 
exists an infinite number of universes in which A takes on all possible values, with observers such 
as ourselves only being able to exist and observe the universe in those particular universes in which 
the value of A permits the emergence of intelligent beings. In a sense such a viewpoint is actually a 
denial of physics, saying that when we can explain something we do, and when we cannot, we appeal 
to the anthropic principle. 48 From a more technical viewpoint it may be the case that even with an 
infinite number of universes, the cosmological constant could be Planck density in all of them, and that 

48 A possibly more accurate name might be the anthropic lack of principle, with it not being clear which phenomena 
are to be attributed to the anthropic principle and which not. 
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while the Planck length itself might vary from universe to universe, the relation of A to it might not. 
Additionally, it would be an extraordinary coincidence if the only possible allowed range of A needed 
for intelligent beings would precisely be within the range needed to solve the early universe Friedmann 
equation fine-tuning problem. Rather, one might expect a somewhat broader range for A, and with the 
Friedmann fine-tuning required range (viz. to one part in 10 60 ) then being a very small fraction of it, 
it is far more likely that the anthropic principle would put us outside the Friedmann required range 
rather than in it. Finally, if the anthropic principle is to be a part of physics, then there has to be some 
falsifiable prediction that it makes which would enable us to test it. 

7.7 Testing dark energy beyond a redshift of one 

However, regardless of any of these theoretical concerns, there is actually a direct, model- independent 
test that one can make to see whether it really is true that the real world cA is indeed equal to 
0.7 p c ~ 0.6 x 10- 8 ergs/c.c. (viz. a = Vt K H 2 /c 2 = 8ttGA/3c 3 ~ 0.4 x 10- 56 /cm 2 ), this being the value 
which would make f2\(to) be equal to 0.7 today. Specifically, since p m redshifts while A does not, p m 
would then dominate over A at earlier times, with the universe then having to be in a decelerating 
phase rather than an accelerating one. For the given values of SIm^o) an d ^a(^o) such a switch between 
acceleration and deceleration would occur before z = 1, and so a monitoring of the Hubble plot above 
z = 1 would then test for such a switch over. As such, the very existence of a switch over near z = 1 
is a reflection of the cosmic coincidence, with initial conditions having to be fine-tuned in the early 
universe so that the switch from p m -dominated deceleration to A-dominated acceleration does not take 
place in the early universe itself but is delayed for no less than ten billion years. In Fig. 3 we have 
therefore plotted the higher z Hubble plot required of an Q M (t ) = 0.3, ^a(^o) — 0.7 universe, and if 
the standard model numbers are right, the luminosity function must be found to follow the indicated 
curve. 

To get a sense of the level of precision which is required to make this test, in Fig. 3 we have 
also plotted the expectations of the (fW(*o)> ^a(^o)) = (0,0) universe with q = and the pure A 
(^Af(^o)j ^a(^o)) = (0,0.37) universe with q Q = —0.37 which we discussed earlier, since the first of 
these two universes represents the coasting borderline case between acceleration and deceleration, while 
the latter represents the best fit permanently accelerating pure A universe. 49 Discriminating between 
the three cases (accelerating/ decelerating, coasting, and permanently accelerating) should eventually be 
achievable since at z = 2 these three typical universes respectively yield apparent magnitudes m = 26.75, 
27.04 and 27.17, while at z = 5 they yield m = 29.14, 30.25 and 30.40, and should thus enable us to 
determine what sort of phase the universe might be in at redshifts above one. 

Now while full discrimination between these various options would require obtaining as many data 
points above (ideally well above) z = 1 as there are below z — 1, it is of interest to note that a few z > 1 
data points have recently been obtained [JUIE] including one (SN 1997ff) as high as z = 1.75. Of these 
new data points nine have been identified as being particularly reliable and so we have extended our fits 
to include them. In order to perform the fits, we note that in j2] the 54 z < 1 data points are presented 
as apparent magnitudes m, while in jJT] the nine new z > 1 data points are presented as distance moduli 
m — M. With these two methods of characterizing the data being related via m = M + 25 + 51og 10 <iL 
where is the luminosity distance as measured in megaparsecs, in the fitting of the 54 z < 1 data 
points of [2. M is treated as a free parameter, with the notion that the supernovae are indeed standard 
candles then being confirmed by the fact that in the very good fits of Fig. 2 a common value for 
M actually is obtained. On now fitting the combined 63 (=54+9) data points in fits which allow M 

49 For the moment our interest in these particular cosmologies is only in having cosmologies which are not decelerating 
above z — 1 so as to provide a contrast with the standard cosmology in that region. However, shortly we shall see that 
such non-decelerating cosmologies are of interest in their own right, as they are relevant to the alternate conformal gravity 
solution to the cosmological constant problem which we present below. 
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to again vary for the 54 points in the sample, the standard (Qm^o), ^a(^o)) = (0.3,0.7) universe is 
found to yield x 2 — 74.5 and M = —19.15, with the best fi fc (t ) = universe fit being found to be 
the slightly lower mass (and thus more cosmically repulsive) (Qm^o)' ^a(^o)) = (0.21,0.79) universe 
with x 2 — 69.9 and M = —19.19. For comparison, the pure A (Qm^o), ^a(^o)) = (0,0.37) universe 
with g = —0.37 yields x 2 — 70.7 and M = —19.14, while the pure curvature (Qm^o), Q\(to)) = (0,0) 
universe with go — yields x 2 = 78.4 and M = —19.06. As we see then, while a standard Om(^o) > 0, 
^a(^o) = 1 _ ^a/(^o) cosmology with q(z > 1) > and an fiM^o) = 0, Rv(to) = 0.37 cosmology with 
q{z > 1) < both outperform the Qm^o) = 0, ^a(^o) = cosmology with q(z > 1) = in the z > 1 
region, statistically the q(z > 1) > and q(z > 1) < cosmologies are equivalent to each other, with 
the current z > 1 region data being just as well fitted by a model which is accelerating above z = 1 
as by one which is decelerating above it. Consequently, for the moment it is not possible to determine 
whether the universe is in an accelerating or a decelerating phase above z — 1, with both options being 
compatible with currently available data. To be able to display the nine z > 1 data points of jH] on 
the apparent magnitude plot given in Fig. 3, we have taken the fitted value of M (= —19.14) and 
used it to extract values for m for the nine points from their reported m — M distance moduli, with 
it being these inferred m values which are then displayed in the figure. While the x 2 fitting itself is 
not at all sensitive to the value used for H(t ) which is required for the luminosity distance d^, the 
extraction of an actual value for M itself from the fitting is. For instance with H(t ) = 72 km/sec/Mpc, 
in megaparsec units the quantity 5\og w (c/H(t )) evaluates to 18.099, while for H(t ) = 65 km/sec/Mpc 
it evaluates to 18.321, a quite substantial difference of 0.222. With an increasing value for the Hubble 
parameter leading to a decreasing 51og 10 (c/if (t )), an increasing value for the Hubble parameter thus 
also leads to an increasing derived M. In its turn an increased M leads to an increased (dimmer) value 
for m as extracted from a given distance modulus. On making this extraction, as we see from Fig. 3, 
the extracted z > 1 m values then precisely straddle the three indicated theoretical curves, to thereby 
prevent current data from being able to differentiate between acceleration and deceleration in the z > 1 
region. 50 With the SN 1997ff supernova at z — 1.75 being known to be lensed by two foreground 
galaxies which lie at z = 0.5 |42j . this particular supernova may actually be dimmer (larger m) than it 
appears. Filling in the z > 1 Hubble plot with a large enough sample of data points and allowing for 
systematic effects such as lensing should thus be very informative. 

Apart from the use of type la supernovae to explore the z > 1 region, a few other techniques have 
also been developed, with Daly and coworkers having used the very powerful extended FRII radio 
sources as standard yardsticks [J3J El and Schaefer having explored the use of gamma ray bursters 
as standard candles [IS]. For the FRII approach some 20 radio sources which go out to z = 2 have 
already been fully analyzed (11 being above z = 1), with there being some more which go out to 
z = 3 or so, while the gamma ray burster technique will be able to go out to z = 4.5. Global fits 
to the 54 z < 1 supernovae studied above and the 20 FRII radio galaxies have been made, with the 
(OmC^o); ^a(^o)) = (0.3,0.7) standard cosmology giving x 2 — 74.4 for the 74 total data points and an 
(f^M^o)) ^a(^o)) = (0,0.38) cosmology giving x 2 = 74.1 [46 . The available radio source data are thus 
seen to be just as compatible with both decelerating and accelerating universes in the z > 1 region 
as the supernovae data are. Further exploration of the z > 1 Hubble plot using all of the available 
techniques will thus be necessary in order to determine how the universe behaves in the z > 1 region 
and ascertain whether in fact ^a(^o) really is 0.7 in the real world. 

Having now described and diagnosed the dark matter and dark energy problems, we turn next to 
an analysis of some alternatives to the standard model in which these problems are resolved. And since 
the cosmological dark energy problem is an intrinsically general-relativistic one, to seek alternatives to 
dark matter and dark energy will oblige us to have to consider alternatives to Einstein gravity itself. 

50 The z > 1 data are however able to exclude the possibility that dust or cosmic evolution may be causing a dimming 
of the z < 1 supernovae, since both of those effects would dim the z > 1 supernovae even more, an effect which as noted 
in 02 the z > 1 data now exclude, with the behavior of the Hubble plot thus indeed being of cosmological origin. 
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8 Alternatives to Einstein gravity 



8.1 Pure metric theories with additional curvature dependent terms 

As we had noted earlier, while it is sufficient to use standard Einstein gravity to obtain the standard 
Schwarzschild metric phenomenology for the solar system, is not in fact necessary. Thus in seeking 
alternatives to Einstein gravity, as well as look for generalizations of it which contain the Einstein 
equations in an appropriate limit, one can also look for generalizations which are able to recover the 
Schwarzschild phenomenology while never reducing to the Einstein theory in any limit at all, with the 
more general rule thus being to generalize the Schwarzschild solution rather than the Einstein equations 
themselves. We shall thus look at both kinds of generalization (either reducing to the Einstein equations 
or bypassing them altogether) and after detailing possible options, shall then look for a fundamental 
principle which would enable us to unambiguously select just one option from amongst them. 

Within the framework of general covariant pure metric theories of gravity, the most straightforward 
generalization of the standard Einstein theory is to augment the Einstein-Hilbert action of Eq. (jlOHjl 
with additional general coordinate invariant pure metric terms, additional terms the smallness of whose 
coefficients or the specific structure of which (vanishing in Ricci flat geometries) would cause them to 
make negligible contributions in the solar system. The choice of possible additional terms is actually 
unlimited, since it not only encompasses appropriate contractions of arbitrarily high powers of the 
Riemann tensor (such as the fourth order actions of Eqs. (jl(J5j) and (jl(J6|) ). 51 provided candidate 
additions are able to avoid violating locality, in principle such additions need not even be expressible 
as a power series expansion of Riemann tensor terms at all. 52 Moreover, within the class of additions 
to the Einstein-Hilbert action, one must also include the addition that Einstein himself made, viz. a 
cosmological constant term, as well. Without a fundamental guiding principle almost nothing can be 
said either in favor or against any candidate such addition, with this freedom being due to the fact that 
it was not an appeal to a fundamental principle which led to the specific choice of the Einstein-Hilbert 
action in the first place. While it is sometimes stated that the choice of the Einstein-Hilbert action can 
be justified on grounds of simplicity, simplicity itself is not a law of nature, and the simplicity of the 
Einstein approach to gravity is not so much in the simplicity of its equations (or the ease with which 
they might be solved) but in the simplicity of its concepts. For Einstein, the simplicity of the concepts 
and the intrinsic elegance and beauty of the theory were paramount, and one could argue that in adding 
in dark matter and dark energy in the ad hoc way that is currently being done, one gives up much of 
the inner beauty of the theory. 



8.2 Additional fields 

Within the framework of general covariant theories of gravity, the next most straightforward general- 
ization of Einstein gravity is to introduce additional macroscopic gravitational fields to accompany the 
metric tensor itself, with the most common choice being additional scalar fields. In and of itself, the 
addition of covariantly coupled scalar fields does not imperil the validity of the equivalence principle, 
since instead of coupling through the It action of Eq. (jHJ), a test particle will now couple though the 
It action as augmented by the iy action of Eq. (fTTjl. to lead to a particle equation of motion 

{mc + ^ = (mc + kS) m + = -*S* U + (145) 



51 In a quantization of the standard Einstein theory an infinite set of such higher order terms is generated via radiative 
correction counter terms - though these terms then come with Planck scale coefficients which are indeed negligible 
compared to the Ti = Einstein term itself. 

52 The addition to the Einstein-Hilbert action of terms which behave as 1/R a a has been considered by 02j in an attempt 
to explain the acceleration of the universe as a gravitational effect rather than one which is generated by dark energy. 
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which still involves the covariant 4-acceleration D 2 x x /Dt 2 . While Eq. f!145|) and its Eq. (JT2j) antecedent 
can lead to departures from geodesic motion (an issue we return to below when we consider covariant 
generalizations of MOND), it still respects the equivalence principle requirement that the Christoffel 
symbol term and the ordinary acceleration term d 2 x x /dr 2 always appear in exactly the combination 
indicated. Within the class of scalar tensor theories the most venerable is the Brans-Dicke theory with 
scalar tensor sector action [48J 



where 5* is the scalar field and w is a constant. 53 In this theory there is no fundamental G, with G 
instead being induced as an effective gravitational constant in solutions to the theory in which 5* is non- 
zero. With S being able to take different values in different epochs, the Brans-Dicke theory therefore 
allows for a varying effective G, and with it then being this epoch-dependent effective G which will 
then normalize Qm in each epoch, while not solving the Friedmann equation fine-tuning problem the 
Brans-Dicke theory does at least allow one to reformulate it as a need to find a dynamics in which S 
could naturally change by 60 orders of magnitude from the early universe until today. While such a 
mechanism is a general one which actually holds for any theory in which an effective G is induced by 
a scalar field, it turns out that solar system constraints on the Brans Dicke theory itself are so tight, 
that there is no room to choose its parameter w in a way which could lead to any potential departure 
from the standard theory. While this is thus a shortcoming of the Brans-Dicke theory, it is not actually 
a generic shortcoming of all scalar tensor theories of gravity, only a shortcoming of the scalar tensor 
theory which is based on the particular action of Eq. (j!46|) . with every individual candidate scalar 
tensor theory of gravity needing to be evaluated on a case by case basis. 

8.3 Additional spacetime dimensions 

The third kind of generalization of the standard Einstein theory is one which has now become ex- 
tremely popular and widespread, namely increasing the number of spacetime dimensions. Such ideas 
originated with Kaluza and Klein in the nineteen twenties who considered the possible existence of a 
fifth dimension, with the concept of extra dimensions remerging much later following the development 
of superstring theories. Characteristic of these higher dimensional approaches is not just the existence 
of extra dimensions, but also of extra fields. The Kaluza-Klein approach is based on a 5-dimensional 
metric, a 15-component quantity whose 4-dimensional decomposition consists of a 10-component rank 
two tensor graviton, a 4-component vector field, and a 1-component scalar field, while the string theory 
approach typically involves scalar dilaton fields and not just rank two tensor gravitons. 54 To account 
for the apparent absence of any dimensions beyond four, it had generally been thought that all but four 
of the higher dimensions would be compactified down to Planck length size. More recently though it 
was recognized by Randall and Sundrum [HI E] that if we lived on a brane (viz. membrane) which was 
embedded in a higher dimensional anti-de Sitter bulk, the anti de Sitter bulk could (for at least some 
choices of brane matter fields) inhibit the propagation of gravitational signals into it no matter how 
large in size it might be. 

While these higher dimensional theories all contain a massless graviton, it does not follow that a 
4-dimensional observer confined to the brane will then see standard 4-dimensional Einstein gravity on 
it, with the gravity on the brane being modified by the very embedding itself, a modification which 
actually opens up a possible new way to address the dark matter and dark energy problems. In order to 
see how such modifications come about, we note first that before imposing any dynamics whatsoever, 

53 With its 1/5* dependence the scalar field kinetic energy is not expressible as a power series expansion in S. 
54 String theory typically also generates Planck scale higher order Riemann tensor action terms. 
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for the Gaussian normal metric 



ds 1 = g MN dx M dx N = dw 2 + q^dx^dx" (147) 



considered in the typical 5-dimensional bulk case, use only of the purely geometric Gauss embedding 
formula entails that the induced Einstein tensor ^G^ v on the brane (viz. the one calculated with the 
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induced metric g M „ seen by an observer at w = 0) is given by 

^Gfj, u = -G a pq a ^ v + -G^q^y — -G A A q flu — 

-K a a K, u + K a ^K av + l -{K a a fq^ - l -K afi K^ q^ (148) 

where Gmn is the 5-dimensional Einstein tensor (viz. the one calculated with the full 5-dimensional 
metric Qmn), E^ u is the particular projection 

E,u = q%q 5 u C%, 5 (149) 

of the 5-dimensional Weyl tensor Cabcd, and K^ u (the (fiu) component of the extrinsic curvature 
Kmn = (g A M ~ nAn M)(g B n ~ n B nN)nB;A °f the brane) is given by 

Kn* = Q a ^u n /3;a • (150) 

Similarly, in terms of the acceleration vector = n A riN-A which lies in the brane {n N = when 
n N n,N = 1), one can derive a second purely geometric relation 

n C K^ v . c -R BD q B ^q D u = -^ ) R^ v + ]^{a^ v + a u ^)-K A A K^-a^a v . (151) 

Consequently, no matter what form we may assume for the brane world energy-momentum tensor 
Tmn, we immediately see that even if we do take Gmn to obey the Einstein equations in the higher 
dimensional bulk space, viz. 

Gmn = —i^Tmn , (152) 

it does not at all follow that ^G^ u has to then obey the standard Einstein equations in the 4-space. 

To see what equation ^G^ does in fact obey, we note that if we take the bulk to possess a 5- 
dimensional cosmological constant A 5 and the brane to possess some energy-momentum tensor t^ v 
which is confined to it, viz. 

T M n = -M9mn + t^ v 8 m 8 v n 8{w) , (153) 

then since every term on the right-hand side of Eq. (j!51|) is continuous at the brane 57 it follows that 
the discontinuities on the left-hand side must cancel, 58 so that on crossing the brane we obtain 



dw d w K^ + k\ (t^ - £flW° tt ) 5 ( w ) 



r0+ 

, (154) 



55 Here we take w to denote the fifth coordinate, M and N to denote 5-dimensional coordinates and /i and v to denote 
the usual four coordinates on the brane, and have chosen the geometry so that the brane is located at w = and the 
normal to the brane is given by um = (0, 0, 0, 0, 1). In this coordinate system even though the components of the induced 
metric q^ v only range over the 4-space indices, <? M „ is nonetheless a function of all five coordinates. 

56 See e.g. where a more detailed treatment of the embedding issues and formulae discussed here may be found. 

57 Both ^ and a M lie in the plane of the brane, and even though the extrinsic curvature has a step function 
discontinuity at the brane, terms quadratic in it do not. 

58 Sincc Ka V is a step function, its covariant derivative is a delta function. 
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with the fully covariant Israel junction conditions 50 j 



K^iw = 0+) - K^(w = Q- 



I [11/ 



-q»v(w = 0)r a a 



(155) 



then emerging. 59 For the brane world set-up a Z 2 symmetry about the brane is imposed so that the 
extrinsic curvatures on the two sides of the brane are related by K^ v {w = + ) = —K^ v {w = 0"), with 
the extrinsic curvature at the brane then being given by 



K^{w = 0+) = -K^(w = 0- 



4 r i 



(156) 



where 



Given this discontinuity and the original 5-dimensional Einstein equations of Eq. (J152j) which fix the 
value of Gmn at the brane, the extraction of the continuous piece of Eq. ()148|) is then direct and yields 

(4) G F = ^4A 5qfiv {w = 0) - «fn^ - E^(w = 0) (157) 

where 

E^(w = 0) = \[E^{w = 0+) + E^{w = (T)] (158) 

and 

n M , = + ^r>^ + ig^(«, = 0)r aP r^ - ^q, u (w = 0)(r Q J 2 . (159) 

Finally, on introducing a brane cosmological constant A and decomposing the brane t^ u as 

T (a , = -\q t „,{w = ti) + Si W (160) 

where S^ v represents all other brane matter field sources, we find that Eq. (J159)) can then be rewritten 

as 

(4) C^ = A 4 g^(w = 0) - nlS^ - Kj-rr^ - E^(w = 0) (161) 

TV = -^,^/ + ^ Q a ^ + ^(^ = o)5 a ^-^( w = o)(^j 2 , 

4 = ^ , A 4 = ^(^A 2 + 6A 5 ) = ^(k 2 A + ^A 5 ) . (162) 

As constructed, Eq. (j!61j) thus represents the generalized Einstein equation seen by an observer on 
the brane, with the E^w = 0) and tt^ terms heralding a departure from a standard 4-dimensional 
gravity with an effective gravitational constant k\ and an effective cosmological constant A 4 . 60 With 
regard to the effective cosmological constant, we note that it is not A which is to serve as the cosmological 
constant for cosmology on the brane but rather A 4 , and since the anti-de Sitter bulk A 5 has to be 
negative, the effective A 4 has to be smaller than the A expected from ordinary 4-dimensional early 
universe phase transition physics. As we had noted earlier, a possible way to quench the contribution 
of the cosmological constant to cosmology is to have two cosmological constants, one being a dynamical 
one coming from the energy-momentum tensor (viz. A) and the second being an a priori fundamental 
one. As we now see, what we had thought of as a fundamental cosmological constant from the point 
of view of 4-dimensional physics might actually be a K5A5/2 term coming from the embedding in a 

59 In a D dimensional bulk and a D-space brane the factor 1/3 is replaced by 1/(D — 2). 

60 From the point of view of 4-dimensional physics, the massless 5-dimensional graviton of Eq. I|152f) decomposes into a 
massless 4-dimensional graviton and an infinite tower of massive 4-dimensional spin two modes jSJ |5] , with these massive 
modes then contributing to Eq. As such, these massive modes could also serve as potential dark matter candidates. 
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higher dimensional bulk. While embedding thus does provide a natural rationale for the existence of 
the K5A5/2 term, at the present time no rationale has yet been found for why the k|A 5 /2 term and 
the terms might be as close to each other as an £l\(to) = 0.7 universe would require. Within 

this extra dimension picture we note that it has also been suggested [52] that the acceleration of the 
universe might be explainable not as a cosmological constant effect at all but rather as a consequence 
of a leakage of gravitational signal into the extra dimensions. 

8.4 The modified Friedmann equations on the brane 

To get a sense of the relevance of the E^ v {w = 0) and 7i^ u terms in Eq. (|161|) . it is convenient to descend 
to the special case where S^ u is taken to be the perfect fluid 

Syj, = (p + p)U^U u + pq flv (w = 0) , (163) 

since in this case tt^ then also takes the form of the perfect fluid 

V, = (R + P)U^U U + Pq^{w = 0) (164) 

with energy density and pressure 

R= U ' P = ^(P 2 + 2 PP) ■ ( 165 ) 

Considering now the cosmological case where g M „(u» = 0) is taken to be a standard Robertson- Walker 
metric with scale factor a(t), the (00) component of Eq. (|161|) is then found to take the form 

- 3 -^^ = -A4-4p-^-E o(w = 0) , (166) 
a z 12 

with it thus being Eq. ()166|) which is to describe cosmic evolution on the brane, to be compared with 
the standard (unembedded) 4-dimensional Friedmann evolution equation 

3(<i 2 + k) , 9 . 
- 1 J- ' = -A 4 - 4p 167 

or 

where the k\p 2 /12 and Eqq(w = 0) terms are of course absent. Now while it is geometrically possible 
to embed a brane with a Robertson- Walker geometry into a 5-dimensional bulk with an anti-de Sitter 
geometry, it is not obligatory, with the bulk being able to have less than maximal 5-symmetry in general. 
Thus while the Eqq(w = 0) term will vanish if the bulk is anti-de Sitter (since the bulk Weyl tensor 
then vanishes), it should not in general be expected to do so. Its non- vanishing would thus entail an 
unavoidable modification of standard cosmology whose impact on the dark matter and dark energy 
problems could be of interest. Moreover, even if the bulk geometry is such that the Eqo(w = 0) term 
does in fact vanish, Eq. ()166|) will still differ from standard cosmology due to the presence of a term 
which is quadratic in the energy density p. Such terms will be of great importance when p is big, namely 
in the early universe, though in permanently expanding cosmologies such terms would redshift away 
much faster than the term which is linear in p and thus be nicely suppressed today. With the p 2 term 
leading to early universe modifications, it would be of interest to see how they might modify the initial 
condition fine-tuning flatness problem. 

While it might be thought that a late time embedded E 00 (w = 0) = cosmology would always have 
to behave the same way as an unembedded one, this is not necessarily the case. Thus in the specific 
case where the effective A 4 is negative, the Robertson- Walker 3-curvature k is zero and E 00 (w = 0) is 
zero, Eq. (|166|) is found to admit of the exact radiation era (p = A/a 4 ) solution 
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where b 2 = — K5A5/6, A 4 = —3H 2 , and k\ = K\{b 2 — H 2 ) 1 / 2 . As we see from Eq. (j!68j) . in this particular 
model a{t) is not in fact late time expanding. Rather, it only expands to a bounded maximum value 
given by 

<^t) = ^[b+(b 2 -H 2 )V 2 } . (169) 
With the ratio of the terms in Eq. (J166j) which are quadratic and linear in p being given by 

4p 2 H 2 

12k\ P ~ 2[b 2 -H 2 + b{b 2 - H 2 ) 1 ' 2 ] 



(170) 



at this maximum, we see that when H is of order b, the right-hand side of Eq. (jl70j) is of then of order 
one. Since the p 2 /p ~ 1/a 4 ratio is at its minimum when a(t) is at its maximum, we thus see that 
in this particular cosmology at no time is the p 2 term ever negligible compared to the p term, with 
this particular embedded cosmology never being able to approach the associated unembedded one at 
all. With Eqs. (|148|) and (|151|) being completely general purely geometric relations which will hold 
whenever there is a higher-dimensional embedding of any sort whatsoever (they do not even require the 
presence of a brane), despite the fact that higher dimensional models do possess a massless graviton, we 
see that such models could in principle yield modified 4-dimensional physics, and thus each candidate 
such higher-dimensional model has to be analyzed independently to see what kind of 4-dimensional 
physics it might lead to. 



8.5 Relativistic MOND theory 



Through the use of additional fields it is also possible to generalize Milgrom's non-relativistic MOND 
theory to the relativistic regime. The first step towards doing this was taken in |53J where a scalar field 
ij) was introduced with action 



8nGL 2 



dH-g) 1/2 f(L 2 9 a % a ^) 



(171) 



(/ is a scalar function and L is a constant), with this action then being added on to the standard 
Einstein-Hilbert one. Additionally, the scalar field coupling to a test particle was taken to be given by 



m f (—g a/3 dx a dx 13 ^ ^ , 



(172) 

with the test particle thus obeying Eq. ()12j) as written with S = e^. With the non-relativistic limit 
of / being taken to be the MOND function, the MOND formula of Eq. (|9fij) was then found to ensue 
with the MOND acceleration parameter being identified as ao = 1/L. While this treatment shows what 
is in principle needed in order to recover MOND (or in general to produce a covariant theory whose 
non-relativistic limit is not the standard second order Poisson equation but a modified version of it), as 
formulated the above specific treatment had two difficulties. One was that ip waves would propagate 
faster than light (the requisite / not being quadratic), and the other being that ip would have no 
effect on the propagation of light waves (the light cone is invariant under the conformal transformation 
9a/3 e 2 ^9ap) leaving the theory unable to explain gravitational lensing without dark matter. To 
resolve the acausality problem a second scalar field A was introduced 
(171) being replaced by 



153] with the action of Eq. 



iy>, a) = --i d*x(-gy/ 2 



if 



g a % a ^ + g aP A ]a A,^ + u(A 2 ) 



(173) 



where r\ was a real parameter and v a real valued function. Building on this idea Bekenstein [56] then 
found a a way to avoid the gravitational lensing problem as well in a generalization which involved 
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both scalar and vector fields. The theory which Bekenstein constructed contains the metric g^ v , a 
timelike 4-vector field U a which obeys the constraint g af3 U a U/3 = — 1, a propagating scalar field <fi and a 
non-propagating scalar field o. The action for the metric is taken to be the standard Einstein-Hilbert 
action, to which is added an action 



cr 



(174) 



for the scalars and an action 
K 



I{U a 



32nG 



dH-g) 1/2 



g^g^ {u a;fi - u Ka ) {u Pil 



u, 



2A 



(175) 



for the vector. (Here k, K, and i are constants, F is a general scalar function, and A is a Lagrange 
multiplier which enforces the constraint g a ^U a Uf3 = —1.) With this formulation, MOND can now 
be regarded as being a fully-fledged, fully consistent relativistic theory which retains all of MOND's 
non-relativistic features, and it will be of interest to see how it might address the dark energy and 
cosmological dark matter problems, problems which we noted earlier fall right in the MOND regime. 
While the above treatment does show how it is possible to embed the MOND formula of Eq. (j^Bj) in 
a fully relativistic setting, a drawback of the approach is that the function F (to which the MOND 
function fi(a/a ) of Eq. is related) is not specified by it, as the only constraint on it is that F{kGa 2 ) 
be a general coordinate scalar. It would thus be of interest to find a dynamics which could constrain 
the function F in a way that could naturally lead to the behavior which is phenomenologically required 
of the MOND function in its large and small argument limits. 



8.6 Modifications of the nature of the spacetime geometry itself 

Beyond the above more or less conventional generalizations of Einstein gravity (more Riemann tensor 
dependent terms, more fields, or more spacetime dimensions), 61 there are also some less orthodox ones 
which are also worthy of consideration, generalizations which modify the role played by the metric and 
the nature of the geometry itself. The first of these is to introduce torsion, with the Christoffel symbols 
then no longer being symmetric in their lower indices. Since such a modification would typically only 
be expected to be of relevance to microscopic physics, it is not thought to have any effect on the dark 
matter and dark energy problems. The second modification is to allow the metric to be non-symmetric 
in its indices, with then being a 16-component tensor with ten symmetric and six antisymmetric 
components. Such a proposal would not affect the line element ds 2 = g^dx^dx 11 at all since the 
symmetric dx^dx v can only couple to the symmetric part of g^. With the antisymmetric part of the 
metric having precisely the same number of degrees of freedom as the 6-component electromagnetic field 
strength (E and B) this generalization also opens up the possibility of unification of electromagnetism 
with gravity, a program which has been developed in detail by Kursunoglu jST] and generalized by 
him to encompass unification with the other fundamental forces as well jHH]- Models based on a non- 
symmetric metric have also been explored by Moffat, enabling him [521 EDI EHJ E2 EH] to provide a 
reasonable explanation of galactic rotation curve systematics and the accelerating universe without the 
need to invoke dark matter or dark energy. 

Perhaps the most radical departure from Einstein gravity was actually the first generalization of it 
ever considered, when only a couple of years after the very introduction of general relativity Weyl [U3] 
actually generalized Riemannian geometry itself. In what was the very first attempt at unification of 
electromagnetism with gravitation through metrication, Weyl introduced a local transformation which 

61 Not that any modification of Einstein gravity can really be described as being conventional. 
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he referred to as a "gauge" transformation under which the metric and the electromagnetic field would 
jointly transform as 

g^{x)^e 2a ^g, v {x) , (176) 

Afj,(x) -> A^(x) - e<9 M a(x) , (177) 

with gravitation and electromagnetism thus unifying by sharing a common a(x). Given such a joint 
transformation, Weyl then departed from the Riemannian geometry of general relativity and replaced 
it with a new geometry, "Weyl Geometry" , in which the Riemannian connection was generalized to the 
A M -dependent 

= \d Xa (d^gav + d v g atl - d a g^) + l -g x ° {g av A^ + g afl A v - g^A a ) , (178) 

so that instead of the covariant derivative of the metric being zero, it would instead be given by 

2A fJi 

ST-* = QusT + V + V = — • (179) 

e 

The utility of this generalized definition of covariant derivative is that Eq. ()179j) has the remarkable 
property of being invariant under the joint transformations of Eqs. ()176j) and (|177|) (transformations 
under which transforms into itself). With the covariant derivative of a tensor transforming as a 
vector, Weyl was thus able to connect the tensor g^ v and the vector A^ 1 in an intricate geometrical 
fashion, though at the price of departing from Riemannian geometry. 62 It was through the development 
of this theory that the notion of a gauge transformation was introduced into physics, with its usage 

62 In passing we note that there is no conflict with general coordinate invariance in having a connection which is not 
simply the standard pure metric based Christoffel symbol r^(C.S.) = (l/2)g Xa {d^g av + d v g alJL — d^g^). In fact, in 
general the connection can even be completely independent of the metric altogether and serve as a totally independent 
dynamical degree of freedom. All that is required of a general connection T x is that it transform the same way under 

a general coordinate transformation as the standard (C.S.), viz. as F x — ► (dx' x /dx p )(dx T /dx' fJ, )(dx a /dx' v )FP. a + 
(dx' x /dx p )(d 2 x p /dx' p dx'"), since then covariant derivatives constructed via its use will still transform as true general 
coordinate tensors. (Because the (l/e)g XrT (g^uA^ + g aih A v — g^vAa) term in Eq. I|178|) is a true tensor, then given the way 
r^„(C.S.) transforms, it follows that the quantity of Eq. i|178|) does indeed transform as a connection from which true 
tensor covariant derivatives such as g pv . v of Eq. I|179|) can be constructed.) However, in geometries with a metric and an 
independent connection, even though the metric will still give the invariant line element and be responsible for raising and 
lowering indices, true tensor quantities such as g^. T — d T g^ v J rF p ra g av +F" rr g <T ^ will not in general be zero, with it then not 
being possible to move the metric in and out of covariant derivatives in the usual Riemannian way. Additionally, in terms 
of a general connection, the Riemann tensor as constructed as iC„ re = (dT x ^/dx K ) — (dT^/dx") + F^F* — F^F^n will 
still transform as a true rank four tensor, with g pK R^\ K still being a general covariant scalar. Using this prescription for the 

Ricci scalar and forming covariant derivatives with the general F x v connection will thus still enable us to build an action 
which is a true general coordinate scalar even when the connection is not taken to be the Christoffel symbol or even not 
specified at all. And while the variation of such a scalar action with respect to the metric will automatically be a true rank 
two tensor, its independent (Palatini) variation with respect to an independent connection (in general a variation which is 
independent of that associated with the metric) will still be a true tensor too, since even though the connection is not itself 
a tensor, its variation is, as the (dx' x /dx p )(d 2 x p /dx' p dx' L ') term drops out in the transformation of the difference between 
any two connections which both transform as F x „ -> {dx' x / dx p ){dx T / dx'^dx" / dx ,l/ )T p Ta + {dx lX / dx p ){d 2 x p / dx'^dx'"). 
Consequently, dealing with a connection which is independent of the metric is just as covariant as taking the connection to 
be the standard Christoffel symbol, and doing so can thus generate departures from standard Riemannian-based gravity 
which are every bit as covariant. (In such cases it should be noted that if the matter field is taken to have an action which 
is independent of the connection altogether - such as the free Maxwell action - because of the matter field equations of 
motion, the covariant derivative of the matter field energy-momentum tensor with respect to the Christoffel-symbol based 
connection will vanish identically. However since this covariant derivative is linear in r^(C.S.), and since the covariant 
derivative of the matter field energy-momentum tensor with respect to F x is necessarily a true tensor, it follows that the 
covariant derivative of the matter field energy-momentum tensor with respect to r^„(C.S.) is a true tensor too, with its 
vanishing in one coordinate system entailing its vanishing in all other coordinate systems too.) 

While there is actually no difference between the Palatini approach and the standard approach when the action is just 
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in Eq. ()176|) entailing a change in the magnitude of g^ v and thus of its size (this being the meaning 
of the term "gauge"). With the Einstein- Hilbert action not being invariant under the transformation 
of Eq. (J176j) . Weyl's theory thus possessed neither the Einstein equations or Riemannian geometry, 
and was not so much a generalization of Einstein gravity as a rather substantial departure from it. 
Moreover, following the development of quantum mechanics with its complex-valued wave functions, 
Weyl was able to redefine his gauge transformation so that instead of being accompanied by Eq. ()176|) . 
the gauge transformation of Eq. (jl77j) would instead be accompanied by a change in the complex 
phase of the wave function (ip — > e ia ip). With Weyl's original invariance of Eq. (j!76)l ) entailing that 
all mass parameters would be zero identically (something which at times prior to the development of 
spontaneous symmetry breakdown was regarded as totally unacceptable), and with the use of Eq. (|177J) 
in accompaniment with ip —>■ e ia ip then becoming so fruitful, Weyl's geometric theory has essentially 
been set aside 63 - so much so in fact that the term "gauge invariance" is now taken to refer exclusively 
to complex phase transformations, with the transformations of Eq. (J176)) now being known as conformal 
or scale transformations instead. 

8.7 Conformal gravity 

In the course of developing his theory Weyl also discovered [HZ] a tensor with a remarkable geometric 
property, the so-called conformal or Weyl tensor 

Cx^ Vfi = Rx^un — 2 (gxuRfiK — gxnRfMu — 9^uR\k + g^Rxu) + g-^a (gxug^n — gxng^v) • (180) 

As constructed, this particular combination of the Riemann and Ricci tensors and the Ricci scalar has 
the property that under the local transformation of Eq. ()176|) it transforms as 

C\ VK {x) - C\ VK (x) (181) 

with all derivatives of a(x) dropping out identically. As such the Weyl tensor C x ^ lUK thus bears the 
same relation to conformal transformations as the Maxwell tensor does to gauge transformations, 
with the kinematic relation g^F^ = having as a counterpart the kinematic g /J,K C x R = 0, with the 
Weyl tensor being the traceless piece of the Riemann tensor. 

Given this particular property of the Weyl tensor, a far less radical version of Weyl's theory is 
suggested, one in which the invariance of Eq. (|176|) is retained but the unification with is dropped, 
so that the geometry is then a standard strictly Riemannian one in which the covariant derivative of the 
metric is zero as usual. We thus retain the metric as the gravitational field, have it couple covariantly 
in the usual way, but endow gravity with an additional symmetry beyond coordinate invariance, viz. 
conformal symmetry. As such, this theory is known as conformal gravity and has been pursued by 
various authors, with (as we shall see below) the interest of the present author being in applying it 
astrophysics and cosmology with a view to solving the dark matter and dark energy problems. The 
great appeal of this conformal symmetry is that its imposition actually leads to a unique choice of 

the Ricci scalar (since the stationary variation with respect to the connection then yields a condition whose solution 

requires the stationary T^ v to actually be r^ !/ (C.S.) after all), in cases where the action contains additional fields and 
cross-terms between these fields and the Ricci scalar, the Palatini approach can lead to equations of motion for the metric 
and the connection which are totally different from the ones in which the connection is taken to be given by the standard 
r^„(C.S.) from the outset. We thus note that taking advantage of this fact has enabled the authors of jSS] to develop a 
Palatini-based scalar-tensor generalization of the approach of |47| which also allows for cosmic acceleration to emerge as 
a purely gravitational effect. 

63 See however H^)] where it is suggested that Weyl's vector field might be become massive by a Higgs mechanism in 
which it combines with the component of the complex doublet scalar field whose vacuum expectation value spontaneously 
breaks the SU(2) x U(l) electroweak symmetry. 
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gravitational action, as there is one and only one action which is invariant under the local conformal 
transformation of Eq. (J176|) . namely the Weyl action 



Iw = -dig J d^x(-g) x/2 Cxn VK C XtJ,VK 

= -a g J d 4 x(- g y/ 2 \Rx, UK R x ^ K - ZR^RP* + ±(R 



a \2 



(182) 



where a g has to be dimensionless, with conformal gravity thus possessing a dimensionless coupling 
constant. 64 Moreover, not only does the conformal symmetry uniquely select out the Weyl action, it 
expressly forbids the presence of any fundamental cosmological term, and is thus a symmetry which 
is able to control the cosmological constant; 65 and as we shall see below, even after the conformal 
symmetry is spontaneously broken (as is needed to generate particle masses), the contribution of a then 
induced cosmological constant to cosmology will still be under control. 66 

However, in addition to forbidding a fundamental cosmological constant, the same conformal sym- 
metry also forbids a fundamental Planck mass, and thereby excludes any fundamental Newton constant 
or any fundamental Einstein-Hilbert action. Nonetheless, the conformal theory cannot be excluded on 
these grounds, since, as we had noted earlier, as long as a theory contains the Schwarzschild solution 
as an exterior solution, that will suffice to recover the standard solar system phenomenology. However, 
such an eventuality does not immediately appear likely for the action Iw since it involves the Riemann 
tensor, a tensor which does not vanish in the Ricci flat Schwarzschild solution. However, it turns out 
that the Lanczos Lagrangian 



,1/2 



Rx^R^ VK - w^et + (R a a ) 



happens to be a total divergence [o9], so that Iw can be rewritten as 

R^ K - \{R 



I w = -2a g J d 4 x(-gy/ 2 



a \2 



to now only involve the Ricci tensor and scalar. Variation of the action of Eq. 
the metric then yields 
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with respect to 
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where WPX and were respectively introduced in Eqs. (jl(J7j) and ()108|) . 67 Then with the variation 
of the matter field action Im with respect to the metric yielding a matter energy-momentum tensor 
T^, the variation of the total action Iw + Im with respect to the metric yields as equation of motion 



4:Q! r , 



;i86) 



The parallel between C and F^ K also carries over to the form of the action with the 



Jd*x(-g)^C\ VK C^« 

action being the conformal analog of the —(1/4) J d 4 x(—g) 1 / 2 F pK FP K action, an action whose coefficient is also dimen- 
sionless. 

65 It was the ability of conformal symmetry to control the cosmological constant which first attracted the present author 
to conformal gravity |68| . with its application to the dark matter and dark energy problems only coming later. 

66 Conformal invariance Ward identities are not modified by a change in vacuum, with the conformal energy- momentum 
tensor having to remain traceless under a change of vacuum, with the induced cosmological constant term then having 
to be neither bigger nor smaller than the other terms in T M „. 

67 Alternatively, one can first make the variation of the action Iw 3 — J d 4 x(— g) 1 ! 2 R\ fiVK R x,1UK to obtain ^3] 
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identical vanishing of variation Sl^/Sg^ of the Lanczos action II = j d A xL[ j yielding the kinematic identity 



2R^ aTp R»< JT P + {l/2)g^R aTpX R" T P x + AW 
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using Eq. I187|) one can then conclude that is not independent of and 



ARP a R va - 2gP v R ar R ar - 2R a a RP v + (1/2) (R a a ) 2 = 0, on 
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with use of the kinematic identity 



R^ v . a = ~{R a a ) wv + R^Rar - R l " a R vu (187) 



allowing one to rewrite Eq. ()186j) in the compact form 68 



C^Q\l\vK C^ VK R^ 



Finally, as we see from the form of and given in Eqs. (|107|) and (jlU8|) . the Schwarzschild 
solution is indeed an exact solution to the conformal theory in the T^ u = region, just as required. 

Now we had noted earlier that within the framework of gravitational theories which did not re- 
duce to the Einstein equations of motion in any limit whatsoever, as long as those theories admit the 
Schwarzschild solution they would suffice for solar system phenomenology. In principle such theories 
could include gravitational actions based on arbitrarily high powers of the Ricci tensor and Ricci scalar, 
with there thus being an infinite number of such options. With there also being an infinite number of 
options for theories which do contain the Einstein equations, the beauty of conformal symmetry is that 
it serves as a principle which chooses one and only one gravitational theory from amongst all possible 
theories, regardless of whether they contain the Einstein equations or not. Since the conformal symme- 
try sets to zero any fundamental cosmological constant, it thus addresses the most severe problem in 
physics head on, and is thus a reasonable starting point for developing a gravitational theory; and so in 
the following we shall look at its observational implications for astrophysics and cosmology. However, 
we note now that the ability of the conformal theory to address the cosmological constant problem so 
directly is that this symmetry simultaneously excludes a fundamental Einstein-Hilbert term as well; 
with the difficulties faced by the standard theory in addressing the cosmological constant problem per- 
haps being due to the fact that in any starting point which permits the Einstein-Hilbert action, there is 
no reason to exclude a cosmological constant term at all, with it being very hard to get rid of something 
which has no reason not to be there. 

In the conformal theory Eq. (|186|) is to replace the standard Einstein equations of Eq. and with 
the rank two gravitational tensor W^ u being kinematically traceless, covariantly conserved, and having 
to transform as 

W^{x) -> e-^W^ix) (189) 

under the conformal transformation of Eq. (|176j) . it follows that the consistency of the conformal theory 
requires its energy-momentum tensor to be traceless, covariantly conserved and transform as 

T^{x) -> e- 6a{x) T^{x) (190) 

too. The consistency of the theory thus sharply constrains the energy-momentum tensors that are 
allowable, with there being no freedom to introduce new ad hoc sources (such as dark matter or 
dark energy) at will, with a typical allowed energy-momentum tensor being the generic one we gave 
earlier as Eq. 1)61)). Moreover, our use of a conformal invariant action for gravity dovetails with 
our earlier discussion of the constraints of an SU(3) x SU(2) x £7(1) invariant theory of massless 
fermions and gauge bosons on the structure of the energy-momentum tensor. Specifically, even though 
the fermion spin connection was only introduced into physics in order to make the fermion kinetic 
energy term — / d i x(—g) 1 ^ 2 i'ijj'y^(x)[d fl + r fl (x)]ij} be general coordinate invariant, nonetheless under 
ip(x) — > e~ 3a ( x ^ 2 ijj(x), V£(x) —> e~ a ^V^(x) (V a M (x) is the vierbein from which ^^(x) = V^{x)^f a 
and Tpix) = [ f y v (x),d fi j v (x)]/8 — [j v (x),'y a (x)]r^ a/ /8 are constructed) this term happens to be locally 

68 Even though C\^ VK {x) itself acquires no derivatives of ce(x) under a conformal transformation, because of the way the 
Christoffel symbols transform its covariant derivatives do, with it being the so-called Bach tensor combination 2C M ' W . A . f ,— 
C^ VK Rx K in which all derivatives of a(x) drop out. 
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conformal invariant as well. Similarly, the general coordinate invariant Maxwell action is conformal 
invariant too. With the full symmetry of the lightcone being the conformal group 0(4, 2), and with the 
conformal group being homomorphic to the Dirac algebra SU(2, 2), we see that requiring massless gauge 
bosons and massless fermions to couple covariantly to geometry forces them to couple to it conformally 
as well. The only exception to this requirement is the scalar field kinetic energy, since the minimally 
coupled scalar field kinetic energy is not invariant under S(x) — > e - a i x ) S(x), g^x) — > e 2a( - x ^g^ v (x). 



Rather it is only the combination — / d^x^—g) 1 ^ 2 
conformal theory which will prove central in the fol 



S^S.J2 - S 2 i^/12j which is, an aspect of the 
owing. With conformal gravity we thus resolve the 



conflict between standard gravity and the tracelessness of T^ v which is required of a theory built out of 
massless fermions and gauge bosons, with a traceless T^ v being coupled in Eq. (|186|) to a gravitational 
rank two tensor W^ v which is itself traceless too. Interestingly, while conformal symmetry actually 
forces the SU(3) x SU(2) x U(l) theory to be second order in the first place, 69 this same symmetry 
forces gravity itself to be fourth order. Thus while conformal invariance requires the gravitational side 
of Eq. (|186|) be fourth order, the same symmetry obliges its matter side to be second order, so that 
with conformal invariance the matter T^ v is of the same order (viz. second) as it is in the standard 
gravitational theory. Having thus singled out and motivated the use of conformal gravity, we turn now 
to see how it fares in addressing the dark matter and dark energy problems, an exercise which is of 
interest not just in and of itself, but also because it will uncover some generic features which could serve 
as a guide for other theoretical approaches which attempt to address these same problems. 



9 Alternatives to dark matter 

9.1 Modifying Newtonian gravity at large distances or small accelerations 

As it started to become clear that there were mass discrepancies in galaxies, a few physicists started 
to entertain the possibility that Newton's law of gravity might need modifying on large distance scales, 
with some early ideas along these lines being made by Finzi (701, Tohline [71], Kuhn and Kruglyak 
|72j and Sanders [75] . with some related work being made by Fahr (Tlj and Gessner [73]. Typical of 
these approaches was the introduction of an extra potential which was to accompany the Newtonian 
one, and while such approaches could in principle even be covariantized via the use of Eq. (j!45|) and 
thereby be made compatible with the relativity principle, by and large approaches such as these in the 
end tended not to fit data well enough to make them compelling. 70 More compelling was Milgrom's 
MOND approach because, without any need for dark matter, it could encompass a large amount of 
data extremely well via the use of a very compact formula which involved a very small number of 
assumptions and a universal acceleration scale criterion which not only explained why departures from 
the luminous Newtonian expectation are to occur for galaxies, but also why none are to be expected for 
the solar system. MOND fits to the data sample of Fig. 1 for instance are spot on [21], as are its fits 
to a large sample of low surface brightness galaxies |77j, with reasonable fitting being found for a quite 
wide range and variety of non-relativistic systems from molecular clouds to clusters of galaxies |78j . 
systems which in first approximation are found to obey the fundamental v ~ M 1//4 relation between 
mass and velocity which the MOND theory possesses. 71 And with MOND having now been successfully 
covariantized, it should be taken seriously since, as noted, it has captured a basic truth of nature, 

An action such as n f ^xi-gY^lF^F^} 2 would still be gauge invariant but would not be conformal invariant, and 
would have to have a dimcnsionful coupling constant k. The absence of terms such as these is simply assumed by fiat in 
the standard SU (3) x SU(2) x U(l) theory, with it being conformal invariance which actually provides a rationale for their 
absence, to thereby secure the renormalizability of the standard model of strong, electromagnetic and weak interactions. 

70 For a more detailed discussion of candidate alternate theories of this type and of the challenges they face see e.g. [76| . 

71 In passing we note that a similar mass- velocity relation was also found to hold |79| in the alternate conformal gravity 
theory which is described below. 
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namely that it is a universal acceleration scale which determines when departures from the luminous 
Newtonian expectation are to occur. 



9.2 Modifying Einstein gravity — the conformal gravity approach 

While one can develop non-relativistic theories by starting with the data themselves, in doing so it 
is not initially clear what covariantized relativistic theory might ultimately ensue or whether such an 
approach would ever be able to address problems such as the cosmological constant problem, a problem 
which is intrinsically relativistic. On the other hand if one starts with a relativistic theory which is 
chosen because it does address the cosmological constant problem and then works down, one does not 
know what non-relativistic limit one might then encounter. Encouraged by the fact that conformal 
gravity could address the cosmological constant problem, Mannheim and Kazanas [SOI EI] set out to 
determine the non-relativistic limit of the conformal theory, and did so at that time with a quite limited 
objective, namely to see whether a theory which did not contain the Einstein equations could nonetheless 
still recover the standard Newton/Schwarzschild phenomenology. 

To this end they endeavored to find an exact conformal gravity analog of the Schwarzschild exterior 
and interior solutions to standard gravity, viz. they tried to solve the equation W^ u = T^ v j^a g for 
a static, spherically symmetric source. While this meant having to handle a fourth order differential 
equation, in turned out that it was possible to greatly simplify the problem through use of the under- 
lying conformal symmetry that the theory possessed. Specifically, they noted that under the general 
coordinate transformation 



p = p(r) , B{r) 
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p 2 (r) 



with initially arbitrary function p(r), the general static, spherically symmetric line element 

ds 2 = -b(p)dt 2 + a(p)dp 2 + p 2 dVt 2 



would be brought to the form 

ds 2 - - 

On now choosing p{r) according to 
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the function A(r) would then be set equal to 1/B(r), with the line element then being brought to the 
convenient form 

p2{r) -B(rW + 



ds' 



< 2 dr2 + r 2 dQ 2 



Br) 



(195) 



While Eq. ()195|) is coordinate equivalent to Eq. ()192|) with the two general functions a(p) and b(p) 
having been traded for two equally general functions p{r) and B(r), its utility lies in the fact that one 
of these two functions appears purely as an overall multiplier in Eq. ()195j) . Consequently, through use 
of Eqs. (j!89J) and ()190|) the function p(r) can be removed from the theory altogether (i.e. gauged away) 
via a conformal transformation, with the full kinematic content of the conformal theory then being 
contained in the line element 

dr 2 

-B(r)dt 2 + ^— + r 2 dQ 2 (196) 
B(r) 



ds 2 



2 Such approaches while phenomenological should not be thought of as being any more ad hoc than dark matter itself. 
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in the static, spherically symmetric case. Comparing with the standard theory, we see that instead of 
having deal with coupled second order equations for two independent functions, we instead have to deal 
with a fourth order equation for one alone, an equitable enough trade-off. 

Evaluating the form that W^ u takes in the line element of Eq. (|196|) is straightforward though 
lengthy, and leads to 
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for its components of interest. Combining Eqs. (|197|) and ()198|) then yields the remarkably simple 

O A TDlll 1 

i- (W° - W r r ) = B"" + — = -(rB)"" = W A B , (199) 
B v ' r r 

so that in terms of the convenient source function f(r) defined via 

3 



/w = i^w( r °»-^) (200) 

the equation of motion of Eq. ()186j) can then be written in the extremely compact form 

V 4 S(r) = f{r) . (201) 

The remarkably simple equation of motion of Eq. ()201j) contains the full dynamical content of the 
conformal theory in the static, spherically symmetric case, an equation of motion which is completely 
exact and which has not required the use of any perturbative approximations whatsoever. While we 
have not advocated a preference for one theory over another on the grounds of calculational simplicity, 
preferring to make assessments on the basis of conceptual simplicity instead, nonetheless from a purely 
calculational standpoint, Eq. (j201|) can claim to be just as simple as (and perhaps even simpler than) 
its standard gravity counterpart given in Eqs. ()44j) and ([45)1 . The solution to Eq. ()201|) is readily given 
as the analog of the solution to Eq. (Jgfij) given earlier, and can be written as 

B(r > R) = -- [ R dr'r ,2 f(r') - — (* dr'r' A f(r') + w - kr 2 , 
2 jo 6r Jo 

B(r < R) = -- I" dr'r' 2 f(r')- — I" dr'r' 4 f(r') 
2 Jo 6r Jo 

_I f dr'r' 3 f(r')-— [ dr'r' f{r') + w - kr 2 , (202) 

2 Jr 6 Jr 

where the w — kr 2 term is the general solution to the homogeneous V 4 B(r) = equation. Finally, on 
defining 

7 = -^ f R dr'r' 2 f(r') , 2(3 = \ f* dr'r' 4 f(r') , (203) 



2 Jo J v ' ' ' 6 Jo 
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on dropping the kr 2 term (as it does not couple to the source) and setting w — 1, we see that up to 
conformal equivalence the conformal gravity metric exterior to a static, spherically source is thus given 
without any approximation at all as 

1 23 

B {r> R) = -g m = — = l-^ + 1 r . (204) 

Qrr f 

In the region where 2(3 /r ^> we thus nicely recover the Schwarzschild solution, while seeing that 
departures from it occur only at large distances and not at small ones, with the standard solar system 
Schwarzschild phenomenology thus being preserved. 73 



9.3 Impact of the global Hubble flow on galactic rotation curves 

With both (3 and 7 being related to moments of the energy-momentum tensor of the source, and with 
neither moment having any reason to vanish, a general conformal gravity source would then furnish 
both terms, with a source such as a star then producing a non-relativistic potential of the form 

nO = ~ + 2 y r > (205) 

to nicely yield none other than the potential of Eq. (|95|) whose merits were discussed earlier. In an 
application of the potential V*(r) to spiral galaxies, we need to integrate the V*(r) potential over a 
disk of iV* stars distributed as S(-R) = £ e~ R ' R °, with use of Eq. (JA15j) then enabling us to generalize 
Eq. (fTT5j) to 
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hM^K W« ) . (206) 



2Rq \2RqJ \2R 
With the large R limit of v 2 um /R then being given by 

R R 2 + 2 ' 1 U ° 

we see that we precisely recover the two iV*-dependent terms present in the formula of Eq. (|125|) that 
was extracted purely phenomenologically from a study of the centripetal accelerations of the last data 
points in each of the galaxies in the eleven galaxy sample of Fig. 1. The potential of the conformal 
theory thus automatically captures much of the essence of Eq. (|125j) . 

However, absent from Eq. (207) is the iV*-independent 7oc 2 /2 term which we saw was crucial 
to the validity of Eq. (j!25j) . and finding its origin is quite subtle and not at all easy. Specifically, 
intrinsic to the use of Newtonian gravity and to the intuition one then acquires through its repeated 
use is a theorem due to Newton, namely that if one has a spherical distribution of matter, the net 



73 Whilc the vanishing of the Ricci tensor entails the vanishing of its derivatives, the vanishing of the derivatives can 
be obtained without the vanishing of the Ricci tensor itself, with solutions other than the Schwarzschild solution then 
being possible. The constructive approach which leads to Eq. 1204(1 shows that not only is the metric of Eq. I|2()4|l with 
its the linear potential an exact exterior solution to the conformal theory, such a metric categorizes all departures from 
Schwarzschild, as every solution to the conformal theory is conformally equivalent to it. In solving the exterior fourth 
order Laplace equation V 4 £>(r) = one would obtain EDI the solution of Eq. I|2U4|) as expressed in terms of what 
would then be two free integration constants (3 and 7, with it being the matching 81 to the interior region via the fourth 
order Poisson equation V 4 i?(r) = f(r) which allows one to express these integration constants in terms of moments of 
the source and establish their physical meaning. 



62 



Newtonian potential produced at a given point depends only on the matter interior to its location, with 
the potentials due to all the points exterior to it mutually cancelling. In fact so central is this notion to 
our thinking, that when a mass discrepancy is found we immediately introduce dark matter in precisely 
the region where the discrepancy is found, a strictly local approach to gravity. However, this theorem 
of Newton only holds for a 1/r potential, and does not hold in any other case. Consequently, in any 
modification to Newton whatsoever, one cannot exclude the matter exterior to the point of interest, 
viz. one cannot exclude the effect of the potentials due to the rest of matter in the universe. Moreover, 
in modifications to Newton which grow with distance, these effects will be quite pronounced, with the 
biggest contributors to the gravitational potential at any given point then being not the matter sources 
which are the nearest but rather those which are most distant, a thus global gravitational effect, one 
which is quite reminiscent of Mach's principle. The contributions of such distant objects would then 
not depend on the mass of the local galaxy of interest, and since those distant objects form the Hubble 
flow, one would even expect their net contribution to come with a cosmologically relevant scale. From 
the rest of the universe then we can thus naturally expect to obtain some form of iV*-independent 
7oc 2 /2 type term with strength of order the phenomenologically obtained value 70 = 3.06 x 10~ 30 cm _1 
cited earlier. 74 However, what does not follow from this discussion is why the net effect of the entire 
rest of the universe would be to explicitly yield a linear potential term of the form V = 7oC 2 i?/2, as 
it is precisely such a contribution and no other which leads to a centripetal acceleration of the form 
v 2 /R = 7 c 2 /2. 75 

To explicitly obtain this global 7oc 2 i?/2 potential, we need to ask how the comoving Hubble flow 
would look to an observer who uses a coordinate system in which the center of any chosen galaxy is taken 
to be at rest. With a Robert son- Walker geometry being homogeneous and isotropic, every point in the 
geometry can serve as the origin of coordinates, and so we can take the center of any given comoving 
galaxy of interest to serve as that center, with it being the presence of the galaxy of interest itself which 
provides this choice for the origin of coordinates. Having fixed the origin of coordinates this way we now 
need to rewrite the comoving Robertson- Walker metric in static Schwarzschild geometry coordinates as 
referred to this same origin. 76 To this end we note that the general coordinate transformation 

P * f dT (208) 
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(209) 



to yield a metric which we recognize as being conformal to a Robertson- Walker metric with spatial 
3-curvature k = — 7 2 /4 (the spatial part of the Robertson- Walker metric being written in isotropic co- 
ordinates here). Now since a Robertson- Walker metric happens to be conformal to flat, the Weyl tensor 

74 For a galaxy containing N* stars, the distance scale at which the galactic potential which they produce becomes 
of order one is given by r ~ 1/JV*7*, a scale which for a typical N* — 10 11 galaxy is numerically precisely of order 
10 _30 cm _1 , with a cosmological scale of order I/70 thus being precisely the scale on which galactic potentials become 
strong. 

75 An application of the purely local galactic Eq. (|206|) alone to the eleven galaxy sample of Fig. (1) was actually found 
|8HI I84| to fit the shapes of the rotation curves extremely well, but with the data wanting the strength of the total galactic 
j*N* to be universal rather than 7* itself. It took the present author quite some time to grasp that the apparent need 
for a universal 7*jV* in the fits was actually due to the fact that in fitting data, Eq. (|206|l was trying to simulate the 
contribution of a universal, -ZV*-independent, 7 contribution coming from the rest of the universe. 

76 The Robertson- Walker geometry is spherically symmetric about every point while the Schwarzschild geometry is only 
spherically symmetric about a single point, with the Robertson- Walker geometry being spherically symmetric about that 
one particular point too. 
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associated with it would not only be zero, but would even remain so under conformal transformations 
such as ^ 

9tM ^ J R 2 (r)(l- 7oP /4) 2 ^ ' (210) 

transformations which we are free to make in the conformal theory; with a comoving Robertson- Walker 
geometry with k = — 7 2 /4 thus being coordinate and conformal equivalent to a static Schwarzschild 
coordinate system metric with none other than a linear potential. As seen by an observer at rest then 
it follows that the comoving Hubble flow acts precisely as a universal linear potential. Moreover, 
not only do we see that this effect is expressly associated with a negative spatial 3-curvature, 77 we shall 
show below that k actually is negative (and necessarily so in fact) in the cosmology associated with the 
same conformal theory. We thus see that despite our familiarity with the local nature of Newtonian 
gravity, global cosmology can have a local observable effect in galaxies, an effect which is of a quite 
general nature, and which should thus be present in every theory of gravity, and not just the conformal 
theory being considered here. 

With the linear potential of Eq. (J2U9)) being centered at the center of any given spiral galaxy, 78 in 
the weak gravity limit we can directly add it on to the luminous contribution of Eq. ()206|) . with the 
total centripetal acceleration seen by a particle in the galaxy being given by 

Vtot _ tot _ lum , lum , 7oC 2 _ N*j3*C 2 R 

-£--0 -9p +9j +-3- - 2R s 



IRq) \2RqJ \2RqJ \2Rq 



^/,f J- WJ-U^ , (21!) 



2Rq \2RqJ \2R 

with the g l a m and g l ^ m terms being due to the luminous material within the galaxy and the 7oC 2 /2 term 
being due to the luminous material in the rest of the universe. With the large R limit of v^ ot /R being 
given by 

vj t N*f3*c 2 N*j*c 2 , 7oC 2 . . 

R R* + 2 + 2 ' 1 ] 

we see that from the above g tot we precisely recover the phenomenologically established Eq. (j!25j) . with 
conformal gravity thus providing a theoretical rationale for its validity. With the 7* and 70 parameters 
already being fixed by the application of Eq. (jl25j) to the furthest data point in each galaxy, and with 
the mass to light ratios of the luminous matter distributions in each galaxy already being fixed by 
the maximum luminous disk prescription for the inner regions of the rotation curves, there is no more 
freedom left in the theory. The application of Eq. (|211jl to the entire 280 or so data points in the eleven 
galaxy sample of Fig. 1 is thus fully prescribed, with it leading to the fits of |2SJ which are displayed 
in Fig. 1. As we see, the conformal theory fully captures the systematics of the data, and just like 
the MOND fits, shows that it is possible to fit galactic rotation curve data without needing to use any 
galactic dark matter at all. 79 



9.4 Comparison of the conformal gravity and MOND fits 

To understand why fits to those of the rotation curves which are flat (viz. the bright spirals) can be 
obtained in a theory which is based on rising potentials, we recall that the bright spirals are all Freeman 

77 If k is zero there is no effect, while if k is positive, the transformation of Eq. H2()8(l would lead to a pure imaginary 
value for 70. 

In a coordinate system in which any given comoving galaxy is at rest, no other comoving galaxy could simultaneously 
be at rest also. Nonetheless, since every comoving galaxy can serve as the origin of the Robertson- Walker coordinate 
system, for the purposes of determining the motion of particles with respect to the center of any particular galaxy, we 
can refer Eq. 1|2U9[) with its universal 70 to the center of the galaxy of interest each and every time. 

79 As with the analogous MOND case, the numerical values obtained for the conformal gravity parameters used in Eq. 
(|211[1 ensure that no modifications to the luminous Newtonian expectation are to occur on solar system distance scales. 
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limit galaxies with a common central surface brightness £q\ and thus a common value for N*/2tcRI. 
In a theory with rising potentials, such potentials make a quite small contribution in the inner regions 
of the rotation curves, with the peak at R = 2.2R Q being controlled almost completely by the luminous 
Newtonian contribution, to thus yield 

% ~ = 0.8n^(3*R (213) 

c 2 Ro 

at the peak. With reference to say the typical galaxy NGC 3198, we match the inner region peak value 
at R = 2.2Rq to the value at the furthest data point at R = 10R where Eq. ()125|) holds, to obtain 



4N*3* c 8*N* 

^— = O.SvrS^^o = + 5R ( 7 *iV* + 7o ) , (214) 



and thus 



fl 06 N* R* 

7 *iV* + 7o = U ' U0 P =0.12^(3* . (215) 
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With 7*iV* and 70 being of the same order for such galaxies, the universality of 70 and are thus 
correlated, to not only automatically enable all Freeman limit galaxies to match their inner and outer 
regions (something which we recall is completely contrived in dark matter fits), but also to suggest 
that the very existence of the universal Eq 7 itself is of cosmological origin, and that it might therefore 
emerge naturally in a cosmologically based theory of galaxy formation. Additionally, we also note from 
Fig. 1 that in the intermediate region of the NGC 3198 rotation curve at around R = 6R0 or so, 
the luminous Newtonian contribution is at about half of its inner region peak value, while the linear 
potential contribution is at about half of its outer region value, with the sum of their contributions 
thus being equal to the value of Vj 2 um obtained in the inner and outer regions themselves. (With one 
contribution falling and the other rising, there has to be some intermediate point where they cross.) 
Matching the inner and outer regions to each other via Eq. (J214)) then automatically takes care of 
matching the outer and inner regions to the intermediate region as well, to thus make the rotation 
curve pretty flat over the entire 2.2i? < R < 10 Rq region (again we recall how completely contrived 
this particular matching is in dark matter fits). While such flatness is thus naturally achieved in the 
conformal theory, it is important to stress that it is only achievable out to some maximum value of R 
before the rise due to the linear potential wins out over the falling Newtonian contribution, an aspect 
of the theory which is immediately evident in the fits to the lower luminosity galaxies shown in Fig. 1. 
The view of the conformal theory then is that the currently rising low luminosity rotation curves will 
continue to rise while the currently flat high luminosity rotation curves will eventually start to rise at 
large enough R. Such an expectation stands in sharp contrast to both the dark matter and MOND 
expectations where the rising low luminosity rotation curves are expected to flatten off and the flat 
high luminosity rotation curves are expected to stay flat. Extending galactic rotation curves to larger 
distances from the centers of galaxies could therefore be quite instructive. 

The essence of the conformal gravity fits presented above lies in the specific role played by cosmology, 
with the view of conformal gravity being that as a test particle orbits a galaxy, it not only sees the 
local field produced by the luminous material within that selfsame galaxy, the test particle also sees 
the gravitational field produced by the rest of the universe. Particles orbiting a galaxy thus probe the 
global cosmology and thereby enable us to measure the global spatial curvature of the universe, which 
from the measured value for 70 is thus given by k = —2.3 x 10 -60 cm -2 . In the conformal theory the rest 
of the universe thus replaces galactic dark matter, with the view of the conformal theory being that the 
presence of dark matter in galaxies is nothing more than an artifact which is engendered by trying to 
describe a global phenomenon in purely local terms. 80 
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This same could of course equally be said of dark matter in even larger systems such as clusters of galaxies where 
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With Eq. ()211j) being writable as 




(216) 



we additionally see that we can write it in the form of Eq. (jHSj), viz. 
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u{x) = 1 + 
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We thus recognize Eq. (j216J) as being in the form of a MOND type equation with 7oc 2 /2 playing the role 
of a universal acceleration, but with the local g^ m being replaced by the entire local g lum = g 1 ^ + g l ^ m 
which is due to both the Newtonian and the linear potentials of the luminous matter within the galaxy. 
Conformal gravity thus provides a rationale for why there should be a universal acceleration in the first 
place, and also for why it should be associated with a cosmic scale. 81 Moreover, while it had been noted 
earlier that the universal acceleration of MOND had a magnitude of order cHq, this was quite puzzling 
since H is epoch-dependent. However, cosmology can actually supply not just one but two natural 
scales, viz. both H Q and k, and unlike H , the spatial 3-curvature A; is a true general coordinate scalar 
which is not epoch-dependent at all. It is thus k which is the natural quantity with which to associate 
a universal scale, with it being k which produces a global imprint on galactic rotation curves. The 
success of the conformal gravity fits thus suggests a role for cosmology in the interpretation of galactic 
rotation curves as well as a role for the spatial curvature of the universe. We believe these to be general 
features which should be sought in any theory which attempts to explain galactic rotation curves, with 
the absence of a fundamental scale such as k being a possible shortcoming of the standard hierarchical 
= 0, ACDM models of dark matter halos. 82 

As a final comment on the conformal gravity fits to galactic rotation curves which we have presented 
here, we note that even more remarkable than the fits themselves is that the conformal theory was never 
constructed for this particular purpose. Rather, the theory was selected for an entirely different purpose 
(controlling the cosmological constant), and when Mannheim and Kazanas set out to determine the non- 
relativistic limit of the conformal theory, they had no idea what they would find (as indicated earlier, 
they had only wanted to see whether the theory could recover the standard Newton/Schwarzschild 
phenomenology), and had no inkling at all as to how the theory might then impact on galactic rotation 
curves. The fact that the theory then does impact so well on the rotation curves suggests that one 
should therefore give the theory serious consideration. With the motivation for the conformal theory 
being the cosmological constant problem, we turn now to a discussion of how the theory treats that 
particular problem, and in the fitting to the Hubble plot which we shall describe below, we shall again 

the natural interplay in the conformal theory between local and global effects would be even more pronounced. With a 
cluster of galaxies being a large local inhomogeneity in a homogeneous global cosmological background, a full treatment 
of clusters in the conformal theory has to await the development of a treatment of the growth of fluctuations in conformal 
cosmology. 

81 While the MOND and conformal gravity theories agree on the need for a universal scale in galaxies, they depart on 
how the theories are to behave in the region where the existence of the scale is of consequence. (In passing we note that 
the numerical relation between the oo and 70c 2 parameters may be obtained by looking at the behavior of the typical 
galaxy NGC 3198 in the intermediate region near R = 6Rq where the Newtonian and linear contributions cross each other 
in Fig. 1. With 7*iV* being of order 70 for this galaxy, at the crossing point we have (c.f. Eq. H212|) ) P*N*/R 2 ~ 70 
and thus v 2 /R ~ 2 7o c 2 . Comparing this with the MOND value (c.f. Eq. (JTOH) ofv 2 /R ~ (a ) 1/2 (/3* 7V*c 2 /,R 2 ) 1/2 then 
yields the relation clq — 4 7 oc 2 , a relation which numerically is closely obeyed.) 

82 Whatever one's views on alternate theories such as MOND or conformal gravity, if the answer is to be dark matter, 
then dark matter theory really has to be able to reproduce the regularities in the galactic rotation curve data which these 
alternate theories so readily capture, and has to do so with the same easy facility. 
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find the theory working remarkably well on data it had not at all been designed for, data which actually 
did not even exist at the time [UHl EEj conformal cosmology was first developed. 

10 Alternatives to dark energy 

10.1 The conformal gravity alternative to dark energy 

In order to identify possible alternatives to dark energy, we first present the conformal gravity theory 
treatment of the issue, and then extract some generic features from it. In applying conformal gravity to 
cosmology we note that since the Weyl tensor vanishes in a Robertson- Walker geometry, in conformal 
cosmology the entire left-hand side of Eq. (jl88j) vanishes identically, with the equation of motion for 
conformal cosmology thus reducing to the extremely simple 

= . (219) 

With this vanishing of the cosmological T^, we see immediately that in the conformal theory the zero 
of energy is completely determined, to thus give us control of the cosmological constant. 83 To appreciate 
the implications of Eq. (|219|) , we consider as cosmological energy- momentum tensor the typical matter 
field conformal energy-momentum tensor given in Eq. (|64j) . With this energy- momentum tensor having 
to transform according to Eq. (jl90|) under a conformal transformation, its vanishing will persist after 
any conformal transformation is made. With the scalar field transforming as S(x) — > e~ a ^S(x) under 
a conformal transformation, in field configurations in which S(x) is non-zero, the full content of the 
theory can be obtained by working in the particular gauge in which the scalar field takes the constant 
value So. In such a gauge, use of the matter field equations of motion allows us to then write the 
energy-momentum tensor in the convenient form 

V = i^{x)[d v + r"(aO]y> - \sl (sr - \g^R a a ) - g^\S 4 . (220) 

An incoherent averaging of i?/>7 M (x) [d u +T U (x)]ip over all the fermionic modes propagating in a Robertson- 
Walker background will bring the fermionic contribution to T^ v to the form of a kinematic perfect fluid 

TZ = -[(P m + Pm)UW+p m gn , (221) 
c 

with the conformal cosmology equation of motion then taking the form |5o"l ISol 13%] 

T" v = TZ - \% {bt - \g^R a ^j - g^XSt = . (222) 

Despite its appearance, the T^ v = condition can actually be satisfied non-trivially, since the van- 
ishing of the full T^ v can be effected by a cancellation of the perfect fluid contribution against the 
— (1/6)Sq (R^ — g^ u R a a /2) term associated with the back reaction of the scalar field on the geometry. 
Thus in the conformal theory there is energy density not just in the matter fields but in the geometry as 
well (conformal invariance requires the presence of the —S 2 R a a /12 term in Eq. (JBTj) with its necessarily 
negative sign), with it being this additional energy density which actually drives the theory. Finally, we 
note that despite the fact that it is only the full T^ v which is covariantly conserved, because the Einstein 
tensor and metric tensor terms in Eq. (J222|) are independently covariantly conserved, it follows that 

83 While the trace of the energy-momentum tensor must in general always vanish in the conformal theory to thereby 
already give us control of its various components, for cosmology the energy-momentum tensor itself must vanish too, to 
give us yet more control. 
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Xj^ is covariantly conserved also. In the S = So gauge then the perfect fluid does not exchange energy 
and momentum with the gravitational field, and thus obeys precisely the same covariant conservation 
condition as is used in the standard theory, with fermionic particles on average moving geodesically in 
the Robert son- Walker geometry. 

With a slight rewriting of Eq. (|222jl as 



we recognize the conformal cosmology Eq. (|223|) as being of the form of none other than the cosmological 
evolution equation of the standard theory, save only for the fact that the standard G has been replaced 
by an effective, dynamically induced one given by 



Q r 3 
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Thus as noted in (HH] , conformal cosmology is controlled by an effective gravitational coupling constant 
which is repulsive rather than attractive, and which becomes smaller the larger So might be. 85 If now 
we define conformal analogs of the standard ^ M {t) and fl\(t) via 

where A = XSq, then in a Robert son- Walker geometry Eq. (|223j) yields 

R 2 (t) + kc 2 = R 2 (t) (p M (t) + n A (t)) , ti M (t) + n A (t) + ti k (t) = i , 

q (t) = U 1 + ^)n M (t)-n A (t) (226) 

2 V Pm ) 

as the evolution equation of conformal cosmology, 86 an evolution equation which only departs from the 
standard Friedmann evolution equation of Eq. ()138|) through the replacement of G by G c g. sr With 
G e ff being smaller the larger So itself might be, we see that the larger the cosmological constant of the 
theory (viz. the larger XSq), the less it will contribute to cosmic evolution, with Cl A (t) having a self- 
quenching property not possessed by the standard fl\(t), as the latter comes with a fixed, pre-assigned 



The energy-momentum tensor of Eq. (|222ll provides us with an explicit example of a phenomenon we discussed 
earlier, namely that covariant conservation of the kinematic fluid energy- momentum tensor is achievable even if the full 
energy-momentum tensor differs from the purely kinematic one. 

85 In the conformal theory local non-relativistic solar system gravity is controlled by the parameter {3 in the metric of Eq. 
(I204|l . With the conformal coupling constant a g not participating in Eq. I|188|l in homogeneous geometries such as the 
cosmological one in which the Weyl tensor is zero, while participating in the static source f(r) of Eq. (|200|) in geometries 
such as the inhomogeneous Schwarzschild one in which the Weyl tensor is non-zero, G c g is completely decoupled from 
the local G = /3c 2 /M associated with a source of mass M. With the sign of the local G being fixed by the sign of a g , 
a negative effective global cosmological G c g is not in conflict with the existence of a positive local G. The fact that the 
dynamically induced G c g is negative in the conformal theory had been thought of as being a disadvantage since it seemed 
to imply that the local G would be given by the same negative G e g and then be repulsive too. However, as we see, a 
repulsive global cosmological G c g and an attractive local G can coexist in one and the same theory, an aspect of the 
theory which can now actually be regarded as a plus given the recent discovery of cosmic repulsion. (Recall that the best 
fit line of Eq. (|141(l to the accelerating universe data would go through the negative £lM(to) = —0.34 if A = 0.) 

86 In Eq. I|226(l the quantity f2fc(t) is the same as the previously defined fife(i) = —kc 2 /R 2 . 

87 With conformal invariance forcing the matter field action to be second order, in geometries in which the Weyl tensor 
vanishes the gravitational equations are only second order, and thus coincide in form with those of the standard second 
order Friedmann equations. Thus even for cosmology a preference for simplicity of the equations would not favor standard 
gravity over the conformal alternative. 
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G. Consequently, in the conformal theory a large unquenched A can still yield a small fl\(t), with it 
thus being possible to have a A as big as particle physics suggests, and yet, as we shall see in detail 
below, nonetheless not be in conflict with observation. The essence of the conformal gravity treatment 
of the cosmological constant then will be not to quench A at all but to instead quench the amount (viz. 
G e tt) by which it gravitates. 

In applying Eq. ()226j) to cosmology we have to recognize three distinct epochs. Based on our notion 
of a universe that goes through grandunification and electroweak symmetry breaking phase transitions 
as it expands and cools, and of each such phase transition having its own scalar field order parameter 
Sq and critical temperature Ty, the universe will start at the highest temperatures above all phase 
transitions (where every Sq will be zero), cool through a grand- unified phase transition temperature 
Ty UT with order parameter S"^ 111 ", and then subsequently go through the electroweak phase transition 
temperature Ty w with order parameter Sq w at which fermion and intermediate vector boson masses 
are generated. For temperatures above Ty w all particles will be massless, and TZ will act as a radiation 
fluid with 

Pm = 3p m = rj^r-r . (227) 

For temperatures below Ty W there will be a net induced cosmological constant which we can represent 
by an equivalent black-body with necessarily negative A = —aT v /c (the free energy having been being 
lowered by the phase transition, not raised) where Ty represents some blended value of Ty UT and Ty W . 
Analogously, Sq VT and Sq W can be represented by some blended value S . Since both of the Ty VT 
and Ty W temperatures are so huge, it will not matter which particular blended Ty we might use since 
compared to current temperatures it will be huge also; and with aT v then being overwhelmingly larger 
then the currently measured p m (to), it will not matter if we ignore particle masses altogether. 88 Thus 
we can define the model as being one whose evolution at temperatures less than Ty obeys Eq. (j226J) 
with p m = 3p m = aT^ and A = —aT v /c (Ty being huge), and which at temperatures above Ty obeys 

TZ = (228) 

with Sq being zero. 89 

At first glance the condition TZ — would appear to possess no non-trivial solution. However, 
that is not the case, since TZ is constructed here as an incoherent averaging of all the matter field 
modes propagating not in a flat space (where the solution to TZ — would of course be trivial) but 
in a curved one, with the equation TZ = being found [86 to actually support a non-trivial solution 
if the spatial 3-curvature of the Robertson- Walker geometry is expressly negative, since there is then 



88 For Ty of order say 10 15o K and p m (to) of order the detected luminous matter density in the universe at the current 
time to, the ratio cA/p m (io) will be of order 10 60 , with p m (t) only being competitive with cA at early universe temperatures 
of order Ty where particle masses are anyway unimportant. 

89 We should perhaps clarify what we mean by So = here. Our view of the field So is that it is a classical order parameter 
field which represents a spontaneously broken vacuum expectation value of a fermion composite (typically quadrilinear for 
grandunification and mass generating bilinear for electroweak) which is only operative in the critical regions below Ty UT 
and Ty W , and thus absent altogether at temperatures above Ty T . However, there could still be fundamental scalar 
fields in the theory as well (as there indeed would for instance be if conformal gravity is given a conformal supergravity 
extension). While such scalar fields would have vanishing vacuum expectation values above the highest critical temperature 
(the vacuum then being normal rather than spontaneously broken), nonetheless in such a vacuum the expectation value 
of S 2 would not vanish, since a quantum scalar field S can still connect a normal vacuum \ to the one-particle state \n) 
created out of it. In such a case the energy- momentum tensor associated with the kinetic energy of the scalar field would be 
given (c.f. Eq. §Q) as T^ v = (2/3)S^S' v -(l/6)g^S' a S- a -(l/3)SS^' u +(l/3)g^SS' a . a -(l/6)S 2 (W v - (l/2)g^R a a ) 
as evaluated in modes which obey S'^„ + (l/6)SR p, fl = 0, with the classical S now denoting the quantum scalar field 
matrix element (Q\S\n). Then even though the vacuum is such that (J7|5|f2) is zero, an incoherent sum over all |n) would 
still generate a non-zero contribution to the scalar field bilinear T^ v which would then be included as part of T£? n . 
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negative energy density present in the gravitational field to effect the needed cancellation. Conformal 
cosmology thus specifies the sign of A;, 90 leading to the same negative value that we previously found 
would impact on galactic rotation curves. With conformal cosmology early universe dynamics thus 
fixing the sign of k (something not the case with the standard Friedmann cosmology), and with k being 
a general coordinate scalar, the global topology of the universe would then not change as the universe 
cools down, with k then still being negative at temperatures below Ty where Eq. ()226|) applies. 



10.2 Conformal gravity and the accelerating universe 

Given the above considerations, the signs of various parameters of the theory are then fully specified, 
with k, A and G e s all being negative and p m {t) being positive. Given such a pattern of signs, the tem- 
perature evolution of the cosmology is then completely determined. Specifically, with the deceleration 
parameter q(t) being given by 

q(t) = n M (t) - n A (t) (229) 

when p m = pm/3, we see that q(t) is automatically negative in every epoch, with the universe thus being 
a permanently accelerating one no matter what the magnitudes of the various parameters. 91 Given the 
various signs of the parameters, conformal cosmology is found to admit of the exact solution 



where 



2 k{p - 1) fc/feinhV/'cQ 
R (t) = (230) 

(When A = this solution reduces to the k < 0, a > solution given in Eq. (jl33|) .) With Eq. (|230|) 
entailing that the initial time R(t = 0) is finite, conformal cosmology is singularity-free (G e s being 
repulsive rather than attractive). The cosmology thus expands from a finite rather than zero minimum 
size -R min = — k{(3 — l)/2a, and thus from a finite maximum temperature T max ~ l/i? mm , with the 
temperature evolution being given by 

2/3sinh 2 (aV2 ct) a (T max + Ty) 

I -'max 1 V! 

tanhV^ct) = [1-^—1 (1 + IL^SS. J . (232) 

\ I maxJ \ I V ) 

In the solution of Eq. (j230j) the temperature evolution of the quantities ^a(£) and fijvf(i) can be 
written in closed form as 



T 2 

max 


= 1 + 


\a l ' 2 ct) 






n A (f) = l 1 "^) ( 1 + ^ff E ) =( 1 -^) tanhV /2 ct) 



Cl M (t) = ~n A (t) , (233) 

with Eq. ()233)1 holding at any T(t) without any approximation at all. Quite remarkably, from Eq. 
()233|) we see that no matter what the actual value of T max , at temperatures which are well below it, the 
late time fl\(t 3> 0) will be given by 

ft A (t >0) =tanh 2 (a 1/2 ct) , (234) 
90 That k would be negative in the conformal theory was first noted in |85| . 

91 With G c g and A being negative and p m (t) being positive, the quantity ^a(£) i s then positive while the quantity 
O m (t) is negative. 
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to thus have to exclusively lie between zero and one, ultimately asymptoting to a bound of one from 
below. Thus no matter how big Ty might be, i.e. no matter how big A might be, in the conformal theory 
the current era 0,^(t ) is automatically bounded from above. With 0, M (t) being completely negligible 
in the current era (^m(*q)/^a(*o) = ~T£/T V = O(10~ 60 )), it follows from Eq. that the current 

era ftk(to) is given by 



with curvature thus contributing to current era cosmic expansion, 92 with the late time deceleration 
parameter q(t ^> 0) being given by the accelerating 



a deceleration parameter which is therefore automatically bounded between minus one and zero. 

Not only does conformal cosmology automatically quench the current era fl\(to) and thereby show 
that it actually is possible to live with a huge cosmological constant, since the cosmology has no initial 
singularity, it also solves the Friedmann early universe fine-tuning flatness problem, with the absence 
of any initial singularity not obliging the initial f2jvf(t = 0) + Cl\(t = 0) to have to be one [HH]. Further, 
the conformal cosmology also solves the cosmic coincidence problem, since CIm^o) is nowhere near close 
in value to fU^o) a ^ a ^ an d with Oat^o) not needing to be of order one (rather it is of order 1CT 60 ) 
the cosmology even has no need for any cosmological dark matter. With the conformal cosmology 
also having been shown to have no horizon problem or universe age problem jHHl EE] > it thus directly 
addresses many of the major challenges that contemporary cosmology has had to face, and does so 
without needing to appeal to an early universe inflationary era. 93 

With R(t) being completely determined in Eq. f)230|) . it is then straightforward to calculate the 
dependence of the luminosity function on redshift, with the late time luminosity function being found 
|3*%1 EH] to be given by none other than the pure A universe formula given earlier as Eq. (jl40|) . viz. 



as expressed in terms of the current era H(t ) and q . With q being treated as a free parameter 
which has to lie between minus one and zero in the conformal case, as described earlier the accelerating 
universe data are then fitted extremely well with q = —0.37 in the fits exhibited in Figs. 2 and 
3, with no fine-tuning of parameters being required. The success of these fits shows that having a 
particle physics sized cosmological constant is not incompatible with the supernovae data at all, it 
is only having a particle physics sized cosmological constant coupled to standard gravity which is 
incompatible. Since the conformal theory was not constructed after the fact in order to explain the 
accelerating universe data, and since it then does fit the accelerating universe data as naturally and as 
readily as it does, it should be given serious consideration. And as we indicated earlier, a. q = —0.37, 
^a(^o) = 0.37, f^M^o) = O(10~ 60 ), flkito) = 0.63 theory is a fully falsifiable one whose continuing 
acceleration predictions for higher redshift exhibited in Fig. 3 can be directly tested. 

Beyond our interest here in conformal gravity as a classical theory which is to be used to fit as- 
trophysical and cosmological data without dark matter or dark energy, the theory is also of interest 
as a quantum theory. Since it is based on an action I\y with a dimensionless coupling constant, as a 

92 With k being negative and 0^(2 ^> 0) being negligible, the effective late time sum rule CtA(t ^ 0) + Ofe(f ^> 0) = 1 
requires that &&(t 3> 0) has to lie below one, with it being the negative 3-curvature of the universe which was fixed in 
the very early universe thus forcing the current era Oa^o) to be bounded from above. 

93 Since a de Sitter geometry is conformal to flat, conformal gravity also admits of a de Sitter based cosmology |fi8| . 
with it not being excluded that there could have been an early universe inflationary de Sitter phase in the conformal 
theory even if the flatness and horizon problems that it addresses can be solved without it. 
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quantum gravity theory the conformal theory is power counting renormalizable, to thus have far better 
behavior in the ultraviolet that the non-renormalizable standard theory, a difficulty which one tries to 
get round in the standard theory by generalizing it to a 10- dimensional string theory. Because it is 
based on fourth order equations of motion, the unitarity of quantum conformal gravity has been called 
into question since it is suspected that the theory might possess on-shell negative norm ghost states. 
However, while this is indeed the case for the theories which involve both second and fourth order com- 
ponents, it has recently been shown jHZ] not to be the case for pure fourth order theories themselves, 
a point we elaborate on in Appendix D. Conformal gravity may thus well be a completely consistent 
quantum gravitational theory, one which is formulated purely in four spacetime dimensions alone. 

10.3 Quenching the contribution of the cosmological constant to cosmic 
evolution 

The essence of the conformal gravity approach to the cosmological constant problem is to have a 
symmetry which forces A to be zero in the exact symmetry limit (viz. above all cosmological phase 
transitions), and to then, as per particle physics, have a huge one induced as the temperature drops 
below the various critical points, but to have the huge induced one gravitate far less than it would 
do in the standard theory. Thus instead of trying to quench A, one can instead try to quench G. As 
such, a quenching of the cosmological G can be thought of as being a generic way to approach the 
cosmological constant problem, one which could apply in theories other than just conformal gravity, 
theories which need not depart so far from the Einstein-Hilbert action as the conformal theory does. 
Thus in the standard theory it could be the case that rather than be fixed, G could be a running 
coupling constant whose magnitude would depend on epoch or on distance, with a difference between 
the cosmological early universe G and the cosmological late universe G possibly allowing one to solve the 
early universe fine-tuning problem, and a difference between the large distance cosmological G and the 
short distance local G possibly allowing one to solve the cosmological constant problem. 94 Also, it might 
be possible to generate an epoch-dependent G through an embedding in a higher dimensional space, 
since as noted in Eq. (|161J1 . the 4-dimensional G is then not fundamental but is instead induced by the 
embedding. In having an epoch-dependent G which could quench the contribution of the cosmological 
constant to cosmology, such a quenching would ordinarily be expected to quench the contribution of 
ordinary matter to cosmology as well. However, in order to solve the cosmological constant problem 
per se, it is actually only necessary to quench the amount by which A itself gravitates, with there being 
no need to modify the amount by which p m (t) gravitates; with models which can do this (by taking 
advantage of the fact that unlike p m (t) the constant A carries zero 4- momentum) currently being under 
consideration Ed EIJ • Given the fact that attempts to quench the cosmological constant itself have 
so far foundered, attempting to instead quench the amount by which it gravitates appears to be an 
attractive alternative. 



11 Future prospects and challenges 

In this review we have examined standard gravity and some of its alternatives, and have identified the 
extension of the standard Newton-Einstein gravitational theory beyond its solar system origins to be the 
root cause of the dark matter and dark energy problems. However, attempting to critique the standard 
theory is not an easy enterprise, in part because the definition of what constitutes the standard theory 
is something of a moving target. While it is agreed upon that the standard theory is to be based on 
the Einstein equations, such a definition leaves a great deal undetermined, as it leaves T^ v completely 
unspecified, and does not make clear whether or not the gravitational side of the Einstein equations 

94 For a possible role for a running G in solving the dark matter problem see [55] , 
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is to also include a fundamental cosmological constant term (to cancel the one induced in by phase 
transitions in T^). 95 During the period between the development of the inflationary universe model in 
1981 and the discovery of the accelerating universe in 1998 the standard model was taken to be one in 
which both a fundamental and an induced cosmological constant were zero and the matter density was 
at the critical density (viz. Qm(0 = 1, ^a(^) = 0). Then, following the discovery of the accelerating 
universe this definition was revised to mean a universe with ^m(^o) — 0.3, and a net ^a(^o) — 0.7. Such 
a universe has a ready test since it predicts deceleration at redshifts above one (as per Fig. 3). However, 
if acceleration were to be established from analysis of the redshift greater than one Hubble plot, the 
response would not necessarily be to declare the standard theory wrong but to instead introduce, say, 
a quintessence fluid (or possibly even two or more such fluids) with a parameter w = p/p which would 
then depend on redshift in whatever way was needed. The weakness of the standard theory is not so 
much that one necessarily would do this, but rather that one could do it, i.e. that even now the theory 
is still not uniquely specified in a way that could enable it to make fully falsifiable predictions from 
which the theory could not subsequently back away. 96 

Nonetheless, even if the road which led to the present Qm^o) = 0.3, Oa(^o) = 0.7 standard paradigm 
was circuitous, it is still the case that these particular numbers work extremely well when the standard 
theory is applied to the anisotropy structure of the cosmic microwave background. In fact so well do 
they work that it is doubted that there could possibly be any second theory which could do as well. 
However, while it is extremely unlikely for lightning to strike in the same place twice, as far as the 
standard theory is concerned one has to ask if lightning has even struck once - namely do ^m(^o) — 0.3 
and ^a(^o) = 0.7 really describe the energy content of the universe, or is the true Pm(^o) in the real 
world the one given by luminous matter alone and the true A the one given by a particle physics scale, 
since if they are, the standard model fits would then be disastrous. 

However, if in the real world Pm(^o) and A are such that £)m (*o) is equal to 0.01 and ft\(to) is of order 
10 60 , one has to ask why does ^m(^o) = 0.3, ^a(^o) — 0.7 then work. As far as the supernovae Hubble 
plot data are concerned an answer can be provided to this question. Specifically, as we noted earlier, 
the best (f2 M (t ), ^a(^o) plane fits to the supernovae data lie along the fl\(to) = I-I^m(^o) + 0.37 curve 
given as Eq. (|141|) . The best fit minimum is thus quite shallow, and if the conformal gravity theory is 
correct and the real world lies close to the Om(^o) = 0, 0,\(to) = 0.37 solution to Eq. ()226|) (a solution 
which as far as purely phenomenological fitting is concerned is equivalent to using the Qm^o) — 0, 
^a(^o) — 0.37 solution to Eq. ()138j0 . the Q M (t ) = 0.3, Q\(t ) = 0.7 standard model solution would 
then just happen to be close to the real solution. Of course, the argument cuts both ways, and it could 
be that it is the conformal theory fit which is the one which is accidentally correct. To resolve this issue 
will require fitting the cosmic microwave background with conformal gravity to see if it works too, 97 
and if it does, to then see whether the success of either of the two theories (or of any other candidate 
theory) could reveal why the successes of the other theories were then only accidental. 

95 The freedom standard gravity has in choosing T^ v is due to its not using the traceless SU(3) x SU(2) x U(l) based 
T pu provided by fundamental physics. And we note that if we were to couple standard cosmology to a traceless T^, for 
k — we would obtain R(t) ~ t 1 / 2 and cLl = cz/H(to), with a best fit to the 54 data points of Fig. 2 yielding \ 2 = 138. 

96 The contrast with conformal gravity theory for instance in this regard is that the conformal theory is a uniquely 
defined theory (with gravitational action Iw and SU(3) x SU(2) x U(l) invariant matter action) which unambiguously 
predicts continuing acceleration above a redshift of one. Showing that the universe is in fact decelerating above a redshift 
of one would thus rule out the conformal theory, while showing that it was accelerating above a redshift of one would 
not necessarily rule out the standard theory. Alternate theories such as conformal gravity and MOND which endeavor to 
describe the universe without dark matter or dark energy do not enjoy the luxury of still having functions to adjust, and 
any major failure would rule them out. 

97 To apply the conformal theory to the cosmic microwave background requires the development of a theory for the 
growth of inhomogeneities in the model, to see what size standard yardstick the fluctuations would impose on the 
recombination sky, and what size they would then appear to us. Since the size of the fluctuation yardstick of the 
conformal theory would be different from the one used in the standard theory, the conformal gravity value Qk (to) = 0.63 
need not be in conflict with the standard model's fife (to) = 0. 
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The challenge to alternate theories to fit the cosmic microwave background data is a daunting one, 
not just calculationally, but also because alternate programs are extremely manpower limited with few 
people working on them. 98 However, this has to be weighed against the fact that the cosmological 
constant problem is an extremely challenging one for the standard theory, and here there are many 
workers and yet little progress." However, while this is all for the future, for the present three instructive 
guidelines have been identified - namely that there is good evidence that it is a universal acceleration 
scale which determines when a luminous Newtonian expectation is to fail to fit data, that a good case 
can be made that there is a global cosmological effect on local galactic motions which galactic dark 
matter might be simulating, and that in attempting to solve the cosmological constant problem one 
does not actually need to quench the cosmological constant itself, but only the amount by which it 
gravitates. 
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A The potential of a thin disk 

In order to determine the weak gravity potential of an extended object such as a disk of stars each 
with gravitational potential V(r) = —/3c 2 /r + / yc 2 r/2, we follow an approach originally developed 
by Casertano jHS] fc> r both thin and thick Newtonian disks and then generalize it to disks with linear 
potentials. For the Newtonian potential of an initially non-thin axially symmetric distribution of matter 
sources with matter volume density function p(R', z') we need to evaluate the quantity 



r 2ir roo P' n( P' ?'\ 

Va(R,z) = -(3c 2 dR! d(j)' dz'— — V ? ^77? ( A1 ) 

PK ' ' H Jo Jo r J-oc {R 2 + R' 2 -2RR'cos<j)' + {z-z') 2 y/ 2 v ; 

where R', <p', z' are cylindrical coordinates of the source and R and z are the only observation point 
coordinates of relevance. To evaluate Eq. (jAl|) it is convenient to make use of the cylindrical coordinate 
Green's function Bessel function expansion 

/ dkJ^k^J^kr'y™^-^-^-*' 1 , (A2) 



T — T JO 
\f I \ 771= — OO 



with its insertion into Eq. (|Aljl yielding 

poo roo roo 

Vb(R,z) = -2nf3c 2 / dk dR' dz'R'p(R',z')J (kR)J (kR')e- klz ~ zl . (A3) 

JO JO J-oo 

For the case of infinitesimally thin disks with p(R',z') = Y,(R')5(z') Eq. (|A3|) then simplifies to 

Vb(R) = -2nl3c 2 / dk / dR!R'T J (R')J (kR)J Q (kR') (A4) 
Jo Jo 

for observation points in the z = plane of the disk. Hence, for a disk with an exponential matter 
distribution S(-R') = Xoe" ^ and total number of stars = 27rSo/o; 2 (Rq = 1/a is the scale length of 
the disk), use of the standard Bessel function integral formulae 

, dR'R'J (kR')e- aR = , (A5) 

o (or + k A yi l 



UkR) R_ 
(a 2 + P) 3 / 2 2a 



k (^) Kl (-)_,(-),„(-) 
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(A6) 



then leads directly to 
Vp{R) 



-2npc 2 Z / dk 



aJo(kR) 
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-7r/3c 2 S i? 
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Finally, on differentiating Eq. ()A7|) and using the modified Bessel function relations 

I (z)=h(z) , I[(z) = I (z) - , K' Q {z) = -Ki(z) , K[{z) = -K (z) - ^ 



we obtain the relation 
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N(3c 2 R 2 a 3 r 



(Ra 
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presented in the text. 

The utility of the above formalism is that it immediately generalizes to the linear potential case 
(and by extension to the cubic potential and so on), since on setting \r — r'\ 
immediately obtain the potential 



r — r') 2 /\f — r'\, we 
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VJR,z) = / dR' d<f)' dz'R'p{R',z')[R 2 + R' 2 -2R&cos<f>' + {z-z') 2 ] l/2 
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(R 2 + R' 2 + (z- z') 2 )J {kR)J {kR' 



-2RR'J 1 (kR)J 1 (kR') 
For an infinitesimally thin disk Eq. (jAlOj) then reduces at z = to 



-k\z-z'\ 



VJR) = 7r 7 c 2 / dk / dR'R'Z{R') (R 2 + R' 2 )J {kR)J {kR') - 2RR'J 1 {kR)J 1 {kR') 



The use of Eqs. (jA5|) and (jA6|) and the additional integral formula 
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then yields 



2 r°° ( aR 2 J Q (kR) 9aJ (kR) 15a 3 J (kR) QakR.h(kR)\ 
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Through use of the modified Bessel function relations given above taken in conjunction with Eq. (|A6 
and its derivatives, Eq. (jA13|) readily evaluates to 

tt7c 2 E 0j R [ ( Ra\ ( 
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(A14) 



Finally, differentiation with respect to R and repeated use of the recurrence relations of Eq. (jA8|) then 
yields the expression 

(A15) 



RV(R) = (fO Kl (%) 



V 2 J 



presented in the text. 
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B The potential of a separable thick disk 



For non-thin disks simplification of Eqs. (jAf j) and (|A10|1 can be achieved for disks with a separable 
matter distribution p(R',z') = T,(R')f(z') where the symmetric thickness function f(z') = f(-z') is 
normalized according to 

/oo roc 
dz'f(z') = 2 / dz'f(z') = 1 . (A16) 
-oo JO 

Recalling that 

e -k\z-z'\ = _ z 'y-Kz-z') + g ( z , _ z y+Hz-z>) ? ( A17 ) 
we find that Eqs. (|A1|) and ()A10|) then respectively yield for points with z = 

/•oo />oo roo 

V {R) = -4vr/?c 2 / dk / dR' / dz'R'Z(R')f(z')J (kR)J (kR')e- kz (A18) 
Jo Jo Jo 



and 
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Jo Jo Jo 

x \(R 2 + R' 2 + z' 2 )J (kR)J (kR') - 2RR'J 1 (kR)J l (kR') 



(A19) 



in the separable case. Further simplification is possible if S(-R') is again the exponential Soe aR \ with 
use of the recurrence relation J[(z) = Jq{z) — J\{z)jz then yielding 
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RVl(R) = 2Npc 2 a 3 R / dk dz 

Jo Jo 



00 j, f°° Jizy-^'kj^kR) 



{a 2 + k 2 fl 2 



(A20) 



and 
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(A21) 

Further simplification is possible if a specific form for f(z') is specified, with two commonly consid- 
ered ones being of the form 
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with both of them falling off very rapidly once z is much greater than the scale height Zq. The thickness 
function of Eq. (|A22|) is found to lead to rotational velocities of the form 
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and 
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Similarly, the thickness function of Eq. ()A23|) leads to 
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In all of these expressions the needed functions of and its derivatives all converge very rapidly 
to their asymptotic values as their arguments increase. Consequently the relevant k integrations all 
converge very rapidly. As a practical matter, galactic scale heights z$ are usually much smaller than 
galactic disk scale lengths Rq. Consequently, the thickness corrections are only of significance in the 
inner galactic region, and thus have essentially no effect on the linear potential contribution. For the 
Newtonian term the thickness corrections of Eqs. (jA24j) and (jA27|) both tend to slightly reduce the 
overall Newtonian contribution, and serve to help ensure that the inner region rotation curves of Fig. 
1 are well described by the luminous Newtonian contribution alone. 



C The potential of a spherical bulge 

For a spherically symmetric matter distribution such as the central bulge region of a galaxy with radial 
matter density a(r) and N = 4tt J dr'r' 2 a(r') stars, the potential is readily found to take the form 



for a Newtonian potential, and take the form 



4ir(3c 2 



r Jo 
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(A30) 



for a linear potential. Both of these expressions can readily be integrated once a specific form for a(r) 
is specified, and it is from Eq. (jA29|) that dark matter halo contributions to galactic rotation curves are 
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calculated. Despite the simplicity of these expressions, because of projection effects, it is unfortunately 
not the 3-dimensional a(r) which is directly measured in spherical astronomical systems. Rather, it is 
only the two-dimensional surface matter distribution I(R) which is measured, with cr(r) having to be 
extracted from it via an Abel transform 



1 f°° I'(R) 



r (r 2 - R 2 ) 1 / 2 



(A31) 



to thus initially lead to double integrals in Eqs. (|A29|) and (|A30[) . 

Reduction of these integrals to one-dimensional ones which only involve the measured I(R) is however 
possible since on introducing the strip brightness S(x) which obeys 



\R 2 -x 2 y/ 2 



a(x) 



S'{x) 
2-kx 



we can rewrite Eqs. (IA29|) and ()A30|) as 



rVfcr) 



-2f3c 2 S{r) 



2(3c 2 



for the Newtonian potential and 
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for the linear one. Then, use of the relations 
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enables us to conveniently reexpress Eqs. (jA33|) and (jA34|) entirely in terms of I(R), to yield 

ARp 2 roc 

rV'^r) = — / dRRI(R) 
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for the Newtonian potential jHZI and 
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for the linear one. 
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D The ghost problem in fourth order theories 

With the fourth order propagator 

D(k 2 ,M 2 ) 



1 



(fcg - k 2 ){kl - k 2 - M 2 ) 
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(A38) 



associated with the prototype action 

i = -\J d A x (M 2 d^sd^s + dpdySdrars) (A39) 

and equation of motion 

(-<9 2 + V 2 )(-<9 2 + V 2 -M 2 ),S = (A40) 
being writable as a sum of two opposite signatured second order propagators, viz. 

D(k\M 2 ) = ^ , (A41) 

M 2 {k 2 -k 2 -M 2 ) M 2 {k 2 -k 2 ) 

it is quite widely thought that when quantized, fourth order theories such as conformal gravity would 
possess ghost states. However, the equation of motion which is given in Eq. (jA40|) is not actually 
the equation of motion of a pure fourth order theory per se, rather it is that of a second plus fourth 
order theory. With the pure fourth order action being the one obtained by setting M 2 to zero in Eq. 
(|A39|) . and with the separate positive and negative signatured second order propagators in Eq. (jA41|) 
becoming singular in this limit, no conclusion about the particle content of the pure fourth order theory 
can be immediately drawn. 

To investigate what the particle content of the fourth order theory does look like, it is sufficient to 
specialize to field configurations of the form S(x,t) = q(t)e lh ' x , configurations in which Eq. (jA40|) then 
reduces to 

^| + K 2 + c 2 )^|+^ 2 g = , (A42) 

where 

uf + oj 2 2 = 2k 2 + M 2 , ufuj = k 4 + k 2 M 2 , (A43) 

with Eq. (|A42|) itself reducing to the equal frequency u)\ = uo 2 when M 2 = 0. The equation of motion 
given in Eq. (jA42|) can be derived by variation of the acceleration-dependent action first introduced by 
Pais and Uhlenbeck [HE] 

Ipu= \S dt ~ ^ + u2 ^ 2 + u ^ q2 ] ( A44 ) 

where 7 is a constant, and we thus seek to quantize the theory based on Ipu to see what happens when 
we take the U\ — > uj 2 limit. 

For the theory associated with Ipu considered first as a classical theory, we cannot treat q and q as 
independent coordinates, since if q is to be an independent coordinate with canonical conjugate dL/dq, 
we could not then use dL/dq as the canonical conjugate of q. In order to reexpress this theory in terms 
of an unconstrained set of variables, we introduce a new variable x to replace q, with the method of 
Dirac constraints then converting the classical theory associated with Ipu into a theory of two coupled 
classical oscillators (q,p q ) and (x,p x ), with the unconstrained classical Hamiltonian 

2 

tt Px 1 ,7/2, 2\ 27222 / a Ar\ 

H= — +p q x+ -(wj +u 2 )x - -u x u 2 q (A45) 

then resulting [H7j. With this Hamiltonian generating the closed set of Poisson bracket relations 

{x,p x } = 1 , {q,p q } = 1 , 

{x, H } = — , {q, H} = x , {p x , H} = -p q - 7O 2 + u 2 )x , {p q , H} = ^\u 2 2 q (A46) 

7 

for Poisson brackets defined by 

dAdB _dA dB OA OB _ OA OB 

dx dp x dp x dx dq dp q dp q dq 
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this Hamiltonian is indeed the correct classical one, and thus the appropriate one for quantization. 
At the classical level, Eq. (jA46|) yields equations of motion of the form 



fx ■ ■ I 2 , 2\ ■ 2 2 / a ao\ 

x = — , q = x , p x = -p q - 7(0^ + lo 2 )x , p q = 7^i^2? > ( A48 ) 

7 

to nicely recover Eq. (|A42|) . with the substitution of this solution into the Hamiltonian of Eq. (|A45|) 
yielding a stationary Hamiltonian 

= \? - \H + - - 7 g^| (A49) 

whose time independence in solutions to Eq. (|A42|) can readily be confirmed. The unequal frequency 
theory has explicit two-oscillator solution 

q(t) = a x e~™ xt + a 2 e~ iuJ2t + c.c. (A50) 

with energy 

i/sTAT^i 7^ ^2) = - ujl){a\axujl - a* 2 a 2 u>l) , (A51) 

while the uj\ = uo 2 = uo equal frequency theory has explicit solution 

q(t) = ae-™ + c 2 te- luJt + c.c. (A52) 

with energy 

-^stat(^i = ^2) = 47a; 2 (2C2C2 + iuic\c 2 — iuc 2 Ci) . (A53) 

As we see, despite the fact that the equal frequency solution involves a temporal runaway, the appro- 
priately defined energy is still time independent, with temporal runaways not being the problem for 
fourth order theories they are often thought to be. 

Before proceeding to quantize the fourth order theory, we note that a comparison of Eqs. (jA51|) 
and (jA53|) reveals a difference between the unequal and equal frequency theories already at the classical 
level, one that will prove to be of great significance in the following. Specifically, while the unequal 
frequency theory is diagonal in the (ai, a 2 ) basis, with one of the oscillators having negative energy (the 
disease which translates into negative norm ghost states in the quantum theory), the equal frequency 
theory is not diagonal in the (ci, C2) basis. Moreover, if we set C2 = we obtain 

#stat(wi = uj 2 ) = , (A54) 

while if we set c\ = we obtain 

H stat (uj 1 = uj 2 ) = 8^yuj 2 c 2 c 2 . (A55) 

Of these two modes, we see that only the c 2 one carries energy. The zero energy result of Eq. (jA54|) is 
reminiscent of the conformal gravity zero energy theorem, with Boulware, Horowitz and Strominger |55] 
having shown that at the classical level the fourth order gravity theory energy would vanish identically 
if the gravitational field solutions were required to be asymptotically flat. With the restriction to 
asymptotic flatness in the gravitational case being equivalent in the Pais-Uhlenbeck case to requiring 
no temporal runaway, the equivalence of Eq. (jA54j) to the zero energy theorem of jD9j is manifest. 
However, in our analysis of the conformal gravity analog of the Schwarzschild solution, we found that 
conformal gravity also possessed the non-asymptotically flat linear potential term of Eq. (J2U4)) . We now 
see from Eq. (jA55|) that when the asymptotic flatness requirement is dropped, the temporal runaway 
solution would then be allowed, and would (for an appropriate choice of the sign of 7) actually have a 
perfectly acceptable positive energy, with the constraints of the zero energy theorem (a theorem long 
considered to be a shortcoming of conformal gravity) then being evaded. With this is mind, we shall 
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now show that the quantization of the equal frequency Pais-Uhlenbeck theory will lead to one oscillator 
which propagates with positive energy and to a second oscillator which does not propagate at all, with 
the would be ghost mode not being an eigenstate of the Hamiltonian. 

To quantize the theory based on the Hamiltonian of Eq. (|A45J) . the introduction of the four Fock 
space operators a\, a\, a 2 and a 2 defined via 

q(t) = ai e- luJlt + a 2 e~ iuJ2t + H.c. , p q (t) = i-yu^^e-^ 1 * + i-fto 2 to 2 a 2 e- iuJ2t + H.c. , 
x(t) = -iuxCLie-™ 1 * - ito 2 a 2 e~ iuJ2t + H.c. , p x (t) = --fu 2 a ie - iuJlt - ju^e'^ + H.c. 

(A56) 

then furnishes us with a Fock space representation of the quantum-mechanical commutation relations 

[x,Px] = [q,Pq] = i , [x,q] = [x,p q ] = [q,p x ] = [p x ,p q ] = (A57) 
provided the Fock space operators obey 

[ai,aj] = - j\ ' [ a 2>4] = o j— 2 2vT ' [ a i>4] = , [a 1 ,a 2 }=0 . (A58) 

In terms of the Fock space operators the quantum-mechanical Hamiltonian takes the form 

H = 1^{uj\ — uj\){uj\a[ai — u 2 a 2 a 2 ) H — (u>x + u 2 ) (A59) 

2 

with its associated commutators as inferred from Eq. (jA46|) then automatically being satisfied. With 
Ux, uj 2 and 7(0;^— all being taken to be positive for definitiveness, we see that the [ai, a\] commutator 
is positive definite while the [a 2 , a 2 ] commutator is negative definite. And with the Hamiltonian being 
diagonal in the a\ai, a 2 a 2 occupation number operator basis, we see that the state defined by 

ai|fi) = , a 2 \tt) = (A60) 

is its ground state, that the one-particle states 

I + 1) = [2 7 ^K 2 - ^ 2 2 )] 1 / 2 a t 1 |fi) , I - 1) = [2 7 ^ 2 K 2 - u 2 2 )] 1/2 a j 2 \Q) (A61) 

are both positive energy eigenstates with respective energies u)\ and uj 2 above the ground state, that the 
state I + 1) has a norm equal to plus one, but that the state | — 1) has norm minus one, a ghost state. 
Thus, as anticipated from Eq. (jA41j) . we find that the Hamiltonian of the unequal frequency theory can 
indeed be diagonalized in a basis of positive and negative norm states. 100 In the presence of interactions 
that would be added on to the Pais-Uhlenbeck action Ipu-, the asymptotic states for the S-matrix that 
such interactions would then generate would be the eigenstates of the free Hamiltonian of Eq. (|A59|) . 
While there are negative norm asymptotic states in the eigenspectrum of the free Hamiltonian, in and of 
itself that would only be a disaster if interactions would connect in and out states of opposite signature, 
and it is thus of interest to note that Hawking and Hertog |l()()j have argued via a path integral treatment 
of the constraints of Ipu that this does not in fact occur, with the unequal frequency theory then being 
viable. Hence even when M 2 is non-zero, assessing fourth order theories on the basis of the structure 
of Eq. (jA41|) is too hasty, and a careful quantization of the theory which takes constraints into account 
first needs to be made. 



100 Despite the presence of ghost states, the propagator of the theory is still causal, with the relative minus sign in Eq. 
(|A4ip not affecting the fact that each of the two second order propagators which appear in D(k 2 , M 2 ) is a standard causal 
second order propagator. 



84 



With the commutation relations of Eq. (jA58|) being singular in the it)\ — > uj 2 limit, we see that we 
cannot infer anything about the structure of the equal frequency Fock space from a study of the unequal 
frequency one. To take the limit we need to find some new set of Fock operators whose commutation 
relations would instead be non-singular. To this end we thus introduce the basis 

1 / , 2bu\ 1 / , 2bu\ 
ai = -\a-b + —— I , a 2 = -\a-b —I , 

a = a 1 (l + ± ; ) + a a (l-±) , b=^-(a 1 -a 2 ) . (A62) 
for the unequal frequency theory where we have set 

u = , e = . (A63) 

These new variables are found to obey commutation relations of the form 

[a, a f ] = A , [a, b f ] = \i , [b, a f ] = \i , [b, b f ] = v , [a, b] = , (A64) 

where 

e 2 (2u 2 — e 2 ) 
X = v = 7 2 2T~3 ' ^ = TpTT^i r^l, ■ ( A65 ) 

In terms of these new variables the coordinate q(t) of Eq. f)A56|) gets rewritten as 



q(t) = e 



-iujt 



. s 2ibu . 
ojcos et sm et 



+ H.c. , (A66) 



and thus has a well defined e — >• limit, viz. 

q(t, e = 0) = e~ iujt (a -b- 2ibut) + H.c. . (A67) 
Similarly, in the new variables the Hamiltonian of Eq. (|A59|) gets rewritten as 

(26 + 6 + a% + 6 f a) + u , (A68) 

and it too has a well-defined limit, viz. 

H(e = 0) = 87w 4 (26 t 6 + a ] b + 6 f a) + 00 , (A69) 

with the following commutators of interest having limiting form 

[H(e = 0),at] =a;(a t + 26 t ) , [#(e = 0), a] = -u(a + 2b) , 
[H(e = 0), 6 f ] = ub ] , [F(e = 0), b] = ~ub , 

[a + 6, a f + 6 f ] = 2/t , [a - b, a f - 6 f ] = -2/i , [a + b, a ] - b ] ] = , (A70) 

where p, = /i(e = 0) = l/(87u; 3 ). 

With the equal frequency theory now being well-defined, we find that for the Fock vacuum 
defined by a\Q) = b\Q) = 0, H(e = 0)\Q) = u\Q), the H(e = 0) Hamiltonian possesses only one and not 
two one-particle states. The state b'\Q) is an eigenstate with expressly positive energy 2u, while the 
state a/\Q) is not an eigenstate at all, a quantum-mechanical one-particle spectrum which thus reflects 
the features of the classical spectrum exhibited in Eqs. (jA54|) and (|A55|) . With this same pattern 
repeating for the multi-particle states, the equal frequency theory is thus a very unusual one in which 
the full Fock space has the same number of basis states as a two-dimensional harmonic oscillator, while 
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the Hamiltonian itself only has the number of basis states associated with a one-dimensional one. The 
(v\Q) states thus only couple off shell (where they can serve to regulate the ultraviolet behavior of 
theory) but do not materialize as on-shell asymptotic states. While the disappearance of these modes 
from the eigenspectrum of the Hamiltonian solves the ghost problem in the theory, this disappearance 
is still quite perplexing and requires further explanation. 

The reason for this highly unusual outcome derives from the fact that while the normal situation for 
square matrices is that the number of independent eigenvectors of a square matrix is the same as the 
dimensionality of the matrix, there are certain matrices, known as defective matrices, for which this is 
not in fact the case. A typical example of such a defective matrix is the non-Hermitian two-dimensional 
matrix 



with c non-zero but otherwise arbitrary, since even though this matrix has two eigenvalues both of 
which are real (despite the lack of hermiticity) and equal to 1 (no matter what the value of c), solving 
the equation 



leads only to q = when c 7^ 0, and thus to only one eigenvector despite the two- fold degeneracy of the 
eigenvalue, with the space on which the matrix M acts not being complete. The defective matrix M 
given above is in the form of a Jordan block matrix, with Jordan having shown that under a similarity 
transform an arbitrary square matrix can be brought to either a diagonal form, or to a triangular form 
such as that of Eq. (|A71|) in which every matrix element on one side of the diagonal is zero. Jordan 
block form matrices have the property that no matter what values the matrix elements take on the other 
side of the diagonal, the secular equation for the eigenvalues only involves the elements on the diagonal 
itself. Then if these diagonal elements are all real, the matrix will have real eigenvalues despite not being 
Hermitian. (While Hermitian matrices must have real eigenvalues, there is no converse theorem which 
would oblige the eigenvalues of non-Hermitian matrices to necessarily be complex, with non-Hermitian 
matrices still being able to possess real eigenvalues in certain cases.) As such, a Jordan block form 
matrix can be thought of as being a diagonal matrix M\ to which has been added a second matrix M 2 
with non-zero elements on only one side of the diagonal. With this second matrix being a divisor of 
zero, while its addition to the diagonal matrix does not affect eigenvalues, its being a divisor of zero does 
affect eigenvectors, to thereby cause the full Jordan block matrix M\ + Mi to have a smaller number of 
eigenvectors than eigenvalues. 

To see how these considerations apply in the case of interest to us here, we need to track the e — > 
limit carefully. As regards first the unequal frequency Hamiltonian of Eq. (|A68|) . we note that its action 
on the one-particle states a^Q), yields 




(A71) 




(A72) 



h<J\q,) 



1 r 



(4lu 2 + e 2 )a j \{l) + (4lu 2 - e 2 p\Q) 



Htf\n) 



2cu L 

1 r 



eV|fi) + (Alu 2 - e 2 )tf\tt) 



(A73) 



2lu 



In this sector we can define a matrix 




(A74) 



whose eigenvalues are given as 2lu + e and 2lu — e. For such eigenvalues, energy eigenvectors which obey 




(A75) 
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are then readily constructed as 



\2uj ± e) = [±ea f + (2uj =f e)6 f ] |f2) . (A76) 

As we see, as long as e 7^ 0, the two one-particle sector eigenvectors of the Hamiltonian H are distinct, 
and H has two eigenvectors to go with its two eigenvalues. However when we now let e go to zero, 
the two eigenvectors in Eq. (jA76|) collapse onto a single eigenvector, viz. the vector the two 

eigenvalues collapse onto a common eigenvalue, viz. 2uj, and the matrix M(e = 0) of Eq. (lA"74|) 
becomes defective, with the H(e = 0) Hamiltonian then being free of negative norm eigenstates. The 
equal frequency theory is thus fully acceptable, to thereby resolve what had been thought of as being 
one of the biggest difficulties for fourth order theories. 

As regards the conformal gravity theory itself, we note that in a g^ v = rj^ + h^ v linearization of the 
theory around flat spacetime, the conformal gravity rank two tensor of Eq. (jl85)l reduces to 



r = -w p ii urj K pa - -nrurKp,, (A77) 



where 

K u» = h K» _ - v ^h a a , = rf v d a d a - d»d u . (A78) 

In the conformal gauge 

d v <r - \g^guA9 up = o (A79) 

(viz. the gauge condition which is left invariant under g lxv {x) — > e 2a( - x ^g lJjU (x)) the source- free region 
gravitational fluctuation equation W^ u = then reduces to 

(-<9 2 + V 2 ) 2 K» U = . (A80) 

With this equation of motion being decoupled in its tensor indices, we see that each tensor component 
precisely obeys none other than the M 2 = limit of Eq. (|A40[) . Consequently, it would be of interest 
to see if the structure we have found for the prototype equal frequency Pais-Uhlenbeck theory can carry 
over to the conformal gravity theory once its gauge and tensor structure is taken into consideration, since 
that would then permit the construction of a fully renormalizable, fully unitary gravitational theory in 
four spacetime dimensions, one which despite its fourth order equation of motion, would nonetheless 
only possess one on-shell graviton and not two. 101 



101 It would also be of interest to see whether the typical linear in time growth given in Eq. (|A52(I might enable conformal 
gravity early universe tensor fluctuations to grow to macroscopic size by recombination. 
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